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EXECUTIVE  SUMMARY 


A  number  of  empirical  (statistical,  regression-oriented)  and  mechanistic 
(process-oriented)  models  are  presently  available  to  examine  the  relationship  between 
air  pollution  stress  and  plant  response. 

However,  both  types  of  models  have  their  strengths  and  weaknesses.  The  experi- 
mental and  modelling  steps  in  an  impact  assessment  program  are  closely  related;  the  data 
collected  often  dictate  the  type  of  model  that  can  be  developed  or  applied,  while  speci- 
fying the  use  of  a  given  model  at  the  outset  will  influence  the  experimentation.  In  air 
pollution  research,  it  is  common  to  find  these  two  steps  treated  separately  and  even 
conducted  by  different  researchers,  a  situation  that  invariably  leads  to  confusion  and 
conflict  between  the  biological  validity  of  a  model  and  its  mathematical  representation. 

Short-term,  acute  exposures  enable  the  development  of  models,  represented  by  a 
response  plane,  while  long-term,  chronic  exposures  enable  the  development  of  models  rep- 
resented by  a  surface  spanning  several  response  planes.  Each  response  plane  represents 
the  potential  impact  of  air  pollution  at  a  defined  plant  development  stage  and  response 
planes  will  differ  in  shape  because  physiologically,  for  example,  a  crop  will  show 
differential  response  to  a  set  or  similar  exposure  at  different  periods  of  its  life. 
Much  of  the  published  literature  on  air  pollutant  impacts  on  agriculture  or  forestry  has 
not  been  able  to  account  for  this  differential  response,  with  scientists  simplifying 
biology  by  developing  a  whole  season,  dose-response  statistic.  In  regression  analyses 
using  published  data  supplied  by  the  original  researchers,  at  least  one  study  has  shown 
that  potentially  more  meaningful  impact  is  assessed  using  models  that  recognize  the 
relative  importance  of  growth  stages.  Many  air  pollutants  occur  in  a  stochastic  manner 
and  it  is  questionable  whether  whole-season,  dose-response  functions  that  are  so  commonly 
used  have  any  interpretive  value.  The  policy  makers  and  regulatory  personnel  prefer  a 
simple  approach  which  would  facilitate  implementation  and  administration  of  ambient  air 
quality  standards.  Mathematical  artifacts  introduced  through  the  computation  of  seasonal 
dose  or  through  long-term  averaging  techniques  do  not  make  much  biological  sense.  Perhaps 
an  approach  which  could  be  meaningful  consists  of  the  use  of  median  and  percentiles 
computed  from  short  duration  pollutant  concentrations  and  used  as  input  parameters  in 
models  which  describe  chronic  effects.  Such  an  approach  is  free  of  the  influence  of  the 
non-normal  distribution  of  the  occurrences  of  pollutant  concentrations.  However,  such  a 
strategy  would  require  the  re-computation  of  much  of  the  data  on  pollutant  exposure 
kinetics  and  plant  response. 

Ecological  effects  assessed  using  either  of  the  two  groups  of  models  defined  in 
this  report  will  have  a  degree  of  uncertainty  associated  with  the  quality  of  data  used 
in  the  modelling.  Empirical-regression  models  are  mainly  plot-level  models,  while 
explanatory-simulation  models  are  mainly  process-level  models.  In  the  strictest  sense, 
a  model  is  only  applicable  at  the  level  from  which  it  is  derived,  yet  we  are  often  faced 
with  making  impact  assessments  at  the  regional  level.  This  possibility  of  incurring  a 
scale  error  is  one  of  the  most  critical  issues  facing  air  pollution  researchers  interes- 
ted in  impact  assessment.  There  is  no  general  agreement  among  researchers  on  how  to 
deal  with  the  scale  problem,  and  while  this  situation  persists,  any  policy  formulated  on 
regional  impact  assessment  must  acknowledge  the  uncertainty. 
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Recommendations 

Based  on  what  is  discussed  in  this  report  regarding  air  pollutant  exposure  and 
vegetation  response,  we  recommend  the  following: 


1.  The  way  to  obtain  a  simple,  single  (or  few)  equation(s)  "summary  model", 
possibly  for  regional  scale  applications,  is  to  start  with  a  "comprehen- 
sive" computer  simulation  of  the  biology  (all  the  relevant  components  and 
processes,  as  they  are  known  to  function).  Presumably,  this  model  would 
be  subjected  to  field  testing,  improved  if  needed  and  re-tested.  Then,  it 
is  used  to  derive  the  "summary  model"  equation(s),  from  the  behaviour  and 
sensitivity  analysis  of  variables  in  the  comprehensive  model. 

2.  Summary  models  should  be  structured  in  terms  of  the  probabilities  of 
responses  derived  from  probability  density  functions  of  the  input  parame- 
ters and  forcing  functions. 

3.  Due  to  the  lack  of  decisive  evidence  demonstrating  direct  vegetation 
effects  induced  by  acid  aerosols  or  acidic  precipitation  under  field 
conditions,  no  models  should  be  considered  at  this  time  for  these  parame- 
ters per  se.  However,  they  should  be  included  as  input  parameters  in 
examining  other  pollutant  effects  -  in  other  words,  examine  joint  effects. 

4.  Present  pollutant  averaging  techniques  for  describing  exposure  kinetics 
make  neither  biological  nor  mathematical  sense.  Parameters  such  as 
frequency  of  exposure,  peak  episodal  concentrations,  and  numerical  area 
under  individual  exposures  should  be  used  as  model  inputs  in  defining 
pollutant  exposure  kinetics  (see  Sections  2  and  3). 

5.  For  the  purposes  of  the  present  study  (ADRP),  "computer  simulation"  is 
required  prior  to  the  design  of  field  experiments,  data  collection,  and 
analysis  of  air  pollutant  exposure  and  vegetation  response. 

6.  For  gaseous  pollutant  effects  on  plants,  the  models  of  Coughenour  (1981), 
SILVA  (Kercher  and  Axelrod  1981),  and  SUCROS  (van  Keulen  and  de  Wit  1982) 
are  good  starting  points.  Summaries  of  these  models  are  described  in  the 
Appendices  of  this  report.  While  these  models  are  mechanistic  in  nature, 
they  can  be  adapted  into  empirical,  regression  models. 
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1 .  INTRODUCTION 

A  model  can  be  described  as  information,  data,  principles,  and  the  like, 
arranged  or  grouped,  usually  mathematically,  so  as  to  represent  or  describe  a  certain 
idea  or  a  condition.  In  evaluating  the  impacts  of  air  pollutants  on  cultivated  and 
native  vegetation,  a  number  of  investigators  have  proposed  empirical  (statistical 
regression  oriented)  models  (Larsen  and  Heck  1976;  Benson  et  al.  1982;  Heck  et  al.  1982; 
Loehman  and  Wilkinson  1983;  Medeiros  et  al.  1983;  Nosal  1983,  1984;  and  Fox  et  al.  1986). 
Others  have  developed  mechanistic  (process  oriented)  models  (Ares  1979;  Andersson  et  al. 
1980;  Haines  and  Waide  1980;  Kercher  1980;  Luxmoore  1980;  West  et  al.  1980;  Coughenour 
1981;  Heasley  et  al.  1981;  Kercher  and  Axelrod  1981;,  Miller  et  al.  1982;  Harwell  and 
Weinstein  1983;  King  et  al.  1983;  and  Mortensen  1984).  These  and  other  similar  models 
have  attempted  to  explain  the  relationships  between  acute  or  chronic  pollutant  exposures 
and  vegetation  response. 

In  the  present  context,  "acute  exposure"  is  defined  as  the  occurrence  of  short- 
term  (few  hours  to  few  days),  comparatively  high  pollutant  concentrations.  On  the  other 
hand  "chronic  exposure"  is  defined  as  the  occurrence  of  long-term  (weeks,  months,  entire 
life  cycle  of  vegetation),  comparatively  low  pollutant  concentrations,  with  periodic, 
intermittent  episodes  or  peaks.  However,  there  are  exceptions  to  this  definition  of 
"chronic  exposure".  For  example,  in  parts  of  the  United  Kingdom  the  geometric  standard 
deviation  of  sulphur  dioxide  concentrations  in  rural  areas  is  fairly  constant  (roughly 
2.0).  This  allows  the  prediction  of  the  duration  of  exposure  to  concentrations  within 
defined  thresholds,  from  the  air  quality  monitoring  data  (Fowler  and  Cape  1982). 

Generally  "acute  exposures"  cause  symptoms  of  injury  on  vegetation.  On  the 
other  hand  "chronic  exposures"  may  or  may  not  result  in  symptom  production.  Depending 
on  the  timing  of  the  stress,  whole  plants  can  frequently  recover  from  the  effects  induced 
by  an  acute  exposure  in  one  part  of  the  plant  (e.g.,  group  of  leaves  which  are  at  a 
sensitive  growth  stage),  through  overall  compensatory  processes  of  the  entire  plant.  On 
the  other  hand,  of  greater  concern  are  the  effects  which  result  from  long-term  chronic 
exposures.  Such  chronic  effects  consist  of  changes  in  growth,  productivity,  reproduction, 
and  quality.  These  effects  may  or  may  not  be  preceded  by  symptoms  of  injury.  In 
addition,  chronic  exposures  can  result  in  changes  in  features  relevant  to  aesthetic  and 
recreational  aspects. 

The  ability  to  understand  and  predict  temporal  and  spatial  patterns  of  air 
pollutant  effects  on  vegetation  continues  to  constitute  an  important  need  in  agricultural 
and  natural  resource  management.  Thus,  in  this  report  we  describe  and  evaluate  the 
state  of  published  mathematical  models  of  air  quality  and  vegetation  response.  We  also 
identify  the  limitations  in  these  models  and  describe  some  basic  biological  and  statis- 
tical concepts  of  plant  response  processes  useful  in  experimental  design  toward  a  new 
and  progressive  approach  to  the  problem  analysis. 
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2 .  CHARACTERISTICS  OF  AMBIENT  AIR  QUALITY 

The  pollutants  in  the  atmosphere  exist  as  gases,  vapours  or  particles  (solid 
and  liquid).  The  mechanisms  for  the  transfer  of  these  to  plant  surfaces  consist  of  "dry 
deposition"  (diffusion,  brownian  motion,  impaction  and  sedimentation)  and  "wet  deposi- 
tion" (such  as  rain,  snow,  and  fog).  Any  resultant  vegetation  effects  observed  can  be 
due  to  direct  stress  on  the  plant  species,  to  indirect  stress  mediated  by  the  soil,  or 
to  a  combination  of  both.  Torn  et  al.  (1987)  and  Mayo  (1987)  have  reviewed  the  published 
literature  on  the  effects  of  acidic  and  acidifying  pollutants  on  crop  and  tree  species, 
respectively.  In  evaluating  the  effects  of  pollutants  under  ambient  conditions,  one  must 
define  the  spatial  and  temporal  variability.  In  addition,  because  of  the  nature  of 
atmospheric  processes,  two  or  more  pollutants  can  exhibit  in  a  24  hour  period,  increases 
or  decreases  in  their  concentrations  simultaneously,  sequentially,  inversely,  or  in  some 
variable  or  poorly  defined  pattern  relative  to  each  other.  For  example,  in  a  typical 
situation  of  photochemistry,  high  concentrations  of  ozone  (O3)  and  nitric  oxide  (NO) 
occur  simultaneously,  while  high  concentrations  of  nitrogen  dioxide  (NO2)  precede  O3 
(NAS  1977).  While  high  concentrations  of  O3  are  generally  observed  during  daylight 
hours  (the  exception  being  stratospheric  intrusion  of  Oa),  high  concentrations  of  fine 
particle  sulphate  (SO4)  are  observed  during  nighttime  hours  (Stevens  et  al.  1978). 
Thus,  under  field  conditions,  vegetation  is  first  exposed  to  high  concentrations  of  fine 
particle  SO4  during  the  early  part  of  the  day  and  subsequently  to  high  concentrations 
of  Oa  during  the  late  afternoon  hours.  The  significance  of  this  has  been  demonstrated 
by  Herzfeld  (1982)  and  reviewed  by  Chevone  et  al.  (1986).  When  plants  were  exposed  to 
fine  particle  SO4  followed  by  Oa  under  controlled  conditions,  the  injurious  effects 
of  Oa  were  significantly  increased,  in  comparison  to  the  effects  induced  by  Oa 
alone.  The  aerosol  by  itself  did  not  produce  detectable  negative  effects.  Thus,  certain 
pollutants  such  as  Oa  and  sulphur  dioxide  (SO2)  are  known  to  be  relatively  more 
phytotoxic  than  others  and  under  appropriate  conditions  exist  at  sufficient  ambient 
concentrations.  Other  pollutants  such  as  NO2  and  SO4  are  known  to  be  relatively 
less  phototoxic  and  do  not  appear  to  exist  at  ambient  concentrations  sufficient  to  cause 
injury.  However,  what  is  critical  is  the  joint  effect  of  the  two  types  of  pollutants. 
In  addressing  this  issue  not  only  must  one  define  the  comparative  relationships  of  the 
occurrence  patterns  of  pollutants  of  concern,  but  also  provide  correctly  and  satisfac- 
torily an  appropriate  numerical  description  of  frequency  distributions  and  time  series 
analysis  of  individual  pollutant  occurrences. 

Many  scientists  assume  that  the  frequency  distribution  of  pollutant  concentra- 
tions follows  a  normal,  "bel 1 -shaped"  distribution.  The  following  is  an  example  of 
normal  distribution:  X  is  said  to  have  the  normal  distribution  if  the  distribution 
function  of  X  is  given  by 

^  1  -l/2[(y-y)/cj]^ 

F(x)  =  P(X<x)  =  J  — — — —  e  dy 
-00      /  if  -ir  o 

where  it  can  be  shown  that  the  parameters  y  and  a  are  the  mean  and  standard  deviation  of 
X.     The  standard  normal  distribution  is  the  normal  distribution  with  p  equal  to  0  and 
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cf  equal  to  1  (Snedecor  and  Cochran  1978).  Long-term  pollutant  averaging  techniques 
assume  this  (e.g.,  Oshima  et  al.  1976;  Heagle  et  al .  1986).  In  reality,  if  frequency 
distributions  of  the  occurrence  of  Oa  or  oxides  of  nitrogen  are  plotted  graphically 
(Figure  1),  their  distribution  is  skewed  with  a  long  tail  at  high  values.  This  is  also 
true  for  SO2  (Table  1).  According  to  Fowler  and  Cape  (1982),  if  the  frequency  classes 
are  plotted  in  terms  of  the  logarithm  of  SO2  concentrations,  however,  then  a  more 
symmetrical  distribution  is  obtained,  with  data  from  UK.  An  example  of  lognormal 
distribution  is: 


A  similar  conclusion  was  also  reached  by  Male  (1982)  in  the  US. 

The  assumption  of  log-normality  has  been  questioned  by  Berger  et  al.  (1982)  in 
Belgium.  These  authors  used  24-hour  averages  of  SO2  concentrations  over  a  two  and 
one-half  year  period  in  the  region  of  Ghent,  Belgium.  Lognormal  distribution  was  fitted 
to  the  data  with  poor  results.  The  differences  between  the  empirical  data  and  the  log- 
normal  curve  were  most  pronounced  at  extreme  values  of  the  95th  and  higher  percentiles. 
For  this  reason,  the  authors  used  the  two-parameter  gamma  distribution: 


The  gamma  distribution  provided  a  much  better  fit  and  the  goodness-of-f i t  was  confirmed 
by  several  tests  of  significance. 

Most  recently,  Buttazzoni  et  al.  (1986),  evaluating  the  sulphur  dioxide  concen- 
trations in  the  area  of  Venice,  Italy  concluded  that: 


1.  Their  results  seem  to  suggest  that  2,  3,  or  4-parameter  lognormal  distribu- 
tions do  not  hold  characteristics  of  generality  or  universality  for  the 
interpretation  of  air  quality  data; 

2.  In  fact,  Weibull  and  gamma  distributions  yield  superior  simulations  in 
Venice; 

3.  Bell-shaped  distributions  cannot  fit  L-shaped  experimental  data  sets;  and 

4.  Gamma  distributions  are  less  convenient  than  Weibull  distributions. 


According  to  the  authors,  as  a  final  comment,  "the  research  for  a  universal  form  allowing 
for  direct  comparisons  among  distributions  underlined  by  differing  sites  or  time  spans, 
should  head  toward  a  generalization  of  such  distributions". 


In  addition  to  the  aforementioned  studies,  Lefohn  and  Benedict  (1982)  and  Nosal 


(1984)  examined  the  frequency  distribution  patterns  of  ambient  O3  concentration  in 
parts  of  the  US.  They  concluded  that  such  distributions  follow  a  Weibull  function 
(Figures  1  and  2).    An  example  of  Weibull  distribution: 


(X>0) 


1 


F(x)  = 


(3'*r(a)  0 


Px^^)  =  {(x-5o)/a}  ^'"'e  -t(x-^o)/"'^''  (^Q<x) 
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Figure  1.  Frequency  distributions  for  ambient  ozone,  nitric  oxide,  and  nitrogen  dioxide.  In  the 
upper  left,  the  frequency  distributions  for  ozone  are  broken  down  by  years,  whereas  in 
the  bottom  left  they  are  broken  down  by  site.  Only  minor  differences  were  observed 
among  sites.  Among  years,  1978,  and  particularly  1980,  had  more  high  Oa  concentra- 
tions. The  upper  right  shows  the  frequency  distribution  for  nitric  oxide  broken  down 
by  site.  The  concentrations  were  similar  at  all  three  sites.  The  bottom  right  shows 
the  frequency  distribution  for  nitrogen  dioxide  broken  down  by  site.  Major  differences 
in  the  concentrations  of  NO2  occurred  at  the  different  sites,  with  a  gradient  extend- 
ing out  from  the  urban  area,    (from  Pratt  et  al.  1983). 
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Table  1.  Frequency  distribution  of  averaged  hourly  SO2  concentrations 
downwind  from  a  large  point  source.  Data  from  7  air  quality 
monitoring  sites  were  examined  in  developing  this  summary 
stati  Stic . 


S02 

Concentration  (ppb) 

Absol ute 
Frequency 

Relati  ve 
Frequency  (%) 

Cumul ati  ve 
Frequency  (%) 

0 

89.192 

94.4 

94.4 

5 

3.446 

3.6 

98.0 

10 

1  ,316 

1  .4 

99.5 

15 

260 

0.3 

99.8 

20 

142 

0.2 

99.9 

25 

37 

0.0 

99.9 

30 

24 

0.0 

100.0 

35 

8 

0.0 

100.0 

40 

9 

0.0 

100.0 

45 

3 

0.0 

100.0 

50 

4 

0.0 

100.0 

70 

2 

0.0 

100.0 

130 

1 

0.0 

100.0 

Source:    Fowler  and  Cape 

(1982) 
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DENSITY  OF  WEIBULL  DISTRIBUTION 


Figure  2.    Theoretical  distributions  illustrating  the  Weibull  family  of 
frequency  distributions. 
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In  addition  to  the  previous  discussions  on  gaseous  pollutants,  the  phenomenon 
of  "acidic  rain"  is  of  major  international,  ecological  concern.  Our  knowledge  of  the 
observed  or  potential  effects  of  "simulated  acidic  rain"  on  crops  and  tree  species  have 
been  summarized  by  Irving  (1983)  and  Jacobson  (1984).  As  far  as  we  know  there  are  no 
demonstrated  cases  of  direct  negative  effects  of  "acidic  rain"  on  vegetation  under 
ambient  conditions.  Almost  all  studies  relate  to  the  use  of  "simulated  rain".  This 
consists  of  treating  plant  species  with  a  "solution  of  constant  chemical  composition, 
applied  artificially  at  constant  or  varying  rates  and  amounts".  In  our  opinion  this  is 
irrelevant  to  the  real  world  where  the  chemistry  of  precipitation  varies  significantly 
within  and  between  individual  rain  events  (Pratt  and  Krupa  1985).  The  observed  temporal 
and  spatial  variability  of  rainfall  properties  is  regulated  by:  (a)  regional  versus 
local  meteorology;  (b)  regional  versus  local  distribution  of  pollutant  sources  and  their 
characteristics;  (c)  the  frequency  of  rainfall;  (d)  rainout  versus  washout;  and  (e)  the 
concentration  of  coarse  particles  in  the  air  prior  to  the  rain  event. 

The  literature  contains  numerous  publications  where  "average"  values  of 
constituents  in  precipitation  have  been  used  to  develop  the  "simulant"  in  plant  exposure 
studies.  As  with  the  gaseous  pollutants  discussed  previously,  the  frequency  distributions 
of  the  concentrations  of  ions  in  precipitation  also  do  not  follow  normal  distribution. 
Such  distributions  are  best  described  by  a  Weibull  function  (Nosal  and  Krupa  1986). 

In  describing  the  relationships  of  ions  in  precipitation,  frequently  linear, 
statistical  regression  methods  are  used.  Such  techniques  assume  that  frequency 
distribution  data  is  normally  distributed,  and  the  residuals  are:  (a)  unbiased  (mean 
zero);  (b)  homoscedasti c  (equal  dispersion);  (c)  normal  (Gaussian  dispersion);  and  (d) 
independent.  Recently,  Nosal  and  Krupa  (1986)  showed  that  the  aforementioned  assumptions 
are  violated  by  precipitation  chemistry  data  (Tables  2-5).  In  addition,  there  are  a 
number  of  other  concerns  introduced  through  artifacts  caused  by: 

1.  Sampling  methodology  (Coscio  et  al.  1982;  Electric  Power  Research  Institute 
1986) ; 

2.  Sample  collection  period  (event  or  daily  versus  weekly  or  monthly)  (de 
Pena  et  al.  1985;  Ridder  et  al.  1985;  Sisterson  et  al  .  1985;  Krupa, 
unpublished);  and 

3.  Methods  of  numerical  analysis  used  in  the  data  analysis. 

Krupa  (unpublished)  found  that: 

1.  There  are  significant  differences  in  the  chemistry  of  precipitation,  as 
derived  from  numerical  analysis,  between  warm  (April  15  -  October  14)  and 
cool  (October  15  -  April  14)  seasons.    This  is  to  be  expected. 

2.  The  concentrations  of  certain  ions  (e.g.,  SO4,  NH4,  Ca)  are  greater  in 
the  warm  season  compared  to  the  cool  season.    This  is  also  to  be  expected. 

3.  However,  the  cool  season,  rather  than  the  warm  season  chemistry  regulated 
the  results  of  data  analysis  on  an  annual  basis  at  many  geographic  loca- 
tions in  the  US.  In  other  words,  results  of  statistical  analysis  of  annual 
data  resembles  more  closely  the  cool  season,  rather  than  the  warm  season 
chemi  stry . 
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Table  2.  Summary  statistics  of  the  distributions  of  wetfa11  parameters.^ 


Location  Parameter       Coeff.  Skewness'^       Coeff.  Kurtosis 


Fernberg                     pH  4.14  0.97 

logio  SO4  4.13  2.01 

Marcel!                       pH  3.94  -0.05 

logio  SO4  2.20  1 .30 

Lamberton                    pH  -2.63  -0.77 

logio  SO4  2.03  2.36 

All  Locations  Pooled                pH  3.30  -3.47 

logio  SO4  4.07  2.03 


Data  for  the  three  Minnesota  sites  were  gathered  by  the  National 
Atmospheric  Deposition  Program  (NADP  1978-1984),  Natural  Resources 
Ecology  Laboratory,  Colorado  State  University,  Fort  Collins,  Colorado. 

Standardized  test  values,  at  95%  significance  level,  critical  value 
for  normalcy  is  +  1.96  for  skewness. 


10 


Table  3.    Distribution    of    some   major   wetfall    chemical    components  at 
Marcell,  Minnesota.^ 


Distribution       Ion      Chi-Square      Significance      K  -  Significance 


Gamma 

S04 

82. 

.551 

1  , 

.58L-13 

.1635 

2. 

,28E-6 

Lognormal 

S04 

19, 

.322 

.013 

.0434 

,999 

Weibull 

S04 

65. 

.115 

1  , 

.37E-10 

.1126 

3. 

.OlE-3 

Gamma 

NOa 

96. 

.000 

0, 

.0 

.1978 

0, 

.0 

Lognormal 

NO3 

13. 

.135 

.0688 

.0566 

.3863 

Weibul 1 

NO3 

60. 

.659 

3, 

.46E-10 

.1110 

3, 

.51E-3 

Gamma 

PH 

23, 

.941 

.1210 

.0688 

.1550 

Lognormal 

pH 

21  , 

.906 

.14624 

.0623 

.2446 

Weibull 

pH 

52, 

.12/ 

1 

.96E-5 

.1007 

8, 

.36E-3 

Data  were  gathered  by  the  National  Atmospheric  Deposition  Program 
(NADP  1978-1984) . 

Kolmogorov  -  Smirnov  distance. 
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Table  4.  Summary  statistics  of  the  multivariate  regression  to  predict 
the  relationship  of  log  SO4  with  other  chemical  constituents 
in  wetf al 1 

Regression  equation: 

log  SO4  =  1.44  f-  0.307  log  Ca  »-  0.04^^  log  Mg 

+  0.346  log  NH4  +  0.174  log  Na  -  0.163  pH 


Variable  Coefficient  St.  Dev.  T-ratio 


Constant  1.4395  0.0545  26.42 

log  Ca  0.3067  0.0370  8.30 

log  Mg  0.0425  0.0350  1 .21 

log  NH4  0.3460  0.0183  18.93 

log  Na  0.1743  0.0154  11 .31 

pH  -0.1633  0.0109  -14.93 


S  =  0.1497;  r""  =  77.5%;  R^  adjusted  for  degrees  of  freedom  =  77.3% 


Analysis  of  Variance 

d.f.                 SS  MS 

5     -          42.9881  8.5976 

558  12.5101  0.0224 
563  55.4981 


Regression 

Residual 

TOTAL 


Minnesota  data  from  the  National  Atmospheric  Deposition  Program 
(NADP) 
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Table  5.  Summary  statistics  of  the  tests  of  normality  of  the  residuals 
in  Table  4. 


Coefficient  of  skewness  =  0.690941 
Standard  value  -  6.6989 
Root  (skewness)  =  2.588 
Lower  crit.  (95%)  =  -0.179 
Upper  crit.  (95%)  =  0.179 

Distribution  of  the  residuals  is  significantly  skewed  to  the  right. 


Coefficient  of  kurtosis  =  7.57087 
Standard  value  =  22.1581 
Lower  crit.  (95%)  =  2.67 
Upper  crit.  (95%)  =  3.37 

Distribution  of  the  residuals  is  significantly  leptokurtic. 


Chi-squared  =  25.44    (12  degrees  of  freedom) 
Significance  of  departure  from  normality  =  0.012 
Distribution  of  the  residuals  is  significantly  non-normal. 


Durbin-Watson  statistic  DW  =  1.66123 

DW  <  DW  upper  (5%)  =  1.78;  reject  rho  =  0 

Residuals  are  autocorrelated  with  respect  to  time;  hence  dependent. 


Stratification  of  residuals  by  pH  range    -  4.8  -  5.1  -  6.2 


Test  of  Homoscedastici ty 


Statistic  Significance 


Cochran's  C  =  Max.  Var  /  Sum.  Var.         0.3704  .0000 

Bartlett  -  Box  F  13.901  .0000 

Bartlett  1.07755  .0000 

Hartley's  Max.  Var.  /  Min.  Var.  2.503  .0000 


Distribution  of  the  residuals  is  significantly  heteroscedasti c . 
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While  it  is  not  clearly  described  in  a  number  of  publications  on  the  effects  of 
simulated  acidic  rain  on  vegetation,  we  wonder  how  many  scientists  used  annual  statistics 
in  designing  their  studies  with  "simulated  rain".  After  all,  plants  achieve  active  or 
pollutant  sensitive  growth  stages  (e.g.,  soybean  flowering  and  pod  filling)  only  at 
certain  times  of  the  year. 

Of  particular  concern  is  the  definition  of  the  time  series  of  the  occurrences 
of  rainfall  of  sufficient  intensity  (rate  and  amount)  and  chemistry  of  concern,  its 
relationship  to  the  occurrence  of  other  forms  of  pollutants,  and  the  timing  of  such 
interactive  pollutant  stress  in  relation  to  the  growth  stage  and  sensitivity  of  the  plant 
species.  This  aspect  of  concern  is  discussed  in  greater  detail  in  the  next  section  of 
this  report. 

As  opposed  to  the  long-term  emphasis  of  research  on  the  effects  of  gaseous  air 
pollutants  on  plants,  and  the  recent  impetus  to  study  the  impacts  of  "acidic  rain", 
little  if  any  effort  on  a  comparative  basis  has  been  directed  to  examine  the  effects  of 
fine  particles  that  contain  a  major  fraction  of  SOa,  NO3,  etc.  (Stevens  et  al. 
1978).  Many  plant  scientists  have  difficulty  in  addressing  this  question.  Fine  parti- 
cles (<  1.5  y)  which  are  monodisperse  and  of  uniform  size  and  concentration  are 
extremely  difficult  and  complex  to  generate  on  a  continuous  basis.  In  addition  to  this 
problem,  in  our  opinion,  many  plant  scientists  have  not  fully  applied  state-of-the-art 
atmospheric  physics  and  chemistry  in  evaluating  air  pollutant  -  plant  response  relation- 
ships. Understanding  the  nature  of  the  cause  is  as  important  as  understanding  the  nature 
of  the  receptor  response.  In  our  opinion,  sound  integration  of  the  two  aspects  is 
critical . 
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3.  THE  CONCEPT  OF  POLLUTANT  DOSE 

"Dose"  is  one  of  the  most  frequently  used  terms  in  air  pollution  -  vegetation 
effects  literature.  "Dose"  can  be  defined  as  the  exact  amount  or  extent  of  a  given 
treatment  to  which  a  given  receptor  is  subjected,  at  one  time,  or  at  stated  intervals. 
The  numerical  expression  of  dose  is  an  area  of  much  controversy  at  this  time. 

Historically,  dose  has  been  expressed  as  ct  -  pollutant  concentration,  c, 
multiplied  with  duration  of  time,  t,  of  exposure  or  as  an  average  of  pollutant  concen- 
tration over  exposure  duration  (c/t).  More  recently,  it  has  also  been  expressed  as 
"integrated  exposure"  (IE),  as  determined  by  summing  the  products  of  concentration  times 
the  exposure  period,  when  the  pollutant  concentration  exceeded  a  set,  minimum  threshold 
(Lefohn  and  Benedict  1982) . 

In  the  numerical  definition  of  dose,  several  considerations  are  critical,  in 
regard  to  both  the  cause  (air  pollutant)  and  the  receptor  (vegetation)  response.  Air 
pollutant  averaging  techniques  assume  normal  distribution  of  ambient  pollutant  concen- 
trations. However,  as  previously  stated  the  distributions  of  ambient  ozone  (O3)  and 
sulphur  dioxide  (SO2)  concentrations  are  best  described  by  Weibull  and  gamma  functions 
(Berger  et  al.  1982;  Lefohn  and  Benedict  1982;  and  Nosal  1984).  A  consequence  of  these 
non-normal  distributions  on  the  application  of  pollutant  averaging  techniques,  for 
example,  is  illustrated  by  the  work  of  Hogsett  et  al.  (1985).  In  this  study  two  exposure 
regimes  were  used  to  compare  a  set  of  episodic  ozone  exposure  patterns  with  daily 
exposures.  An  episodic  ozone  pattern,  applied  to  alfalfa  over  a  133  day  growth  period, 
elicited  a  greater  growth  response  compared  to  a  daily  peak  exposure  pattern.  The  two 
profiles  had  similar  integrated  exposure  values,  but  different  7-h  season  means.  The 
7  -h  (0900-1559  h)  means  for  the  two  profiles  were  inversely  correlated  with  the  observed 
growth  response.  The  episodic  profile  had  a  smaller  seasonal  mean,  and  yet  caused  a 
greater  yield  reduction. 

Another  example  of  the  artifacts  induced  by  pollutant  averaging  techniques  per- 
tains to  sulphur  dioxide  concentrations  downwind  from  a  large  point  source  in  Minnesota 
(Krupa  and  Gardner,  in  preparation).  Ambient  ground  level  sulphur  dioxide  concentrations 
were  monitored  continuously  over  a  10  year  period  at  seven  locations  relative  to  the 
source  plume  dispersion.  When  hourly  pollutant  averages  were  examined,  SO2  concentra- 
tions were  reported  to  be  zero  during  approximately  90%  of  the  monitored  period  (Table 
1).  However,  when  the  continuous  data  were  examined,  within  the  90%  time  there  were 
numerous  instances  when  the  plants  were  exposed  to  SO2  concentrations  as  high  as  0.50 
ppm  for  five  minutes.  Even  though  in  this  particular  case  no  negative  effects  were 
observed,  it  is  important  to  note  that  it  is  well  accepted  that,  with  some  exceptions, 
chronic  effects  are  the  result  of  the  cumulative  impact  of  intermittent  episodes.  Such 
a  cumulative  impact  is  not  a  straightforward  additive  function,  because  the  final 
biological  response  is  a  product  of  the  stress  minus  the  ability  of  the  plants  to 
compensate  (Zahn  1970).  This  brings  forth  the  concept  of  frequency  of  occurrence  of 
pollutant  episodes  and  the  time  interval  between  such  episodes.  During  the  time  interval 
between  episodes,  gaseous  pollutants  can  and  do  exist  in  low  concentrations.  Exposure 
of  plants  to  such  low  pollutant  concentrations,  can  either  predispose  them  to  subsequent 
episodes  or  can  confer  tolerance  (Godzik  and  Krupa  1982).  These  phenomena  are  confounded 
by    the    influence    of    other    environmental    variables    and    the    biological    time  clock 
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(discussed  in  the  following  paragraph).  Thus,  in  our  opinion,  under  field  conditions 
plants  do  not  respond  to  a  single  pollutant  exposure  threshold  value  for  injury,  but  a 
range  of  such  values.  Therefore,  we  believe  that  any  single  "cut-off"  value  as  "thresh- 
old" makes  little  biological  sense.  It  is  a  simplistic  approach  to  our  understanding  of 
pollutant  exposure  -  plant  response  relationships. 

From  the  perspective  of  receptor  (vegetation)  response,  plant  physiologists  and 
plant  pathologists  have  clearly  shown  the  importance  of  plant  growth  stage  in  the  end 
response  of  vegetation  to  biotic  or  abiotic  stress  (Teng  and  Gaunt  1980;  Benson  et  al. 
1982).  This  introduces  the  concept  of  biological  time  clock  and  stress  (Figure  3).  For 
example,  defoliation  or  foliar  injury  to  soybean  by  hail  during  the  early  part  of  the 
plant's  growth  would  have  no  effect  on  the  final  yield.  On  the  other  hand,  such  a  stress 
during  flower  set  to  pod  filling  would  result  in  significant  yield  loss  (Lockwood  et  al. 
1977). 

In  this  context  the  influence  of  variables  which  regulate  the  rate  of  plant 
growth  or  its  sensitivity  to  the  pollutant  stress  is  of  great  importance.  Fundamentally 
the  genetic  characteristics  of  the  plant  species  in  question  are  a  major  factor.  A  number 
of  investigators  have  shown  that  within  a  given  crop  species,  under  similar  growing 
conditions,  different  cultivars  of  that  crop  respond  differently  to  similar  pollutant 
exposure  regimes  (Unsworth  and  Ormrod  1982).  Because  of  the  genetic,  variability  or 
heterogeneity,  natural  or  managed  populations  of  trees  also  exhibit  varied  behaviour 
among  individual  trees  within  the  population.  This  overall  influence  of  the  inherent 
genetic  properties  of  the  plant  species  is  further  confounded  by  the  extraneous 
influence  of  atmospheric,  edaphic,  and  additional  biological  variables. 

Some  atmospheric  factors  include  light,  temperature,  relative  humidity,  and 
carbon  dioxide.  Fdaphic  factors  such  as  moisture  and  nutrient  availability  concurrently 
regulate  the  overall  plant  growth  and  behaviour.  The  large  temporal  and  spatial  vari- 
ability of  all  these  factors  is  a  major  contributor  to  the  uncertainty  in  regional  scale 
crop  loss  estimates.  These  environmental  factors  not  only  regulate  the  rate  of  plant 
growth,  but  also  many  of  them  influence  the  stomatal  behaviour.  The  stomata  serve  as  the 
entrance  in  plants  to  atmospherically  derived  chemical  energy  (carbon).  They  also  serve 
as  the  predominant  pathway  for  the  entrance  of  gaseous  pollutants  into  the  leaf.  Majernik 
and  Mansfield  (1972)  have  shown  that  in  bean  plants  stomata  are  opened  and  closed  to 
varying  degrees  by  the  influence  of  the  relative  concentrations  of  carbon  dioxide  to 
SO2.  Similarly,  plants  subjected  to  temporary  drought  followed  by  the  availability  of 
sufficient  moisture  exhibit  an  increased  sensitivity  to  O3  exposures.  Under  ambient 
conditions  such  factors  become  even  more  complex  with  the  influence  of  pollutant 
mixtures . 

A  number  of  investigators  have  developed  crop  loss  models  in  evaluating  the 
effects  of  pathogens  and  pests  (Teng  and  Krupa  1980).  In  this  field  of  research,  many 
process-oriented  models  include  the  role  of  environmental  variables.  However,  the 
influence  of  air  pollutants  has  not  been  included  in  such  models.  Plant  response  to 
either  the  pathogens  or  pests  and  to  air  pollutants  is  influenced  by  interaction  of  these 
variables  in  a  number  of  complex  and  yet  to  be  fully  understood  ways  [refer  to  Unsworth 
and  Ormrod  (1982)].  It  is  beyond  the  scope  of  this  report  to  provide  a  detailed  discus- 
sion of  this  subject. 
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Figure  3.    Conceptualization  of  the  response  surface  of  yield  loss  to 
pollutant  stress  as  influenced  by  crop  development  stage. 
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In  the  context  of  air  pollutant  exposure  and  plant  response,  a  satisfactory 
numerical  expression  of  dose  should  at  least  consider,  for  example: 

1.  the  artifacts  of  pollutant  averaging  techniques; 

2.  the  episodicity  of  pollutant  occurrence  and  exposure; 

3.  the  time  intervals  between  episodes;  and 

4.  the    relationship    between    the    pollutant    stress    and   the   biological  time 
clock  (growth  stage) . 

According  to  Runeckles  and  Brown  (1986),  one  model  which  does  not  provide  a 
simple  dose  statistic,  but  which  minimizes  the  uncertainties  of  overall  dose-response  is 
that  of  Nosal  (1983) . 

The  model  consists  of  a  mixed,  multivariate,  polynomial,  Fourier  regression: 

R  -  I      I  a.(")  x"^  ^    I      I    b.^"^)    sin  x. 
i'l    n    '         '      U^    i  ' 

where: 

R      =     plant  response   (crop  yield  parameters;  annual  basal  area  increment  of 

trees,  and  so  forth); 
xi    =      number  of  pollutant  episodes; 
X2    =      cumulative  integral  of  the  pollutant  episodes; 
tl 

•f    c(t)  dt  =  integral     (exceeding  1)  of  the    concentration    and  the  duration  of 

t2  each  episode,  cumulative  through  the  duration  of  the  study  (ti  -  t2); 

c(t)=-      pollutant  concentration  at  time  t; 

X3    =      peak  pollutant  concentration,  max  c(t)  in  ppb; 

tl  <  t  <  ti;  and 

ai('^),  bi('^)  =  model  parameters. 

The  model  parameters  are  estimated  using  the  least  squares  approach: 

3  3 
min  SSQ  -  (R  -  Z    I  a.^"^     x^.    -      I    Z  b.^"^^  sino.  x.)^ 
i=l  n  i=l  m 

The  actual  degree  of  the  polynomial  (range  of  the  subscript  "n"  in  the  first  double  sum) 
and  the  number  of  sin  terms  (range  of  the  subscript  "m"  in  the  second  double  sum)  vary 
depending  on  the  plant  response  type  and  the  pollutant  type.  However,  in  general  the 
polynomial  degree  n  <  2  and  sin  degree  m  <  1  and  u)^  =  1  appear  to  be  appropriate. 

Thus,  the  model  accounts  for  the  number  of  pollutant  episodes,  highest  peak 
pollutant  concentration,  and  the  cumulative  numerical  integral  (concentration  over  time) 
of  the  exposure.  The  weakness  of  this  model,  however,  is  that  it  does  not  provide  a 
complete  numerical  explanation  of  the  importance  of  the  individual  episodes  or  its 
relationship  to  the  biological  time  clock.  However,  it  considers  the  concept  of  epi- 
sodicity and  fluctuations  in  ambient  air  quality  over  time  with  response,  rather  than 
attempting  to  define  a  single  dose  statistic. 
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While  the  previous  narrative  pertains  to  gaseous  air  pollutants,  significant 
attention  has  been  directed  in  recent  years  to  the  phenomenon  of  "acidic  rain"  (US 
National  Research  Council  1983).  A  number  of  investigators  have  used  linear  or  linear- 
izable  statistical  regression  techniques  in  defining  the  temporal  and  spatial  aspects  of 
rainfall  chemistry  to  which  vegetation  is  exposed.  Recently,  Nosal  and  Krupa  (1986) 
showed  that  the  treatment  of  data  by  linear  or  linearizable  regression  methods  violates 
the  assumptions  required  for  the  use  of  such  techniques  (namely,  normal  distribution  of 
data,  normal  distribution  of  the  residuals  in  the  regression,  homoscedasticity ,  equal 
dispersion,  unbias  or  mean  zero,  and  independence).  Such  data  also  fail  to  qualify  under 
"robust  regression".  These  considerations  are  very  critical,  since  almost  all  studies 
on  the  impact  of  acidic  rain  on  vegetation  relate  to  simulation  experiments  which  are 
based  on  the  assumptions  found  to  be  violated  by  Nosal  and  Krupa  (1986). 

This  overall  discussion  points  to  the  difficulties,  artifacts,  and  errors 
associated  with  the  numerical  description  of  dose.  Unfortunately,  in  the  real  world  the 
vegetation  is  subjected  to  stress  in  a  different  manner.  Any  description  of  dose  must 
make  biological  and  statistical  sense.  If  the  two  requirements  can  be  achieved  simul- 
taneously, such  an  output  would  be  a  desirable  solution  to  the  problem. 

One  aspect  of  dose  which  has  not  been  considered  is  the  question  of  "exposure 
dose"  versus  "effective  dose"  (Runeckles  1974).  Exposure  dose  may  be  defined  as  the 
pollutant  regime  to  which  the  plant  species  is  exposed  as  defined  by  appropriate  ambient 
or  atmospheric  measurements.  Effective  dose  may  be  defined  as  the  actual  quantity  and 
rate  of  the  pollutant  which  is  absorbed  by  the  plant.  It  is  the  "effective  dose"  that 
results  in  the  observed  or  measured  effect.  However,  since  the  numerical  definition  of 
"exposure  dose"  itself  is  an  area  of  much  controversy,  no  effort  is  made  here  to  examine 
the  details  of  "effective  dose".  In  addition,  at  present  we  are  not  aware  of  techniques 
to  accurately  measure  the  "effective  dose"  under  field  conditions.  Possible  approaches 
include  the  use  of  tracers  such  as  stable  or  radioactive  isotopes.  Such  methods,  however, 
are  highly  complex  and  are  limited  by  the  time  of  integration  required  to  measure  their 
changes  in  concentration  within  the  plant  tissue.  Another  approach  is  the  computation 
of  pollutant  absorbed  dose  (PAD)  as  suggested  by  Fowler  and  Cape  (1982).  PAD  =  mean 
pollutant  concentration  (yg/m^)  x  exposure  time  (s)  x  canopy  conductance  (m/s). 
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4.  DEFINITION  OF  THE  TYPES  OF  MODELS 

An  "empirical  model"  as  used  in  the  literature  on  the  subject  matter,  has 
usually  meant  a  mathematical  equation  that  has  been  derived  statistically  from  a  set  of 
data.  According  to  Woodmansee  (1974)  a  statistical  model  can  be  defined  as  "an  equation 
in  which  the  dependent  variable  is  estimated  by  combining  one  or  more  independent 
variables  each  of  which  is  multiplied  by  an  appropriate  coefficient,  in  various  additive 
combinations".  In  contrast,  a  "mechanistic,  process  model"  is  defined  by  Woodmansee  as 
"one  which  attempts  to  represent  a  biological  system  in  terms  of  basic,  well  defined 
laws  or  relationships". 

With  regard  to  plant  response,  we  frequently  use  simplified  mathematical  models 
of  a  few  variables,  usually  describing  how  much  of  something  exists,  rather  than  at  what 
rate  it  is  changing.  Usually  the  description  consists  of  one  equation,  rather  than  a 
set  of  equations,  as  an  empirical  model.  An  example  is  the  crop-weather  model  used  for 
regional  and  inter-regional  forecasting.  At  the  field  level,  we  believe  that  the  empiri- 
cal, statistical,  "regional  scale  variable-based"  models  can  suffer  from  a  lack  of 
realism  and  generality  (ability  to  be  applied  in  conditions  outside  of  those  for  which 
it  was  developed).  As  a  result  we  choose  to  emphasize  in  this  review  the  other  types  of 
growth-production  models  which  generally  use  a  series  of  equations  to  describe  rates  of 
changes  of  various  physical  and  biological  processes,  along  with  environmental  con- 
straints which  can  change  those  rates.    Kercher  and  King  (1985)  also  take  a  similar  view. 

A  useful  framework  for  viewing  various  plant  growth,  development  and  production 
models  is  one  which  is  ecological  in  nature  and  refers  to  the  resources  that  control  or 
regulate  rates  of  growth.  This  framework  consists  of  four  increasingly  complex  levels 
of  production  (de  Wit  and  Penning  de  Vries  1982;  Penning  de  Vries  1982a, b).  Van  Keulen 
and  de  Wit  (1982)  presented  the  same  framework  in  the  context  of  more  comprehensive 
farming  systems  analysis  and  land  use  management. 

Crops  process  models  tend  to  be  in  one  of  the  following  categories: 

Production  level  1:    Plant  growth  is  limited  only  by  light  and  heat  (weather); 
Production  level  2:    Plant  growth  is   limited  by  shortage  of  water  availability,  during 
part  of  the  time; 

Production  level  3:    Plant  growth  is   limited  by  nitrogen  shortage,   during  part  of  the 
time; 

Production  level  4:    Plant  growth  is  limited  by  shortage  of  phosphorus  and  other  needed 
mi  neral s . 

A  production  level  which  should  be  included  is  the  "plant  growth  limited  directly  by 
industrial  pollutants".  Superimposed  on  all  these  would  be  models  of  the  interactive 
effects  of  plant  pathogens  and  pests.  We  know  of  only  one  review  which  covers  the  topic 
of  models  for  studying  the  effects  of  toxic  substances  on  crop  growth  (Kickert  and 
Benenati,  In  Press).  In  addition,  it  is  not  clear  to  us  the  extent  to  which  nutrient 
fertilizer  -  plant  growth  models  include  growth  retarding  effects  caused  by  fertilizer 
overdose. 

Table  6  provides  a  summary  of  the  advantages  and  limitations  of  the  two  major 
classes  of  models  (empirical  and  mechanistic)  reviewed  in  this  report. 
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Table  6.    Advantages  and  limitations  of  the  types  of  models. 


Mechanistic  Process  Models 

Empi  rical 

Feature  or 
Characteri  stic 

Statistical 

Regression  Preliminary 
Models  Model 

Comprehensi  ve 
Model 

Summary 
Model 

Relative  accessabi 1 i ty 
of  input  requirements 

High  (+)i 

Medium  - 
Difficult  ( 

-) 

Amount  of  input  data 
required  to  use/enter 

Low  (+) 

Medium  - 
High  (-) 

"Simplicity"^ 

"High"  (+)  "Medium" 

.  "Low"  (-) 

"Medium" 

Time  required  for 
computer  to  run 

Low  (+) 

Medi  um^ 

Amount  of  device  memory 
requi  red 

Low  (+) 

Medium  - 
High3 

Precision  of  responses 
(output) 

High  (+)^ 

Low  - 
Medium  (-) 

Abi 1 i ty  to  estimate 
response  probabilities 

High  (+) 

Low  (-)5 

"Predictive  value" 

"Medium"            "Low"  (-) 

"Medi  um" 

"High" 

Real i  sm 

LOW    \  ) 

Medium  - 
n  1  g n  \^ } 

"Scientific  value" 

"Low"  (-)           "High"  (+) 

"High"  (f) 

"Low" 

Applicability  to  other 
locations  (generality) 

Low  (-) 

(case  history)    High  (+) 

Medium  - 

Typical  application 
(geographic  area  size) 

Regional 
or  local 

Local 

"Instructive  value" 

"Low"  (-)  "Medium" 

"Low"  (-) 

"High" 

^  (+)  is  an  Advantage;  (-)  is  a  Limitation. 

2  Quoted  entries  are  from  Penning  de  Vries  (1983). 

3  Not  very  relevant  if  using  80286  microcomputers,  or  better. 
If  no  extrapol i tion  outside  the  range  of  original  database. 

5  Most    of  these    models  are    deterministic,  but    could  be    modified  by  an 
experienced  modelling  specialist  to  run  in  a  probabilistic  mode. 
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5.  MATHEMATICAL  MODELS  FOR  CHARACTERIZING  PLANT  RESPONSE  TO  AIR  POLLUTANT  STRESS 

Plant  response  can  be  defined  as  an  observable  or  measurable  change  in  the  con- 
dition of  a  plant  or  a  vegetation  community  due  to  the  exposure  to  some  air  pol 1 utant( s) . 
As  previously  discussed,  such  responses  can  be  classified  as  "acute"  or  "chronic". 
Independent  of  these  two  types,  we  are  particularly  interested  in  the  final  plant 
response  that  can  be  interpreted  by  natural  resource  and  agricultural  economists  in 
cost-benefit  assessment.  Such  an  assessment  is  required  in  formulating  and  implementing 
policies  relevant  to  air  quality.  In  the  present  context  there  are  two  distinct  steps 
involved  in  developing  cost-benefit  analysis: 

1.  Development  and  application  of  assessive  and  predictive  models  which  will 
provide  information  on  plant  response  to  air  pollutant  stress;  and 

2.  Development  and  application  of  economics  models  which  can  be  applied  to 
the  information  derived  in  step  (1). 

This  report  is  restricted  to  a  discussion  of  step  (1),  air  pollutant  exposure  - 
plant  response  models.  lables  7A  and  7B  summarize  the  characteristics  of  selected 
empirical,  regression  models.  A  similar  analysis  is  provided  for  the  mechanistic,  process 
models  in  lables  8A  and  8B.  The  actual  mathematical  equations  relevant  to  models  having 
sufficient  documentation  are  presented  in  Appendices  8.1  and  8.2.  The  information 
contained  in  the  Appendices  can  be  applied  to  other  appropriate  situations  by  investi- 
gators other  than  the  original  authors. 

Ihis  review  identifies  the  advantages,  limitations,  applicability,  and  suggested 
modification  of  the  models  (lables  7A  and  7B;  8A  and  88).  To  some  extent  the  properties 
described  in  these  lables  are  inherent  to  either  class  (empirical  or  mechanistic)  of 
model s  (lable  6) . 

5.1  ACUIE  POLLUTANl  EXPOSURE  AND  PLANl  RESPONSE 

There  are  a  number  of  models  for  describing  acute  pollutant  exposure  and 
subsequent  plant  response  (O'Gara  1922;  Thomas  and  Hill  1935;  Zahn  1963;  van  Haut  and 
Stratmann  1970;  Larsen  and  Heck  1976;  and  Umbach  and  Davis  1986;). 

Of  the  recent  studies,  Larsen  and  Heck  (1976)  presented  a  statistical  model  of 
acute  response  caused  by  Os  and  SO2,  in  terms  of  visible  injury  symptoms  on  the  leaves  of 
a  number  of  crops.  Their  equation  (refer  to  the  Appendices)  is  unusual,  somewhat  awkward 
to  use  for  evaluating  plant  response,  and  is  based  solely  on  pollutant  concentration  and 
duration  under  optimal  plant  growth  conditions.  Umbach  and  Davis  (1986)  used  essentially 
the  same  model  in  their  studies. 

Benson  et  al.  (1982)  provided  a  statistical  model  (refer  to  the  Appendices)  for 
acute  response  of  alfalfa.  While  it  is  simple  to  use,  it  should  only  be  used  for  very 
preliminary  assessments,  especially  in  locations  outside  Minnesota. 

A  statistical  model  for  describing  N02-induced  acute  injury  was  developed  by 
Heck  and  lingey  (1979)  (refer  to  the  Appendices).  As  with  most  regression  models  used 
to  describe  the  acute  effects  of  air  pollutants,  it  does  not  consider  sub-optimal 
conditions,  and  therefore,  should  only  be  used  with  great  caution  for  preliminary  field 
assessments . 
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As  with  the  three  regression  models  described  previously,  in  this  report  we 
refer  to  only  three  process-oriented  models  for  acute  response.  Kercher  (1980)  described 
a  comprehensive  plant  physiological  model.  However,  the  model  requires  improvement  in 
describing  the  pollutant  exposure  kinetics.  At  this  time  it  can  be  considered  for 
research  purposes  only. 

Harwell  and  Weinstein  (1983)  presented  a  preliminary  paper  for  the  design  of  a 
forest  model  capable  of  displaying  acute  effects.  Attempts  to  obtain  updated  information 
from  either  of  the  authors  were  unsuccessful.  At  present  it  cannot  be  used  in  other 
assessment  programs. 

Luxmoore  (1980)  described  an  elaborate  set  of  process  models  at  the  Oak  Ridge 
National  Laboratory  (ORNL).  This  model  is  useful  for  evaluating  acute  effects  in  eastern 
deciduous  forests  relative  to  bioaccumulation  of  sulphur  and  heavy  metals.  Other  air 
pollution  research  programs,  especially  those  addressing  local  point  source  emissions  of 
these  pollutants,  should  investigate  the  accessibility  of  the  documentation  for  this  set 
of  models.    To  our  knowledge,  the  required  documentation  has  not  been  published. 


5.2  CHRONIC  POLLUTANT  EXPOSURE  AND  PLANT  RESPONSE 

Many  models  have  been  published  to  describe  plant  response  to  chronic  pollutant 
exposure.  Empirical,  regression  models  for  crops  include  those  of  Oshima  et  al.  (1976); 
Stevens  and  Hazelton  (1976);  Benson  et  al.  (1982);  Heck  et  al.  (1982,  1984a, b);  Loehman 
and  Wilkinson  (1983);  Medeiros  et  al.  (1983);  Nosal  (1983);  Rowe  and  Chestnut  (1985); 
and  Heagle  et  al.  (1986).  Appropriate  information  on  any  particular  model  can  be 
obtained  by  crop  type  from  Tables  7A  and  78. 

For  forests,  only  the  statistical  models  of  Nosal  (1984),  and  Fox  et  al.  (1986) 
were  found  to  be  sufficiently  described  that  they  could  be  considered  for  replication 
elsewhere.  Clearly,  the  Nosal  approach  provides  better  results,  probably  because  of  a 
more  realistic  definition  of  the  pollutant  exposure. 

In  the  last  few  years,  a  number  of  mechanistic,  process-oriented  and  ecosystem 
level  models  have  been  published  to  analyse  response  of  vegetation  to  chronic  exposure 
to  air  pollutants.  Some  of  these  relate  to  crop  responses  (Kercher  1980;  King  et  al. 
1983;  Mortensen  1984;  and  King  1986).  The  model  of  King  (1986)  is  rather  simple  and 
easy  to  understand,  but  is  in  its  early  stages  of  development.  It  has  a  few  limitations 
(refer  to  Tables  8A  and  88)  which  need  to  be  rectified  to  make  it  useful  for  assessment. 

There  are  also  a  number  of  process-oriented  models  developed  for  forests, 
rangelands,  and  desert  shrub  land  exposed  to  pollutant  stress.  These  include:  Andersson 
et  al.  (1980);  Haines  and  Waide  (1980);  Kercher  et  al.  (1980);  Luxmoore  (1980);  West  et 
al.  (1980);  Coughenour  (1981);  and  Heasley  et  al.  (1981).  Kercher  and  Axelrod  (1984) 
provided  a  relatively  realistic  and  especially  well  documented  treatment  that  should  be 
used  as  a  foundation  on  which  to  improve  our  knowledge  of  forest  assessment.  However, 
as  with  most  other  models,  it  requires  a  more  realistic  mathematical  definition  of 
pollutant  exposure  in  terms  of  frequency  (episodicity) ,  duration,  and  concentration  of 
the  pol lutant( s) . 

For  semi-arid  rangelands,  Coughenour's  (1981)  comprehensive  model  is  well 
designed  and  documented  (refer  to  the  Appendices). 
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6.  DISCUSSION  AND  RECOMMENDATIONS 

6.1  REGIONAL  SCALE  APPLICATIONS 

The  use  and  applications  of  models  at  the  regional  scale  can  be  separately  con- 
sidered for  photochemical  oxidants  and  for  wet  and  dry  acidic  or  acidifying  pollutants. 

However,  there  are  a  number  of  uncertainties  attached  to  the  assessment  of  air 
pollutant-induced  ecological  effects.  Uncertainty  is  an  accepted  phenomenon  when 
modelling  biological  systems.  In  the  present  context,  "this  uncertainty  arises  from  the 
stochastic  nature  of  many  biological  stimulus-response  relationships  and  is  a  reflection 
of  the  inherent  stochasticity  of  both  the  biological  system  itself  and  the  environment 
that  drives  that  system  .  .  .  (this)  makes  it  almost  impossible  to  predict  a  response 
that  is  deterministic,  i.e.,  without  variation  (Krupa  and  Teng  1982). 

Within  the  context  of  actual  investigation,  uncertainties  originate  from  errors 
in:  (a)  monitoring  the  air  quality;  (b)  experimental  and  numerical  definition  of 
pollutant  exposure  kinetics;  (c)  experimental  design  to  characterize  the  cause  and 
effects  relationships;  and  (d)  modelling  the  pollutant-response  relationships. 

Further,  any  transfer  of  results  from  unit  level  models  to  regional  level  leads 
to  "scaling  error".  At  present  there  is  no  general  agreement  among  researchers  on  how 
to  deal  with  the  scale  problem. 

6.1.1  Air  Quality  Models  for  Photochemical  Oxidants 

There  are  a  number  of  problems  that  must  be  considered  in  applying  air  quality 
models  at  the  regional  scale.  The  unavailability  of  air  quality  data  at  remote  rural 
sites  is  one  problem.  The  limitations  of  "averaging"  air  quality  data  (such  as  the 
popular  7-h  daily  average)  also  presents  a  problem  in  that,  while  it  is  simple,  as 
discussed  previously,  it  is  unrealistic.  An  alternative  to  simple  averaging  is  the 
three-parameter  exposure  kinetics  described  by  Nosal  (1983,  1984).  Given  the  sparsity 
of  rural  oxidant  data,  the  three  parameters  of  exposure  can  be  estimated  for  sites  in  a 
limited  geographic  region  using  the  "Krieging"  process  of  spatial  interpolation  (Lajaunie 
1984).  However,  extreme  caution  should  be  used  in  applying  this  process,  particularly 
if  the  number  of  actual  data  points  in  the  interpolation  is  very  limited. 

There  are  also  a  number  of  assumptions  made  in  the  use  of  "Krieging".  Such 
assumptions  need  to  be  tested  at  the  outset,  before  the  procedure  can  be  applied  to 
develop  geographic  predictions  of  photochemical  oxidants. 

6.1.2  Photochemical  Oxidants  and  Regional  Scale  Plant  Response  Modelling 

At  least  in  the  US,  many  believe  that  regional  scale  applications  and  air 
quality  standard  setting  require  a  rather  simple,  one-independent-variable,  single 
equation  model.  This  model  is  best  defined  on  the  basis  of  statistical  analysis  of 
field  data  collected  just  for  this  purpose.  A  number  of  regression  models  discussed  in 
Tables  7A  and  78  were  developed  with  this  view.  These  models  appear  to  provide  high 
values  of  the  coefficient  of  determination  (R^)  because  of  the  use  of  data  sets  of 
small  sample  size. 

A  group  of  counterpart  scientists  regard  the  mechanistic,  process-oriented,  or 
ecosystem  models  as  suitable  only  for  local  site  or  point  source  applications.    We  agree 
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with  this  if  one  were  thinking  only  of  comprehensive  models  as  described  by  Penning  de 
Vries  (1982b,  1983). 

A  number  of  scientists  outside  the  US,  and  not  in  the  field  of  air  pollution  - 
vegetation  effects  research,  have  argued  that  the  way  to  obtain  a  simple,  single  (or 
few)  equation(s)  "summary  model",  possibly  for  regional  scale  applications,  is  to  start 
with  a  "comprehensive"  computer  simulation  of  the  biology  (all  the  relevant  components 
and  processes,  as  they  are  known  to  function).  Presumably,  this  model  would  be  subjected 
to  field  testing,  improved  if  needed,  and  re-tested.  Then,  it  is  used  to  derive  the 
"summary  model"  equation(s),  from  the  behaviour  and  sensitivity  analysis  of  variables  in 
the  comprehensive  model  (Holling  1978;  Penning  de  Vries  19B2b,  1983).  This  is  the 
approach  that  we  recommend. 

In  addition  to  those  process  models  reviewed  in  Tables  8A  and  88,  there  are 
quite  a  number  of  preliminary  and  comprehensive  process  models  of  crops  and  tree  species 
not  yet  adapted  to  application  in  air  pollutant  -  vegetation  response  research  (Kickert 
1984a) . 

To  apply  this  strategy,  one  must  evaluate  the  results  of  sensitivity  analysis 
on  the  comprehensive  process  models  to  derive  one  or  more  summary  models.  The  process 
models  for  Oa  effects  (Tables  8A  and  88)  on  crops  and  tree  species  are  not  sufficiently 
complete  to  allow  their  use  in  the  approach  described  here.  As  a  result,  summary  models 
for  03  exposure  and  plant  response  might  be  derived  from: 

1.  Using  the  aforementioned  strategy  on  the  appropriate  comprehensive  models 
of  SO2  effects  (Tables  8A  and  88)  "retrofitted"  for  Oa  exposure  and 
response;  or 

2.  Using  available  non-pollutant  summary  crop  models  and  adapting  them  to 
respond  to  O3  exposure,  through  further  research;  or 

3.  Using  available  non-pollutant  comprehensive  crop  models  with  aforementioned 
strategy  for  deriving  summary  models,  and  then  adapting  those  models  to 
Oa  exposure. 

In  each  of  these  cases,  we  also  recommend  that  the  summary  models  be  structured 
in  terms  of  the  probabilities  of  responses  derived  from  probability  density  functions  of 
the  input  parameters  and  forcing  functions.  This  would  enable  the  eventual  users  of 
these  tools  to  judge  the  risks  they  choose  to  take  in  managing  control  technologies  and 
strategies  in  their  efforts  to  protect  against  initiating  deleterious  responses  in 
vegetation . 

Option  (1)  is  recommended  for  obtaining  a  regional  photochemical  oxidant  - 
forest  response  summary  model,  by  starting  with  the  SILVA  model  (Kercher  and  Axelrod 
1981).  In  this  effort,  the  SO2  pollutant  exposure  sub-models  could  be  replaced  with  a 
three-parameter  (Nosal  1983,  1984)  exposure  model  for  oxidants.  One  or  more  summary 
models  should  then  be  derived  and  applied  on  a  regional  basis  using  "Krieging"  (with 
strict  regulations)  to  obtain  the  required  input  data. 

For  option  (2),  there  is  a  published  summary  model  called  "SUCROS"  for  small 
grains  including  wheat  (van  Keulen  and  de  Wit  1982)  which  includes  phenological  events 
(growth  stages).     This  model  could  be  used  for  Alberta.    However,  this  model  would  also 
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need  to  be  adapted  for  O3  exposures.  For  regional  scale  applications,  we  recommend, 
with  the  restrictions  described  previously,  Krieging,  three-parameter  exposure  kinetics, 
and  the  input  data  for  the  "SUCROS"  summary  model. 

For  option  (3),  from  the  evidence  provided  by  Kickert  (1984b)  for  the  sensi- 
tivity analysis  of  results  of  some  crop  models,  the  most  sensitive  factors  (Kickert, 
unpublished)  could  be  used  to  derive  simplified  summary  models,  for  example,  from  the 
CORNGRO  (corn  growth)  model.  It  must  be  realized,  however,  that:  (a)  the  kinds  of  crops 
for  which  sensitivity  data  were  described  by  Kickert  (1984b)  are  far  fewer  than  the 
number  of  different  crops  for  which  published,  comprehensive  process  models  exist  (there 
are  few  or  no  sensitivity  analysis  data  for  wheat,  barley,  oats  or  alfalfa  models, 
although  such  models  are  amenable  to  more  sensitivity  analysis  -  computer  simulation 
experiments);  and  (b)  corn  is  not  a  dominant  crop  in  Alberta,  Canada. 

If  existing  models  of  small  grains  and  alfalfa  were  to  be  used,  then  they  would 
first  have  to  be  enhanced  to  respond  to  inputs  of  ozone  exposure  data,  and  sensitivity 
analyses  would  then  have  to  be  performed.  For  analysing  potential  crop  response  to 
regional  scale  pollutant  exposure,  this  type  of  crop  plant  phenological  model  would  then 
have  to  be  reduced  to  a  phenology -based  summary  model  which  would  still  contain  a 
conceptual  foundation  based  on  processes  (rather  than  being  solely  statistical,  as  are 
many  dose-response  models). 

With  any  of  these  approaches,  a  hybridization  between  comprehensive  mechanistic 
process  models  of  crops  and  forest  trees,  and  more  simplified  summary  models  could  be 
obtained.  Such  models  would  be  capable  of  responding  to  pollutant  inputs  on  a  regional 
scale  and  for  rural  areas,  in  ways  that  would  be  more  meaningful  than  collecting  a  lot 
of  field  data  and  using  routine  statistical  regression  methods  only. 

In  any  of  these  options,  for  regional  scale  application,  we  recommend  that 
researchers  consider  using  the  Krieging  method  in  the  strictest  sense  to  derive  data  for 
three  or  four  input  parameters  needed  for  a  Mixed  Multivariate  Polynomial  -  Fourier  type 
exposure  model.  Such  information  should  be  coupled  with  other  meteorological  and  crop 
data  as  input  parameters  to  a  summary  model  of  vegetation  response  on  a  regional  basis. 

6.2  DRY  AND  WEIFALL  ACIDIC  OR  ACIDIFYING  POLLUTANTS 

Some  scientists  feel  that  dry  (acid  aerosols)  and  wetfall  acidic  pollutants 
(acidic  rain)  are  only  important  as  more  than  additive,  response  inducing  agents  with 
other  air  pollutants.  Clear  evidence  for  causation  of  direct  effects  of  exposure  to 
increased  hydrogen  ion  concentration  on  crop  and  tree  surfaces  in  actual,  natural 
terrestrial  field  environments  is  lacking.  There  are  two  mechanistic  process  models 
(Andersson  et  al.  1980;  Haines  and  Waide  1980)  and  one  statistical  model  (Medeiros  et 
al.  1983)  for  acidic  precipitation,  but  these  are  not  sufficiently  developed  to  lead  us 
to  recommend  their  use  in  other  projects.  Until,  and  unless,  evidence  for  field  effects 
on  vegetation  is  established  by  observations  using  physiological  ecology  and  pathology, 
truthful,  reliable,  and  useful  vegetation  response  models  cannot  be  assembled.  Unfor- 
tunately, standard  precipitation  meteorology  observations  do  not  include  leaf  wetness 
measurements  as  influenced  by  dew  condensation  or  fog  interception.  These  processes  of 
surface  wetting  followed  by  exposure  to  dryfall  appear  to  be  a  significant  mechanism  in 
inducing  vegetation  effects. 
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6.2.1  Local.  Point-Source  Oriented  Applications 

For  the  assessment  of  pollutant  effects  on  crop  and  tree  species  in  local  areas 
around  point  sources  of  sulphur,  oxides  of  nitrogen,  and  pollutant  mixtures,  comprehen- 
sive mechanistic  process  models  of  vegetation  growth  could  be  used.  In  these  studies, 
only  a  few  to  several  field  study  sites  would  be  needed  for  application,  in  contrast  to 
regional  scale  applications. 

6.2.2  Sulphur  Dioxide  or  Hydrogen  Sulphide 

For  gaseous  sulphur  pollutants  and  responses  of  small  grain  crops,  we  recommend 
that  Coughenour's  (1981)  model  (Table  8B)  be  adapted  to  the  needs  of  an  assessment  pro- 
gram. Such  an  approach  should  include  enhancements  to  the  pollutant  exposure  interaction, 
to  account  for  the  discontinuous,  fluctuating  air  quality  regimes  of  varying  durations 
and  peaks. 

In  forest  environments  with  SO2  exposures,  we  recommend  the  adaptation  of  the 
Kercher  and  Axelrod  (1981)  SUVA  model.  This  needs  to  be  retrofitted,  however,  with  more 
realistic  air  pollutant  exposure  kinetics  such  as  those  in  the  Nosal  (1983,  1984)  models. 

6.2.3  Oxides  of  Nitrogen  and  Pollutant  Mixtures 

Unfortunately,  we  cannot  recommend  any  biological  process  models  of  field 
effects  caused  by  oxides  of  nitrogen  or  pollutant  mixtures  for  crop  or  tree  species 
because  we  have  not  found  any.  The  statistical  models  for  nitrogen  dioxide  (Heck  and 
Tingey  1979)  and  for  pollutant  mixtures  (Medeiros  et  al.  1983;  Rowe  and  Chestnut  1985) 
are  the  only  ones  in  this  category  of  which  we  are  aware. 
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8.  APPENDICES 

The   appendices   contain   the  mathematical    equations   for  some  of  the 
models  referenced  in  Tables  7A,  B  and  OA,  B. 


Equations  for  empirical  statistical  models  shown  in  Table  7: 


Benson  et  al.  (1982) 

Heck  and  Tingey  (1979) 

Heck  et  al.  (1984b) 

Loehman  and  Wilkinson  (1983) 

Nosal  (1983) 

Oshima  et  al.  (1976) 

Stevens  and  Hazelton  (1976) 


Heagle  et  al.  (1986) 
Heck  et  al.  (1982) 
Larsen  and  Heck  (1976) 
Medeiros  et  al .  (1983) 
Nosal  (1984) 

Rowe  and  Chestnut  (1985) 


Equations  for  mechanistic  process  models  shown  in  Table  8: 

Coughenour  (1981)  Kercher  (1980) 

Kercher  and  Axelrod  (1981)  King  et  al .  (1983) 
Mortensen  (1984) 


Computer   program   equations    for   the   SUCROS   Summary  Model    of  small 
grains  (non-pollution): 

van  Keulen  et  al .  (1982) 


Those  models  for  which  we  could  not  find  published  mathematical 
equations,  or  statistical  regression  equations  for  which  coefficients  were  not 
published,  are  not  included  in  the  appendices. 


For  specific  literature  citations  included  within  any  given  model,  the 
reader  is  referred  to  the  original  publication  of  the  model  itself,  as  cited 
in  this  report. 
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8.1  EQUATIONS  FOR  EMPIRICAL  STATISTICAL  MODELS 

Benson,  F.J.,  S.V.  Krupa,  P.S.  Teng,  and    D.E.  Welsch.    1982.    Economic  assessment  of  air 

pollution  damage  to  agricultural  and  si  1 vicultural  crops  in  Minnesota.  Final 

Report  to  Minnesota  Pollution  Control  Agency,  Roseville,  Minnesota.    270  pp. 


In  general,  the  effect  of  any  stress  factor  at  one  point  in  a  crop's  growth  may 
be  expressed  in  terms  of  a  proportion  of  response  relative  to  a  check.  This  check  may 
be  the  crop's  response  under  conditions  of  zero  stress  or,  in  the  case  of  ozone  and  sul- 
phur dioxide,  conditions  with  background  ambient  concentrations.  The  general  empirical 
model  is  expressed  as: 

Y  =  f(X^  )  ...  [1] 

where 

Y  =  proportion  of  reduction  in  response  parameter,  e.g.,  dry  matter  or  yield, 

and 

X-j      stress  dosage  parameter  at  time  t-j . 

With  linear  least  squares  approximation,  an  example  of  a  functional  relationship 

may  be: 

Y  =  bo  +  biX  4-  b2X2  ...  [2] 

where 

Y  -  proportion  of  yield  loss,  and 

X    =  dose  at  a  crop  development  stage. 


To  develop  a  method  for  estimating  crop  response  over  a  growing  season,  it  is 
essential  that  the  dynamic  "episodal"  nature  of  pollutant  concentrations  be  accounted 
for.  The  approach  adopted  is  a  dynamic  matrix  model,  where  for  each  specified  time 
interval,  dt,  a  functional  relationship  was  developed,  as  in  equation  [2].  A  proportion 
of  yield  reduction.  Ay,  is  incurred  during  each  time  interval.  The  net  crop  response 
over  a  growing  season  is: 

n 

J  y  •  dt 

1 

where 

n  =  maximum  number  of  growing  days. 


The  time  step  for  integration,  dt,  was  dependent  on  the  specific  crop  and  the 
raw  data  used  in  the  model  development.  For  example,  with  alfalfa,  where  pollutants 
have  a  direct  effect  on  the  economic  biomass,  dt  ^  1  day.  With  a  cereal,  where  pollutants 
have  an  indirect  loss  effect,  i.e.,  the  plant  part  affected  by  pollution  is  different 
from  the  part  with  economic  value  (the  grain),  dt  was  a  larger  time  step. 


continued . 
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Three  summary  statistics  were  tested  for  characterizing  pollutant  dosage  on  a 
daily  basis.    These  were: 

n 

I  hi 

i=^l         the  mean  dose  per  hour  when  dose  exceeds  ambient  background, 
n 

n 

I  hi 

i=^l         the  mean  dose  per  hour  over  a  standardized  24  hour  period. 


24 
n 

z  h-j  the  total  dose  per  day. 
i-1 


Depending    on    the    size    of    the    time    interval,    these    daily    statistics  were 
combined  over  the  appropriate  number  of  days. 
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Heagle,  A.S.,  V.M.  Lesser,  J.O.  Rawlings,  W.W.  Heck,  and  R.B.  Philbeck.  1986.  Responses 
of  soybeans  to  chronic  closes  of  ozone  applied  as  constant  or  proportional 
additions  to  ambient  air.     Phytopathology  76:51-56. 


Equations  for  the  regression  of  seed  yield  of  cultivar  Davis  soybeans  on  chronic 
doses  of  O3  added  in  constant  or  proportional  amounts  to  the  O3  in  ambient  air. 

Polynomial  models:^ 

Constant  addition;!^ 

y  =  549  -  2723X 
(16)  (226) 

Proportional  addition: 

y  =  426  +  3637X  -  59,814x2 
(28)     (1174)  (10,482) 

Ignoring  method  of  addition  (reduced  model): 

y  =  357  +  9514X  -  200,060x2  ^.  945,529x3 
(64)     (3987)     (68,869)  (352,771) 

Full  (segmented)  polynomial  models:^ 

y  =  469  -  44x  (for  x  <  0.044  ppm) 

(33)  (891) 

y  =  469  -  3630(x  -  0.044)  (for  x  >  0.044  and  constant 

(33)  (362)  addition) 

y  -  469  -  4112(x  -  0.044)  (for  x  >  0.044  and 

(33)    (467)  proportional  addition) 

Weibull  models:*^ 

Constant  addition: 


y  _-  490  e-<x^0-126)'-' 
(20,  0.007,  0.44) 

Proportional  addition: 

y  .  479  e  -(x/0-103)3. 
(12,  0.003,  0.67) 

Ignoring  method  of  addition: 

y=492e-(^/0-^2^)'- 
(21 ,  0.006,  0.43) 


continued. 
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y  =  yield  (g)  of  seeds  per  metre  of  row;  x  -  O3  dose  characterized  as  the  seasonal 
7  h/day  mean  concentration  (ppm).  Standard  errors  for  each  parameter  estimate  are 
shown  in  parentheses. 

The  quadratic  equation  (with  standard  errors  in  parentheses)  for  constant  additions 
was  : 

y  =  500  (31)  -  694x  (1063)  -  15.901x2  (8145). 

The  segmented  polynomial  model  provides  for  a  linear  response  (b-|  =  -43. 8x)  to  Oa 
levels  <  0.044  ppm  and  a  different  linear  response  above  0.044  ppm  for  constant 
addition,  -    -3650,    and    for   proportional    addition,    bp   =^    -4112,    and    x   =^  Oa 

concentration. 


^     For  the  Weibull    model  (y  =^  a  e'-  ) ;  y  =  estimated    yield  of    seeds  (g)  per  metre 

of  row;  X  =  Oa  dose  as  seasonal  7  h/day  mean  concentration  (ppm);  a  =  maximum  seed 
yield  at  0  ppm  Oa;  <?  =  Oa  concentration  at  which  a  is  reduced  by  63%;  and  c  is  a 
dimensionless  shape  parameter.  Standard  errors  in  parentheses  are  for  a,  o,  and 
c,  respectively. 


Heck,  W.W.  and  D.T.  Tingey.  1979.  Nitrogen  dioxide:  time-concentration  model 
to  predict  acute  foliar  injury.  U.S.  Environmental  Protection  Agency, 
EPA-600/3-79-057.     Corvallis,  Oregon.  16  pp. 


Time-concentration  response  equations  for  a  selected  group  of  plants  exposed 
to  nitrogen  dioxide.^ 


Concentrations  (ppm)  to  produce  the 
I  (in  %)  in  T  (h) 


Plants  by  Name 
(Common,  Cultivars) 


Equation 
(C  -Aq  ^-A-|  I  fA2/T)b 


1-5,  1=5, 
T=l  1=8 


1  =  50, 
T=l 


1  =  50, 
T=8 


Susceptible 

Oats 

C=l  , 

.45+0. 

,13 

1+2. 

,39/T 

0. 

,76 

4. 

,5 

2, 

,3 

10, 

,3 

8. 

,3 

(Clintland  64) 

Radish 

C=2, 

.43f0. 

.14 

I+l  . 

,02/T 

0. 

,83 

4, 

,1 

3, 

,2 

10. 

,5 

9. 

.6 

(Cherry  Bel le) 

Oats 

C=l  , 

.75+0. 

.15 

1+3. 

,24/T 

0. 

.56 

5. 

,7 

2. 

,9 

12, 

,5 

9. 

,7 

(329-80) 

Bromegrass 

C=2 

.49+0. 

,16 

I+l  , 

.90/f 

0. 

.71 

5, 

,2 

3, 

,5 

12, 

,4 

10. 

,7 

(Sac  Smooth) 

Begonia 

C=2, 

.45+0. 

,15 

1+2, 

,99/T 

0, 

.63 

6, 

,2 

3. 

.5 

12, 

.9 

10. 

,3 

(Thousand  Wonders 

White) 

Chrysanthemum 

C=3, 

.16+0. 

,16 

1+2, 

,14/T 

0. 

,72 

6, 

.1 

3. 

,7 

13. 

.3 

11  , 

.4 

(Oregon) 

Oats 

C  =  2 

.79+0. 

,14 

1+2. 

,88/T 

0. 

,50 

6, 

,4 

3 

.8 

12. 

.7 

10, 

.2 

( Pendek) 

Wheat 

C  =  2 

.80+0, 

,13 

1+2, 

.94/1 

0. 

,52 

6. 

,4 

3, 

.8 

10. 

.2 

7, 

.7 

(Wei  Is) 

Sultana 

C  =  3 

.93+0, 

,13 

I+l  . 

.73/T 

0, 

,67 

6, 

.3 

4. 

.8 

12. 

.2 

10, 

.7 

(White  Imp) 

Broccol i 

C  =  3 

.07+0. 

,20 

1+2. 

.94/1 

0, 

,53 

7. 

.0 

4. 

.5 

16. 

.0 

13, 

.4 

(Calabreese) 

Periwinkle 

C=2 

.92+0, 

,23 

1+3. 

.02/T 

0, 

,55 

7, 

.1 

4. 

.5 

17. 

.4 

14, 

.8 

(Bright  Eyes) 

Intermediate 
Cotton 

( Paymaster) 
Cotton 

(Acala  4-42) 
lobacco 

(Bel  B) 
lobacco 

(Bel  W3) 
lobacco 

(White  Gold) 


C=2, 

.97+0. 

,23 

I+l  . 

.94/1 

0. 

.58 

7. 

.1 

5. 

,3 

17. 

,4 

15, 

.7 

C=3. 

.68+0. 

,22 

1+3. 

.15/1 

0. 

.50 

7. 

.9 

5, 

,2 

17. 

,8 

15, 

.1 

C=3. 

.62+0. 

.21 

1+3. 

.98/1 

0. 

38 

8. 

.7 

5. 

,2 

18. 

,1 

14, 

.6 

C=3 

.65+0. 

,18 

1+4. 

.40/1 

0. 

.31 

9 

.0 

5. 

,2 

17. 

,1 

13. 

.2 

C=4 

.03+0. 

,30 

1+3 

.56/1 

0. 

.40 

9 

.1 

6. 

,0 

22. 

,6 

19. 

.5 

Tolerant 

lobacco  None 

(Burley  21) 
Corn  None 

(Pioneer  509-W) 
Corn  None 

(Golden  Cross) 
Azalea  None 

(Alaska) 

Sorghum  (Martin)  None 
Cucumber  (Long  Marketer)  None 


1  of  36  plants  injured; 

0.5  h,  26  ppm,  33  percent. 
4  of  36  plants  injured; 

all  injuries  were  1.6  percent, 
1  of  36  plants  injured; 

7.0  h,  6  ppm,  7  percent. 
1  of  36  plants  injured; 

1.0  h,  17  ppm,  5  percent. 
0  of  36  plants  injured. 
0  of  36  plants  injured. 


conti  nued . 


Equations  were  developed  from  exposures  limited  in  time  (0.6  -  7.0  h)  and 
denote  acute  injury  symptoms  to  the  plants.  Concentrations  used  ranged 
from  1  to  20  ppm  of  nitrogen  dioxide.  Plants  are  grouped  in  three 
susceptibility  categories. 

C  is  nitrogen  dioxide  concentration  in  ppm;  I  is  percent  injury;  T  is  time 
in  h;  and  Aq,  A-] ,  and  A2  are  constants  (partial  regression  coef- 
ficients) specific  for  pollutant,  plant  species,  and  environmental 
conditions  used. 

r2,  multiple  correlation  coefficient  squared  which  represents  the  percent 
variation  explained  by  the  model. 


Heck,  W.W.,    O.C.  Taylor,    R.  Adams,    G.  Bingham,  J.  Miller,  E.  Preston,  and  L. 

Weinstein.  1982.  Assessment  of  crop  loss  from  ozone.  Journal  of  the 
Air  Pollution  Control  Association  32(4):  353-361. 


Linear  regression  models  of  crop  yield  as  a  function  of  Oa  concentration. 


Predicted  Yield^ 
in  g/plant  (+SE) 


Linear  Functions^  0.025       0.06  0.10 

Crop  (±SE)'^  ppm  ppm  ppm 


Corn 

y 

= 

247 

8 

-    260  X 

0 

65 

241  .3 

232.2 

221 

8 

(Coker  16) 

(±6 

3) 

(±68) 

(+5.2) 

(+3.2) 

(±3 

4) 

Soybean 

y 

= 

23 

5 

-    126  X 

0 

96 

20.3 

15.9 

10 

9 

(Corsoy) 

(±0 

7) 

(±10) 

(+0.5) 

(+0.3) 

(±0 

5) 

Soybean  -  1977 

y 

- 

96 

6 

-    385  X 

0 

90 

87.0 

73.5 

58 

1 

(Davis) 

(±7 

5) 

(±78) 

(+5.9) 

(+4.0) 

(+3 

6) 

Soybean  -  1978 

y 

95 

3 

-    309  X 

0 

99 

87.6 

73.8 

64 

4 

(Davis) 

(±9 

3) 

(+109) 

(+6.9) 

(+4.4) 

(±4 

7) 

Kidney  Bean 

y 

= 

16 

8 

31  X 

0 

84 

16.1 

14.9 

13 

7 

(Calif.  Light  Red) 

(±0 

6) 

(±7) 

(+0.5) 

(+0.4) 

(+0 

3) 

Lettuce,  Head 

y 

1065 

7 

-  5978  X 

0 

94 

916.2 

707.0 

467 

9 

( Empi  re) 

(+66 

4) 

(±701) 

(+50.3) 

(+30.5) 

(+24 

2) 

Peanut^ 

y 

= 

173 

2 

-  1046  X 

0 

99 

147.1 

110.4 

68 

6 

(NC-6) 

(±3 

6) 

(±42) 

(+2.6) 

(+1.6) 

(±1 

8) 

Spinach 

y 

22 

7 

-    106  X 

0 

98 

20.1 

16.3 

12 

1 

(America) 

(±2 

1 ) 

(±23) 

(+1.6) 

(±1.1) 

(±1 

1) 

Spinach 

y 

42 

1 

-    193  X 

0 

93 

37.3 

30.5 

22 

8 

(Hybrid  7) 

(±4 

9) 

(+53) 

(+3.8) 

(+2.6) 

(±2 

6) 

Spinach 

y 

46 

1 

-    238  X 

0 

94 

40.1 

31  .8 

22 

3 

(Virof lay) 

(±4 

5) 

(±27) 

(+3.5) 

(+2.4) 

(+2 

4) 

Spinach 

y 

23 

3 

-    121  X 

0 

996 

20.2 

16.0 

11 

2 

(Winter  Bloomsdal 

e) 

(±8 

5) 

(+49) 

(+1 .9) 

(+1.3) 

(±1 

3) 

Turnip^ 

y 

12 

9 

94  X 

0 

86 

10.6 

7.3 

3 

5 

(Just  Right) 

(±0 

5) 

(+8) 

(+0.4) 

(+0.2) 

(±0 

4) 

Turnip^  (Purple 

y 

7 

2 

49  X 

0 

94 

6.0 

4.3 

2 

3 

Top  White  Globe) 

(±0 

3) 

(±5) 

(+0.2) 

(+0.2) 

(±0 

3) 

Turni  p 

y 

5 

3 

-      36  X 

0 

89 

4.4 

3.1 

1 

7 

(Shogoin) 

(±0 

3) 

(±4) 

(+0.2) 

(+0.1) 

(±0 

2) 

Turnipd 

y 

18 

1 

-    116  X 

0 

75 

15.2 

11.1 

6 

5 

(Tokyo  Cross) 

(±1 

5) 

(±24) 

(±1-1) 

(+0.8) 

(±1 

3) 

Wheat 

y 

6 

6 

18  X 

0 

92 

6.2 

5.5 

4 

8 

(Blueboy  II) 

(±0 

3) 

(±3) 

(+0.2) 

(+0.1) 

(±0 

1) 

Wheat 

y 

5 

.8 

-      21  X 

0 

98 

5.3 

4.5 

3 

7 

(Coker  47-27) 

(±0 

•  l ) 

(±2) 

(+0.2) 

(+0.1) 

(+0 

1) 

Wheat^ 

y 

5 

.7 

16  X 

0 

82 

5.3 

4.7 

4 

1 

(Holly) 

(±0 

•3) 

(±3) 

(+0.2) 

(+0.1) 

(±0 

1 ) 

Wheat 

y 

4 

.9 

12  X 

0 

88 

4.6 

4.2 

3 

7 

(Oasis) 

(±0 

.2] 

(±2) 

(+0.2) 

(+0.1) 

(±0 

1 ) 

continued . 


Linear  regression  model  of  the  form  y  =  bo+b-|X,  where  y  =  yield  in  g/plant 
and  X  =  seasonal  7-h/day  mean  Oa  concentration. 

The  standard  error  term  was  calculated  using  the  mean  square  error  from 
the  analysis  of  variance. 

Values  may  differ  slightly  from  those  predicted  by  the  equations  due  to  a 
rounding  of  numbers  in  the  functions  shown.  The  0.06  ppm  seasonal  7-h/day 
mean  Oa  concentration  is  the  maximum  mean  concentration  expected  in  many 
parts  of  the  US  when  the  one-hour  standard  of  0.12  ppm  is  just  met.  The 
0.10  ppm  seasonal  value  was  arbitrarily  used  as  a  maximum  concentration  to 
consider  and  approximates  the  NF-3  treatment. 

More  complex  models  account  for  significantly  more  of  the  variation  than 
these  simple  linear  models. 
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Plateau-linear  models  of  crop  yield  as  a  function  of  O3  concentration.^ 


Predicted  Yield*^ 

(g/plant)  Predicted 

  %  Yield 

Reduction  at 

0.025  ppm    0.06  ppm    0.10  ppm  x  =  0.06  ppm 
Plateau-Linear  Functions^    r2         (+$£)         (±SE)         (±SE)  (±SE) 


Peanut  (NC-6) 

y=142.3  if  x<e=0.037 


184.6-1160X  if  x>0.037 


Turnip    (Just  Right) 

y=10.7  if  x<e  =0.038 
y=15.5-127x  if  x>0.038 


0.99 


0.96 


142.3 
(±3.3) 


10.7 
(±0.4) 


115.0 
(±3.6) 


7.8 
(±0.6) 


68.6 
(±5.2) 


2 

(±1 


19.2 
(±1.1) 


26.6 
(±4.1) 


Turnip    (Purple  Top  White  Globe) 

y=6.0  if  x<e  -0.034  0.99  6.0  4.4  2.0  26.0 

y=8.1-60x  if  x>0.034  (±0.3)        (±0.4)        (±0.6)  (±3.4) 


Turnip    (Tokyo  Cross) 

y=l 4.8  if  x<e  -0.054  0.94         14.8  12.7  4.4 

y=27.0-226x  if  x>0.054                    (±1-2)  (±1.2)  (±2.9) 

Wheat  (Holly) 

y=4.9  if  x<e  =0.087  0.99           4.9  4.9  4.4 

y=8.2-38x  if  x>0.087                       (±0.1)  (±0.1)  (±0.2) 


13.9 
(±3.8) 


0 

(-) 


3  This  includes  only  those  data  sets  for  which  the  fit  was  significantly 
improved  over  the  linear  model.  The  O3  is  expressed  as  the  seasonal 
7-h/day  mean  concentration. 

b  y=bo  if  x<e;  y=( bg-b-i e) x  if  x>e.  The  e  is  the  threshold  concentration 
(ppm)  below  which  yield  (y)  is  unchanged  at  a  value  bg;  b-)  is  the  slope  of 
the  linear  portion  if  the  O3  concentration  (x)  is  greater  than  e. 

c  Values  may  differ  slightly  from  those  predicted  by  the  equations  due  to  a 
rounding  of  numbers  in  the  functions  shown.  The  0.06  ppm  seasonal  7-h/day 
mean  Oa  concentration  is  the  maximum  mean  concentration  expected  in  many 
parts  of  the  US  when  the  one-hour  standard  of  0.12  ppm  is  just  met,  while 
0.10  ppm  was  arbitrarily  used  as  a  maximum  7-h/day  concentration  for 
comparison. 
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Percent  yield  reduction  as  a  function  of  Oa  concentration. 


Predicted  %  Yield 
Reduction  at  Two  O3 
Concentrations 
(±SE) 


Crop  Linear  Function^*  0.06  ppm        0.10  ppm 


Corn    (Coker  16) 

Y 

-2.7  +  108  X 

3, 

.83  H 

h  1  .0 

8 

.1 

2. 

.1 

Soybean  (Corsoy) 

Y 

-15.5  +  621  X 

21  . 

.7 

+ 

1 

.4 

46 

.5 

+ 

2. 

.9 

Soybean  -  1977  (Davis) 

Y 

-11.1  +  443  X 

15, 

.5 

+ 

2 

.4 

33 

.2 

+ 

5, 

.1 

Soybean  -  1978  (Davis) 

Y 

-8.8  +  353  X 

12. 

.3 

3 

.6 

26 

.4 

+ 

7. 

.7 

Kidney  Bean    (Calif.  Light  Red) 

Y 

-4.8  +  193  X 

6. 

.8 

1 

.3 

14 

.5 

+ 

2. 

.8 

Lettuce,  Head  (Empire) 

Y 

-16.3  +  652  X 

22. 

.8 

+ 

1  , 

.7 

48, 

.9 

+ 

3. 

.6 

Peanut  (NC-6) 

Y 

-17.8  +  711  xC 

24. 

.9 

+ 

0 

.7 

53 

.3 

1  . 

.5 

Spinach  (America) 

Y 

-13.2  +  527  X 

18, 

.5 

+ 

3 

.0 

39 

.6 

6. 

.4 

Spinach  (Hybrid  7) 

Y 

-13.0  +  517  X 

18. 

,1 

3, 

.7 

38 

.9 

8. 

,0 

Spinach  (Viroflay) 

Y 

-14.8  +  594  X 

20, 

,7 

3. 

.1 

44 

.4 

+ 

6, 

.6 

Spinach  (Winter  Bloomsdale) 

Y 

-14.9  +  599  X 

20, 

,9 

+ 

3, 

.4 

44. 

.8 

+ 

7, 

.3 

Turnip  (Just  Right) 

Y 

-22.1  +  887  xC 

31  , 

,0 

+ 

2, 

.1 

66. 

.4 

4, 

.5 

Turnip  (Purple  Top  White  Globe) 

Y 

-20.5  +  817  xC 

28, 

,6 

+ 

2, 

.4 

61  . 

.4 

5, 

.0 

Turnip  (Shogoin) 

Y 

-20.5  +  818  X 

28, 

,6 

+ 

2, 

.6 

61  . 

.4 

+ 

5, 

,5 

Turnip  (Tokyo  Cross) 

Y 

-19.1  +  763  xC 

26. 

,8 

4, 

.6 

57. 

.4 

+ 

9. 

.8 

Wheat  (Blueboy  II) 

Y 

-7.4  +  290  X 

10. 

.4 

1  , 

,4 

22. 

,3 

-«- 

3. 

,0 

Wheat  (Coker  47-27) 

Y 

-10.1  +  405  X 

14. 

.2 

+ 

1  . 

,2 

30. 

,4 

2. 

.6 

Wheat  (Holly) 

Y 

-7.5  +  304  xC 

10. 

.6 

1  . 

.5 

22. 

,6 

3. 

,3 

Wheat  (Oasis) 

Y 

-6.2  +  250  X 

8. 

7 

+ 

1  . 

.7 

18. 

.6 

+ 

3. 

.7 

This  includes  all  data  sets.  The  O3  is  expressed  as  the  seasonal  7-h/day 
mean  concentration. 

100  b] 

Y  -    [0.025-x] 

a 

where  Y  =  predicted  percent  yield  reduction,  b]  =^  the  regression  coef- 
ficient from  the  yield  model,  a  =  the  predicted  yield  at  a  seasonal  7-h/day 
mean  O3  concentration  of  0.025  ppm  and  x  ^  seasonal  7-h/day  mean  Oa 
concentration. 

These  models  must  be  considered  with  caution.  More  complex  models  account 
for  significantly  more  of  the  variation  than  these  simple  linear  models. 
Compare  the  predicted  percent  yield  reductions  for  the  plateau-linear  model 
at  0.06  ppm  to  get  an  idea  of  the  differences. 

Values  may  differ  slightly  from  those  predicted  by  the  equations  due  to  a 
rounding  of  numbers  in  the  functions  shown.  The  0.06  ppm  seasonal  7-h/day 
mean  Oa  concentration  is  the  maximum  mean  concentration  expected  in  many 
parts  of  the  US  when  the  one-hour  standard  of  0.12  ppm  is  just  met;  the 
0.10  ppm  was  arbitrarily  used  as  a  maximum  7-h/day  mean  Oa  concentration 
for  comparison. 
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Heck,  W.W.,      W.W.  Cure,      J.O.  Rawlings,      L.J.  Zaragoza,      A.S.  Heagle,  H.E. 

Heggestad,  R.J.  Kohut,  L.W.  Kress,  and  P.J.  Temple.  1984.  Assessing 
impacts  of  ozone  on  agricultural  crops:  II.  Crop  yield  functions 
and  alternative  exposure  statistics.  Journal  of  the  Air  Pollution 
Control  Association  34:  810-817. 


Estimates  of  the  parameters  and  the  residual  mean  squares  (RMS)  from 
fitting  the  Weibull  model  using  four  different  exposure  statistics  to  describe 
Os  concentration. 


Parameter  Estimates  for  Ozone  Exposure  Statistics^ 

Estimated   

Parameters 

and  RMS^  M7  P7  Ml  PI 


Barley 

'Poco'C,  1982  -  Shafter,  CA 

a  (se)           1.988  (0.051)  1.991  (0.049)  1.990  (0.050)  1.989  (0.052) 

a  (se)           0.205  (0.669)  0.322  (0.744)  0.246  (0.651)  0.466  (1.562) 

c  (se)           4.278  (17.15)  4.231  (12.66)  4.385  (15.26)  3.724  (13.51) 

RMS  (21)        0.0248  0.0245  0.0246  0.0248 


Bean,  kidney 

'Calif.  Light  Red'd, 


(Full  Plots  -  FP),  1982  -  Ithaca,  NY 


a  (se) 
o  (se) 
c  (se) 
RMS  (7) 


2878  (416) 
0.120  (0.013) 
1.171  (0.489) 
30354 


3586  (1444) 
0.171  (0.071) 
0.960  (0.641) 
30779 


3394  (1134) 
0.130  (0.048) 
0.865  (0.533) 
31489 


3776  (2025) 
0.183  (0.100) 
1 .000  (0.779) 
28870 


'Calif.  Light  Red'^,  (Partial  Plots  -  PP) 

i  (se)  2447  (183)           2437  (122) 

a  (se)  0.108  (0.007)      0.187  (0.007) 

c  (se)  3.942  (1 .733)      4.727  (1 .456) 

RMS  (7)  100130  116257 


1982  -  Ithaca,  NY 

2444  (122)  2432  (127) 

0.140  (0.006)  0.242  (0.011) 

3.88  (1 .164)  4.620  (1 .839) 

103095  215898 


Corn 


'PAG  397 
a  (se) 
a  (se) 
a  (se) 
RMS  (17) 


1981  -  Argonne, 
13968  (322) 
0.160  (0.004) 
4.280  (0.723) 
2644143 


IL 

13928  (311) 
0.244  (0.004) 
6.912  (1  .161) 
2391941 


13971  (324) 
0.182  (0.005) 
4.349  (0.730) 
2463065 


13899  (31  5) 
0.248  (0.004) 
7.401  (1.293) 
2505320 


'Pioneer  3780',  1981  -  Argonne,  IL 

a  (se)  12533  (323)  12524  (325)  12552  (329)  12539  (330) 

a  (se)  0.155  (0.004)      0.220  (0.004)  0.176  (0.005)  0.244  (0.004) 

c  (se)  3.091  (0.461)      4.811  (0.734)  3.142  (0.467)  4.987  (0.771) 

RMS  (17)        774833  697778  773206  740201 

Cotton 

'Acala  SJ-2',  1981  -  Shafter,  CA,  ( Irrigated-I) 

a  (se)  5546  (356)  5699  (489)  5526  (379)  5746  (506) 

a  (se)  0.199  (0.020)      0.253  (0.026)  0.211  (0.020)  0.268  (0.027) 

c  (se)  1.228  (0.402)      1.182  (0.415)  1.315  (0.419)  1.214  (0.423) 

RMS  (9)         80907  117382  85315  113000 


'Acala  SJ-2',  1981  -  Shafter,  CA,  (Droughted  -  D) 

a  (se)  3323  (278)  3323  (292)  3317  (272)  3323  (387) 

a  (se)  0.269  (0.076)      0.351  (0.101)      0.280  (0.073)  0.366  (0.128) 

c  (se)  1.969  (1.308)      1.980  (1.286)      2.131  (1.389)  2.027  (1.706) 

RMS  (9)         232623  267524  237305  261563 
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Estimated 
Parameters 
and  RMSa 


Parameter  Estimates  for  Ozone  Exposure  Statistics^ 


M7 


P7 


Ml 


PI 


'Acala  SJ-2'e,  1982  -  Shafter,  CA,  (Irrigated  -  I) 
a  (se)  5872  (555)  6192  (990)  5901  (580) 

a  (se)  0.088  (0.007)      0.174  (0.023)      0.108  (0.009) 

c  (se)  2.100  (0.796)      1.374  (0.700)      2.084  (0.798) 

RMS  (9)         417881  436672  427681 


6321  (1038) 
0.211  (0.032) 
1.197  (0.587) 
524417 


'Acala  SJ-2'6,  1982  -  Shafter,  CA,  (Droughted  -  D) 

a  (se)  4938  (412)  4934  (403)  4947  (417) 

a  (se)  0.094  (0.008)      0.187  (0.018)  0.115  (.009) 

c  (se)  3.333  (1.537)      2.860  (1.504)  3.355  (1.551) 

RMS  (9)         443277  395639  441216 

'Stoneville  213'^',  1982  -  Raleigh,  NC 

a  (se)  3686  (140)  3982  (252)  3732  (154) 

a  (se)  0.112  (0.004)      0.161  (0.006)  0.127  (0.004) 

c  (se)  2.577  (0.416)      2.005  (0.382)  2.638  (0.440) 

RMS  (21)        52171  63823  50211 


4926  (404) 
0.213  (0.020) 
2.997  (1  .572) 
404471 


3828  (220) 
0.175  (0.005) 
3.061  (0.580) 
93827 


Peanut 

'NC-6' 
a  (se) 
o  (se) 
c  (se) 
RMS  (17) 


1980  -  Raleigh, 
7485  (231) 
0.111  (0.002) 
2.249  (0.225) 
109603 


NC 


7705  (299) 
0.167  (0.004) 
2.272  (0.254) 
168540 


7586  (258) 
0.128  (0.003) 
2.180  (0.228) 
122330 


8381  (501) 
0.186  (0.006) 
2.025  (0.264) 
251656 


Sorghum 

'DeKalb-28' 
a  (se) 
o  (se) 
c  (se) 
RMS  (17) 


1982  -  Argonne,  XL 


8137  (218) 
0.296  (0.135) 
2.21  7  (1  .229) 
195875 


8122  (201  ) 
0.305  (0.075) 
3.984  (2.098) 
197090 


8156  (235) 
0.361  (0.171) 
2.161  (1.228) 
197626 


8138  (211) 
0.349  (0.093) 
3.924  (2.174) 
205486 


Soybean 

' Corsoy ' 
a  (se) 
a  (se) 
c  (se) 
RMS  (7) 

'  Davis ' ^ 
a  (se) 
a  (se) 
c  (se) 
RMS  (43) 


1980  -  Argonne, 
2785  (179.3) 
0.133  (0.010) 
1  .952  (0.587) 
32514 

1981  -  Raleigh, 
5593  (863) 
0.128  (0.019) 
0.872  (0.284) 
91452 


IL 


2817  (201) 
0.197  (0.013) 
2.122  (0.663) 
37589 


NC 


5729  (1124) 
0.205  (0.035) 
1 .000  (0.378) 
141738 


2808  (198) 
0.167  (0.012) 
1  .917  (0.595) 
33667 


6094  (1366) 
0.139  (0.038) 
0.778  (0.305) 
97635 


'Davis'g.h^  1982  -  Raleigh,  NC,  (Constant  add.  -  CA) 

a  (se)  4931  (227)  5308  (544)  4968  (250) 

a  (se)  0.128  (0.008)  0.248  (0.023)      0.151  (0.009) 

c  (se)  2.144  (0.500)  1.369  (0.501)      2.129  (0.578) 

RMS  (7)  141799  415352  150747 


2799  (196) 
0.242  (0.013) 
2.572  (0.798) 
28457 


6620  (3626) 
0.201  (0.111) 
1 .000  (0.736) 
236429 


5473  (737) 
0.231  (0.015) 
1 .628  (0.711) 
247265 


'Davis'g,  1982  -  Raleigh,  NC,  (Proport.  add.  -  PA) 

a  (se)  4805  (106)  4917  (146)           4810  (108)  5037  (232) 

a  (se)  0.103  (0.003)  0.233  (0.009)      0.137  (0.004)  0.267  (0.011) 

c  (se)  4.077  (0.648)  2.554  (0.412)      3.763  (0.619)  2.540  (0.601) 

RMS  (7)  124407  113383                  128122  144331 
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Parameter  Estimates  for  Ozone  Exposure  Statistics^ 

Estimated   

Parameters 

and  RMS^  M7  P7  Ml  PI 


'Essex'^  1981  -  Beltsville,  MD 

a  (se)  4562  (608)           4487  (544) 

a  (se)  0.187  (0.124)      0.194  (0.063) 

c  (se)  1 .543  (1 .505)      2.676  (2.128) 

RMS  (25)  232455  232643 

'Williams'^  1981  -  Beltsville,  MD 

a  (se)  4992  (659)           4954  (592) 

h  (se)  0.211  (0.120)      0.214  (0.067) 

c  (se)  1 .100  (0.901 )      1 .740  (1 .206) 

RMS  (25)  127464  124577 


4622  (600) 
0.168  (0.077) 
2.025  (1.69) 
226924 


4969  (569) 
0.174  (0.067) 
1 .729  (1  .1  71  ) 
138543 


4538  (565) 
0.207  (0.057) 
2.997  (2.184) 
227660 


5037  (619) 
0.228  (0.065) 
1  .936  (1 .290) 
114483 


'  Forrest 'i^'"' -J,  1982  -  Beltsville,  MD,  (Irrigated  -  I) 
a  (se)  4333  (366)  4391  (371)  4369  (388) 

a  (se)  0.171  (0.041)      0.224  (0.041)      0.190  (0.044) 

c  (se)  2.752  (1.872)      3.455  (2.244)      2.967  (2.067) 

RMS  (11)        374947  363968  401359 

'Forrest'^. i.j,  1982  -  Beltsville,  MD,  (Droughted  -  D) 


a  (se) 
a  (se) 
c  (se) 
RMS  (11) 


5384  (706) 
0.205  (0.041) 
1.000  (0.538) 
204273 


5669  (905) 
0.271  (0.045) 
1 .000  (0.577) 
177953 


5661  (944) 
0.216  (0.035) 
1 .000  (0.586) 
177994 


'Williams'^'.i'j,  1982  -  Beltsville,  MD,  (Irrigated  -  I) 

a  (se)  5884  (389)  5826  (326)  5913  (383) 

a  (se)  0.162  (0.017)  0.213  (0.017)  0.176  (0.016) 

c  (se)  1.577  (0.504)  2.226  (0.649)  1.865  (0.565) 

RMS  (11)  124318  132923  145116 

'Williams'^. i 'j ,  1982  -  Beltsville,  MD,  (Droughted  -  D) 

a  (se)  4480  (739)  5126  (942)  5114  (985) 

o  (se)  0.217  (0.059)  0.287  (0.064)  0.229  (0.050) 

c  (se)  1.000  (0.654)  1.000  (0.702)  1.000  (0.715) 

RMS  (11)  189004  169920  171588 

'Hodgson'^,  1981  -  Ithaca,  NY,  (Full  plots  -  FP) 

a  (se)  2590  (439)  2851  (1025)  2685  (875) 

a  (se)  0.138  (0.026)  0.188  (0.042)  0.174  (0.030) 

c  (se)  1.000  (0.609)  1.000  (0.865)  1.000  (0.903) 

RMS  (7)  41226  49037  49268 


4316  (356) 
0.244  (0.041) 
4.709  (3.518) 
422091 


5893  (1205) 
0.312  (0.046) 
1 .000  (0.648) 
181555 


5769  (302) 
0.233  (0.015) 
2.893  (0.803) 
131361 


5339  (1259) 
0.328  (0.061) 
1  .000  (0.781  ) 
159911 


3059  (1817) 
0.225  (0.100) 
1 .000  (1 .180) 

74253 


Tomato 


'Murrieta 
a  (se) 
o  (se) 
c  (se) 
RMS  (28) 


^j.k^  1981  -  Tracy,  CA 

32.9  (1  .10)  33.2  (1  .58) 

0.142  (0.026)  0.209  (0.048) 
3.807  (1 .938)  2.832  (1  .500) 
11.95  13.16 


'Murrieta'^.f ,j ,  1982  -  Tracy,  CA 

a  (se)  32.3  (1.385)        35.5  (11.44)  32.7  (1.656) 

a  (se)  0.082  (0.025)      0.589  (0.626)  0.145  (0.052) 

c  (se)  3.050  (1.871)      1.000  (1.733)  2.668  (1.825) 

RMS  (28)        10.28  10.36  10.36 


33.1  (1.36) 
0.241  (0.044) 
3.447  (1 .594) 
13.75 


36.8  (16.20) 
0.690  (0.787) 
0.977  (2.063) 
11  .01 


continued. 


Parameter  Estimates  for  Ozone  Exposure  Statisticst> 

Estimated   

Parameters 

and  RMSa  M7  P7  Ml  PI 


Wheat,  winter 

'Abe' 1 ,  1982  -  Argonne,  IL 

a  (se)           5363  (246)  5381  (242)  5392  (267)  5508  (295) 

a  (se)           0.143  (0.010)  0.202  (0.010)  0.174  (0.012)  0.228  (0.010) 

c  (se)           2.423  (0.666)  3.148  (0.824)  2.351  (0.671)  3.285  (0.834) 

RMS  (17)        161418  165453  165621  190052 

'Arthur  71',  1982  -  Argonne,  IL 

a  (se)           4684  (156)  4715  (154)  4709  (169)  4837  (190) 

a  (se)           0.148  (0.007)  0.209  (0.008)  0.180  (0.009)  0.234  (0.008) 

c  (se)           2.154  (0.411)  2.711  (0.491)  2.119  (0.418)  2.886  (0.510) 

RMS  (17)        64939  76907  63462  94348 


' Roland ' 
a  (se) 
a  (se) 
e  (se) 
RMS  (17) 


1982  -  Argonne, 
5479  (312) 
0.113  (0.005) 
1  .633  (0.288) 
91726 


IL 

5331  (229) 
0.171  (0.005) 
2.368  (0.349) 
97698 


5598  (373) 
0.136  (0.006) 
1  .572  (0.296) 
91113 


5517  (279) 
0.194  (0.005) 
2.669  (0.397) 
160598 


'Vona' ,  1982  -  Ithaca,  NY 
a  (se)  7857  (2111) 

a  (se)  0.053  (0.016) 

c  (se)  1 .000  (0.341 ) 

RMS  (17)  294084 


8750  (3286) 
0.082  (0.034) 
1 .000  (0.409) 
453900 


11133  (6890) 
0.052  (0.033) 
1 .000  (0.487) 
351597 


13998  (12737) 
0.072  (0.064) 
1.000  (0.619) 
430977 


The  estimated  parameters  (a,  a,  c)  are  given  with  standard  errors  in 
parentheses;  the  residual  mean  square  (RMS)  has  the  degrees  of  freedom  in 
parentheses . 

All  c  estimates  are  in  ppm.  With  the  following  exceptions,  all  estimates 
of  a  are  in  units  of  seed  yield  (kg/ha):  barley--seed  weight  (g  per  head); 
tomato  (both  years) --fresh  weight  (kg  per  plot);  cotton — lint  -«-  seed  weight 
(kg/ha);  peanut — pod  weight  (kg/ha).  In  cases  where  the  estimated  c  param- 
eter is  exactly  1.0,  it  has  been  bounded  from  below  to  obtain  convergence 
in  the  nonlinear  model  fitting  routine.  The  parameters  were  estimated 
from  data  not  showing  the  expected  Weibull  form.  Caution  should  be  used 
in  interpreting  these  Weibull  models.  Other  models  might  better  describe 
the  behaviour  observed  in  these  experiments.  The  a  and  c  estimates  may 
be  slightly  different  than  values  in  Heck  et  al.  because  individual  plot 
data  were  used  to  calculate  the  Oa  statistics. 

The  exposures  started  late  in  this  experiment  so  the  estimate  of  loss 
related  to  Oa  concentration  may  be  conservative. 

"Full"  plots  are  chambers  in  which  the  plants  were  left  undisturbed  for 
the  entire  growing  season  and  "partial"  plots  are  those  from  which  some 
plants  were  periodically  removed  to  obtain  physiological  data. 

Weather  during  the  1982  growing  season  in  California  was  radically  differ- 
ent from  that  in  1981;  rainfall  was  higher  and  temperatures  were  lower  in 
1982. 
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'  Sulphur  dioxide  was  part  of  these  experimental  designs.  A  bivariate 
Weibull  model  was  fit  to  the  1981  Davis,  the  1981  Essex,  and  the  Williams 
soybean  data  sets  while  for  the  other  data  sets  the  SO2  effect  was 
essentially  linear  and  was  accounted  for  by  estimation  of  a  separate  a 
parameter  to  accommodate  the  linear  effect. 

9  These  two  data  sets  are  from  the  same  experiment  and  utilized  the  same  CF 
and  NF-1  treatment  plot  yields. 

h  The  PI  exposure  statistic  for  the  data  set  would  not  permit  fitting  of  the 
Weibull  model;  therefore,  the  second  highest  peak  was  used. 

^  These  four  data  sets  are  from  a  split-plot  design  with  the  two  cultivars 
grown  in  the  same  treatment  plots.  The  homogeneity  test  of  the  proportion- 
al response  portion  of  the  four  models  using  the  M7  statistic  indicated 
that  they  were  homogeneous  and  a  combined  proportional  response  with  a 
single  a  and    c  could  be  substituted  for  the  four  separate  models. 

j      These    factorial    experiments    included    SO2    as    part   of   the  experimental 
design  and   lacked  true  replication.     The  error  term  used  in  scaling  the 
standard  errors  and  in  statistical   testing  was  the  residual  mean  square 
obtained   from  the  best  fitting  polynomial   model   of   form:   yield  =  Pq 
PlOa  +  I322S02  +  (33032  +  £. 

^     The  Ml  statistic  was  not  available  for  this  data  set. 


1 


These  three  data  sets  are  from  a  split-plot  design  with  the  three  culti- 
vars grown  in  the  same  treatment  plots. 
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Larsen,  R.  I.  and  W.  W.  Heck.  1976.  An  air  quality  data  analysis  system  for  interrelat- 
ing effects,  standards,  and  needed  source  reductions:  Part  3.  Vegetation 
injury.    Journal  of  the  Air  Pollution  Control  Association  26:  325-333. 


In  general  mathematical  terms,  the  relations  can  be  expressed  as 
z  =  -  In  mghr  /  1"        -    p  In  t  /  In  Sg    +    In  c  /  In  Sg 

where 

z       is  the  number  of  standard  deviations  that  the  percent  leaf  injury  is  away 
from  the  median,  (z  =  0  for  50%  injury), 

mghr       the  geometric  mean  concentration  of  the  gas  for  1  hour  of  exposure, 

s       is  the  standard  geometric  deviation  for  concentrations  producing  a  stand- 
ard deviation  of  percent  leaf  injury, 

p       is  the  slope  of  the  percent   leaf   injury   line  on   logarithmic  paper  for 
concentration  and  duration, 

t       is  the  exposure  duration  in  hours,  and 

c       is  the  pollutant  concentration. 
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Loehman,  E.  and  T.  Wilkinson.  1983.  Ozone  damage  to  field  crops  in  Indiana.  Station 
Bulletin  Number  426,  Agricultural  Experimental  Station,  Purdue  University,  West 
Lafayette,  Indiana.    38  pp. 

Regression  Equations  for  Percent  Yield  Reduction  From  Ozone 

Soybean,  cv.  Corsoy: 

Z  =  619.2  X  -  15.5 

cv.  Hodgson: 
Z  =  414.7  X  -  10.4 

cv.  Davis: 
Z  =  349.9  X  -  8.7 

cv.  Williams: 
Z  =  266.6  X  -  6.6 

cv.  Essex: 
Z  =  133.4  X  -  3.3 

Z  =  percent  soybean  yield  reduction, 

X  =  average  7  hour  daily  maximum  ozone  concentration  above  0.025  ppm. 


Field  Corn:  Z  = 
Winter  Wheat:  Z  = 


1533.7  X2  -  144.6  X  +  26.6;    X  >  0.069  ppm 

2077.8  X2  -  73.8  X  +  0.54 
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Medeiros,  W.H.,  P.O.  Moskowitz,  E.A.  Coveney,  and  H.C.  Thode,  Jr.  1983.  Oxi- 
dants and  acid  precipitation:  a  method  for  identifying  and  modeling 
effects  on  United  States'  soybean  yield.  Paper  presented  at  76th 
Annual  Meeting,  Air  Pollution  Control  Association,  Atlanta,  Georgia. 
1983,  June  19-24. 


Summary  of  regression  results.^ 


Modelb 


Model  Multiple 
Significance  r2 


1 .  Univariate 


a.  Y  =  29.8  -16580. 3[H+] 

Std.  error:  0.24  6665.7 
P:  <0.01  0.01 


0.01 


0.003 


b.  Y  =  40.8  -2.1[0XID] 

Std.  error:  1.17  0.21 
P:  <0.01  <0.01 


<0.01 


0.05 


Bi variate 


a.  Y  =  40.9 

Std.  error:  1.17 
P:  <0.01 


-9110. 3[H+] 
6553.1 
0.17 


-2.0[0XID] 

0.21 
<0.01 


<G.01 


0.05 


Bivariate  with  Maturity  Zone  Indicator 


a.  Y  =  31 .61  +Zi  -0.82[0XID] 

Std.  error:  1.20  (0.38-0.79)  0.23 
P:  <0.01  (<0. 00-0. 51)  <0.01 


-2589.60([H+]  [OXID]) 
1038.72 
0.01 


<0.01 


0.49 


b.  Y  =  32.24  +Z1     -16565. 5[H+]  -0.93[0XID] 

Std.  error:  1.23  (0.38-0.79)  5550.09  0.23 
P:  0.01     (0.01-0.59)    <0.01  <0.01 


<0.01 


0.49 


Y  =  bushels/acre,  H"*"  =  moles/litre,  OXID  =  7  hour  seasonal  average  (pphm), 
and  Z19  =  zone  intercept  adjustments. 


The  number  of  residual  mean  square  degrees  of  freedom  used  in  determining 
the  significance  of  the  regression  model  coefficients  is  1871,  minus  the 
number  of  coefficients  estimated  in  the  model.     The  number  of^  degrees  of 


freedom  is  based  on  the  number  of    independent  residuals 


where  Y-j 


is  the  model  estimate  of  Y-j .  Although  the  model  estimates  are  also  based 
on  estimated  data  (and  are  therefore  dependent),  this  does  not  affect  the 
independencies  of  the  residuals. 
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The  following  is  a  synopsis  of  the  numerical   regression  models  obtained  after 
least  square  fitting.    The  meaning  of  the  symbols  is  as  follows: 

y    =  specific  plant  response  parameter 
x-|  =  number  of  pollutant  episodes 

X2  =  cumulative  integral  of  the  pollutant  episodes  (concentration  and  duration) 
X3  =  peak  pollutant  concentration 

A.  Fumigation  Type:    Open-Top  Chamber  Dispensed  SO2  Treatment 


1.  Plant  response:    number  of  pods  per  plant 

y  =  290  -  18. 7x^  -  0.0043x2  -  O.OlSx^  -1-  0.354x^2 

2.  Plant  response:    seed  weight  per  plant 

y  =  82.6  -  5.80x^  +  O.OOlBx^  +  O.OHx^  ^  0.106x^2 

3.  Plant  response:    number  of  seeds  per  plant 

y  =  626  -  39. 3x^  -  0.012X2  "  O-O^SOx^  -1-0.753x^2 

4.  Plant  response:    thousand  seed  weight 

y  =  -174  +  13. 4x^  -1-0.047x2  -•-  0.497x2  -  0.326x^2 

B.  Fumigation  Type:    Open-Top  Chamber  Dispensed  Oa  Treatment 


1.  Plant  response:    number  of  pods  per  plant 

y  =  0.142  +  1.42X1  -(-  0.0002x2  -  0.273x3  -  5.64  sinx2 

2.  Plant  response:    seed  weight  per  plant 

y  =  28.0  -  0.161x1  +  O.OOOIX2  -  2.05  sinx2  -  1.79  sinx3 

3.  Plant  response:    number  of  seeds  per  plant 

y  =  43.3  -t-  1.45X1  +  O.OOOIX2  -  0.289x3  -  10.5  sinx2 

4.  Plant  response:    thousand  seed  weight 

y  =  331  -  1.68x1  -  0.0001X2  "  ^-^^  sinx2  -  16.6  sinx3 

C.  Fumigation  Type:     Residual-Ambient.   Unfiltered.  and  Filtered  Open-Top  Chamber 

Treatment 


1.  Plant  response:    pods  per  plant 

y  =  438  -  8.30x1  "  0.154x2  -  0.213x3  -  0.0154x5 

2.  Plant  response:    seed  weight  per  plant 

y  =  180  -  3.36X1   -  0.0863x2  "  0.0847x3  -  O.OII6X5 

3.  Plant  response:    number  of  seeds  per  plant 

y  =  35.5  -  0.0069X2  -  0.0421x3  -1-  O.OOIOX5  -  0.037x6 

4.  Plant  response:    thousand  seed  weight 

y  =  1601  -  26.9X1  "  0.581X2  "  O.63IX3  -  O.OOI6X5 


The  meaning  of  the  independent  variables  is  as  follows 

xi  =  number  of  episodes  of  SO2 

X2  =  integral  of  the  SO2  concentration 

X3  =  peak  of  the  SO2  concentration 

X4  =  number  of  episodes  of  O3 

X5  =  integral  of  the  Oa  concentration 

X5  =  peak  of  the  Oa  concentration. 
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Nosal,  M.  1984.  Atmosphere-biosphere  interface:  analytical  design  and  a  computerized 
regression  model  for  lodgepole  pine  response  to  chronic  atmospheric  SO2 
exposure.  RMD  Report  83/26  and  83/27,  Research  Management  Division,  Alberta 
Environment,  Edmonton,  Alberta.  194  pp. 


3  3 
y    =  Z    E  a:j(n)  x^n  +  j;    i  (bT(ni)sin  w^x^  +  cA^)  cos  t^^x^) 

i=l  n  i=l  m 

where 

y    =  annual  basal  area  increment, 
x-|  =  number  of  SO2  episodes, 

X2  =  cumulative  integral  of  SO2  concentrations  with  respect  to  time, 
X3  =  peak  SO2  concentration,  and 

a.^"^),  b.('^)  =  model  parameters 


Oshima,  R.J.,  M.P.  Poe,  P.K.  Braegelmann,  D.W.  Balding,  and    V.  Van  Way.  1976. 

Ozone  dosage-crop  loss  function  for  alfalfa:  a  standardized  method 
for  assessing  crop  losses  from  air  pollutants.  Journal  of  the  Air 
Pollution  Control  Association  26(9):  861-865. 


Summary  of  the  regression  of  Moapa  69  alfalfa  yield  (Y)  and  leaf  to  harvest 
weight  ratios  (L)  with  seasonal  ozone  dosage  (D),  average  daily  maximum 
temperature  (TH),  average  daily  minimum  temperature  (TL),  and  average  daily 
relative  humidity  (RH). 


Y  =  162.4  -  (0.01503  x  D)  L  =  0.6398  -  (1.941x10-5  x  D) 

df  df 


Matrix  Reg.  Error  r  Matrix  Reg.  Error  r 

Y  X  D           1               7  0.827**  L  x  D  1  7  0.890** 

Y  X  TH         1               7  0.644  L  x  TH  1  7  0.373 

Y  X  TL         1               7  0.536  L  x  TL  1  7  0.089 

Y  X  RH         1              6  0.692  L  x  RH  1  6  0.610 

df  t-values  of  regression  coefficients 

Matrix  Reg.  Error  0  TH            TL  RH  F 


Y 

X 

D 

X 

TH 

2 

6 

-3.949** 

-2, 

,374* 

15. 

.449** 

Y 

X 

D 

X 

TL 

2 

6 

-5.364** 

3, 

.062* 

21  . 

.381** 

Y 

X 

D 

X 

RH 

2 

5 

-3.695* 

0.085 

15. 

.419** 

Y 

X 

D 

X 

TH 

X 

TL 

3 

5 

-4.645** 

-0. 

,882 

1  . 

.624 

13. 

.987** 

Y 

X 

D 

X 

TH 

X 

RH 

3 

4 

-3.129* 

-0. 

,720 

-0.594 

9. 

.463* 

Y 

X 

D 

X 

TL 

X 

RH 

3 

4 

-4.162** 

1  , 

.640 

-0.447 

14 

.654* 

Y 

X 

D 

X 

TH 

X 

TL  X 

RH  4 

3 

-3.447* 

-0, 

,549 

1  . 

.375 

-0.680 

9. 

.148* 

L 

X 

D 

X 

TH 

2 

6 

-4.458** 

-0, 

.444 

12. 

.032** 

L 

X 

D 

X 

TL 

2 

6 

-4.785** 

-0. 

.091 

11  , 

.573** 

L 

X 

D 

X 

RH 

2 

5 

-3.625* 

0.449 

11  , 

.949* 

L 

X 

D 

X 

TH 

X 

TL 

3 

5 

-4.051** 

-0. 

.642 

-0. 

.497 

7, 

.097* 

L 

X 

D 

X 

TH 

X 

RH 

3 

4 

-3.371* 

0. 

.283 

0.550 

6, 

.956* 

L 

X 

D 

X 

TL 

X 

RH 

3 

4 

-3.500* 

-0, 

.159 

-0.744 

7, 

.441* 

Significant  at  0.05  level. 
Significant  at  0.01  level. 
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Calculated  regression  equations  for  significant  interactions  of  Moapa  69  alfalfa  yield 
(Y)  and  leaf  to  harvest  weight  ratios  (L)  with  ozone  dosage  (D),  average  daily  maximum 
temperature  (TH),  average  daily  minimum  temperature  (TL)  and  average  daily  relative 
humidity  (RH). 


Matrix 


Equation 


Y  X  D  X  TH 

Y  X  D  X  TL 

Y  X  D  X  RH 

Y  X  D  X  TH  X  TL 

Y  X  D  X  TH  X  RH 

Y  X  D  X  TL  X  RH 

Y  X  0  X  TH  X  TL  X 


RH 


Y=395.0889-( .0125  x  D)-(2.1864  x  TH) 

Y^=-166.3028-(. 01406  x  D)-i-(5.8919  x  TL) 

Y-168.0667-( .0179  x  D)+(.1107  x  RH) 

Y=25.6715-( .0132  x  D)-(1.2566  x  TH)+(4.3117  x  TL) 

Y=544.3926-( .0167  x  D)-(3.457  x  TH)-(1.703  x  RH) 

Y=-39.7543-(  .0174  x  D)-i-(4.2767  x  TL)-(.5351  x  RH) 

Y=237.1657-(  .0166  x  D)-(2.4218  x  TH)-t-(4.0032  x  TL)-(1.7644  x  RH) 


L  X  D  X  TH 

L  X  D  X  TL 

L  X  D  X  RH 

L  X  D  X  TH  X  TL 

L  X  D  X  TH  X  RH 

L  X  D  X  TL  X  RH 


L=.6976-(l  .8787x10"^  x  D) -( 7  . 001  7x1 0'''  x  TH) 

L=.6549-(l .9455x10"^  x  D)-(2 . 7253x1 O""  x  TL) 

L-.6893-(2. 3155x10"^  x  D) -( 7  . 730x1 0"'*  x  RH) 

L=.8788-(l .8455x10"^  x  D)-(l .4653x1 O"^  x  TH) -( 2 . 1 1 50x1 O"^  x  TL) 

L=.3014-(2. 4361x10"^  x  D)+(l .0964x10"^  x  TH)+( 3 . 5631 xl O"^  x  RH) 

L=.8398-(2. 3472x10"^  x  D) -( 3 . 0520x1 0""*  x  TL)-(3. 0976x10"^  x  RH) 


Temperatures  were  calculated  in  degrees  Fahrenheit. 

Equations  with  negative  intercepts  do  not  produce  negative  yields  because  the  mini 
mum  temperature  values  were  never  lower  than  52°F. 


Rowe,  R.  D.  and  L.  G.  Chestnut.  1985.  Economic  assessment  of  the  effects  of 
air  pollution  on  agricultural  crops  in  the  San  Joaquin  Valley  (SJV). 
Journal  of  the  Air  Pollution  Association  35:728-734. 


Crop 
(Variety) 

Equation 

I.    SJV  field  data  regression  results^ 

Dry  beans 

Y  = 

1.3255  -  0.00802  (OsAVE).  N=87,  DF=73,  r2=0.7682 
(-2.01) 

Y  = 

1.1043  -  0.000871  (03GE10),  N=87 ,  DF=73,  r2=0.7710 
(-2.24) 

Y  = 

1.044  ton/acre  OaAVE  =  69.2  OsGElO  =  69.12 

K^U  I,  UUll 

Y  - 
I  — 

n  ci'^/m  -  n  nn^^iQ  ^noA\/F^    m-i?    rf-^q    R2~n  oiiq 

u.joHj       u.uu^jj   \U3nv  L. )  ,   n—  1  c  f   ur  —  jj,   r\  — u.oi/j 
(-1.39) 

Y  = 

0.4624  -  0.000319  (OaGElO),  N=72,  DF=59,  r2=0.8209 
(-1.73) 

Y  = 

0.440  ton/acre  OaAVE  =  36.5  OaGElO  =  70.1 

Grapes 

Y  = 

7.629  -  0.00000166  (03GE6)2,  N=72,  DF=59 , r2=0 . 8209 
(-1.72) 

Y  = 

7.628  tons/acre  O3GE6  =  494.7 

Potatoes 

YP  = 

=  1.0872  -  0.00185  (O3GEIO)  -  0 . 0068( SO2GEI 0) , 
(-1.17)  (-1.30) 

N=26,  DF=15,  r2=0.7613 

YP  = 

=  1.004  O3GEIO  =  26.9  SO2GEIO  =  4.9 

II.    NCLAN  regression 

results'^ 

Corn  (pooled) 

Y  = 

12277.3  -  301727  (x2) 

Wheat  (pooled) 

Y  - 

4852.5  -  169239  (x^) 

Grain  sorghum 
(Dekalb,  A28+) 

Y  = 

8157.6  -  72388.3  (x^) 

Lettuce  (Empire) 

Y  = 

1065  -  5978  (X) 

Tomatoes 
(Murietta) 

Y  = 

67.1   -  38.85  (SO2)  -  1703  (x^) 

III.    Brewer  and  Aschroft^ 

Alfalfa    %A  Y  -  11.5  -  0.00677  (O3GEIO)  -  0.00133  (SO2GEIO), 

(-4.6)  (-2.3) 
N-11,  DF=8,  R2-0.796,  F=10.43 
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IV.  Legend 


Y  =  yield  per  acre 

YP  =  yield  adjusted  by  national  productivity  index 

O3AVE  =  average  hourly  O3  per  month  summed  over  growing  season  in  pphm 

O36EIO  =  number  growing  season  daytime  hours  where  Oa  >  10  pphm 

O3GE6  =  number  growing  season  daytime  hours  where  O3  >  6  pphm 

SO2GEIO  =  number  growing  season  daytime  hours  where  SO2  >  10  pphm 

X  =  average  daytime  hourly  O3  over  the  growing  season  in  ppm 

SO2  =  average  daytime  SO2  over  the  growing  season 


^     Rowe   and   Chestnut   T-ratios   reported   in  parenthesis   under  air  pollution 
coefficients.    One-tailed  statistical  tests  are  appropriate. 

b     Heck  et  al.  (1984) 

^     Rowe  and  Chestnut  regressions  on  Brewer  and  Ashcroft  results. 


Stevens,  T.  H.  and  T.  W.  Hazelton.  1976.  Sulfur  dioxide  pollution  and  crop 
damage  in  the  four  corners  region:  a  simulation  analysis.  New  Mexico 
Agricultural  Experimental  Station  Bulletin,  647.  Las  Cruces,  New 
Mexico. 

The  effect  is  described  mathematically  by: 
Xs  =  f(Rs.  (C  -  Co)) 
Vi  ,s  =  ^i  ,s  ~  t)i  ,s^s 

where 

R       is  the  species  resistance  factor  for  the  air  pollutant, 

C       is    ambient    pollutant    concentration    producing    visible  foliar 
injury, 

Cq     is  the  theoretical  threshold  concentration  that  can  be  continu- 
ally tolerated  without  foliar  injury  occurring. 


X 

is 

percentage  leaf  area  injury. 

a 

is 

the  assumed  potential  percentage  crop  yield, 

b 

is 

a  yield-leaf  destruction  coefficient. 

y 

is 

the  percentage  crop  yield  under  polluted  air 

i 

is 

growth  stage,  and 

s 

is 

crop  plant  species. 

This  simple  model  is  limited  in  application  because  it  does  not 
include  processes  and  effects  of  soil  water  availability,  air  temperature, 
humidity,  and  solar  radiation. 
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8.2  EQUATIONS  FOR  MECHANISTIC  PROCESS  MODELS 

Coughenour,  M.B.    1981.    Relationship  of  SO2  dry  deposition  to  a  grassland    sulfur  cycle. 
Ecological  Modelling  13:1-16. 


The  resistances  of  diffusion  to  each  of  the  various  sinks  are  implemented  to 
yield  an  absolute  flux  occurring  over  the  solution  interval  of  the  simulation  (1  day). 
The  relationship 

C 

F  =  - 
r 

is  used  to  facilitate  this  calculation,  where  F  is  flux  rate  (yg  SO2  cm"^  s~^),  C  is  con- 
centration (yg  SO2  cm"3)  and  r  is  total  resistance  from  the  point  at  which  C  was  measured 
to  the  sink  (s  cm~^).  When  C  is  measured  at  Zp,  rg  must  be  included  in  r,  but  when  C  is 
measured  in  the  canopy  (as  is  the  case  here),       need  not  be  included  in  r. 

Flux  to  leaf  surfaces  or  interiors  is 

C 

F  =  -  LAI 
r 

where  LAI  is  leaf  area  index. 

Passive  deposition  onto  leaf  surfaces  is  governed  by  the  sink  strength  of  the 
surface.  This  sink  strength  is  represented  here  as  a  constant  leaf  surface  resistance 
(rju).  The  value  of  r^i^  was  found  by  trial  and  error  to  most  closely  match  observed 
and  predicted  values  of  sulphur  accumulation  on  dead  leaves.  The  total  resistance  for 
deposition  onto  live  or  dead  leaf  surfaces  is  then 

r  =  ra  +  rb  +  rsu 

In  this  model,  flux  to  the  soil  is  found  in  a  similar  fashion  by  use  of 
r  =  ra  +  re  +  r^g 

where  r^  and  are  the  atmospheric  resistances  to  diffusion  above  and  below  the 
canopy,  and  r^g  expresses  the  surface  resistance  of  the  ground. 

Over  the  logarithmic  profile  of  region  1,  the  resistance  to  diffusion  is 

ra  =  Up  /  U*2 

where  ra  is  aerodynamic  resistance. 

Friction  velocity  (U*)  is  then  estimated  with  the  stability-corrected  loga- 
rithmic windspeed  profile  relationship 

Up  1  r  (zp  -  d)]  _ 

_  =  _    In    (1  -  oRi)-l/2 

U*     k  L         ZQ  J 

where  Zp  is  reference  height  (cm),  a  is  an  empirical  constant  (=10)  and  Up  is  reference 
height  windspeed.  Zero  plane  displacement  d  =  0.63h  where  h  is  the  height  (cm)  of  the 
plant  canopy,  and  roughness  length  (zq)  follows  as 

zo  =  k  (h  -  d) 

where  k  is  Von  Karman's  constant  (0.41). 


conti  nued . 
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Because  the  windspeed   profile   is   affected   by  atmospheric   stability,   a  mean 
Richardson  Number  is  introduced 


gz  (T  -  Tq) 


where  g  is  the  gravitational  constant  (cm^  s~^),  T  is  reference  height  temperature  (k), 
Tq  is  the  temperature  of  the  surface  losing  heat  (canopy  temperature),  Up  is  reference 
height  windspeed  (cm  s~^),  and  z  is  height. 

/l  \l/2 
rb  =  0.12(  _)  D-2/3 

\V 

where  D  is  the  diffusivity  of  the  gas  and  Uq  is  canopy  windspeed.  The  characteristic 
length  (1)  is 

4A 

1  =  - 
P 

where  A  is  planar  area  (cm^)  of  the  leaf  and  P  the  perimeter  (cm).  A  logarithmic 
profile  is  resumed  in  region  3  and  the  atmospheric  resistance  from  Zc  to  zero  ground 
displacement  may  be  approximated  by 

rc  =  U(zc)  /  [U*(Zc)2] 
Eddy  velocity  at  is 

U(Zc)  =  Uh  e  -a(l-Zc/h) 
where  a  is  defined  by  the  relationship 

/h-d\\ 

Iw  =  -    h  /  In   


The  windspeed  at  canopy  height  h  is  then 

U*     /h  -  d\        _  _i/2 

Uh  =  —  Inl   (1-aRi) 

k       \  zo  / 

An  estimate  of  the  centre  of  the  canopy  layer  is 
Iw 

Zc  =  h  -  _ 
2 

The  depth  of  the  canopy  eddy  is  1^  =  0.473  (h  -  d)  and  extends  down  from  approximately 
h.    Within  this  layer,  friction  velocity  varies  with  height;  therefore,  at  z^ 

U*(Zc)  =  U*e-  a(l-Zc/h) 


SCq  =  SCrt  •  I(phen) 

^10  =  F-10  •  SC  •  Srt/Sroot 

where  SCq  is  the    determined  shoot  S/C    ratio,  SCrj  is  the  actual  root  S/C  ratio,  I(phen) 
is  the    predefined  ratio  of    shoot  S/C  to  root  S/C  as  a  function  of  phenology,  f]o  is  the 
translocation    rate  from  one    root  layer  (g  S  m"^  day"^),    F-|o  is  the  daily    shoot  carbon 
growth    increment  (g  C  m'^  day~i),  Srj  is  root  sulphur  in  a  layer,  and  Srqqt  "^^^^ 
phur  summed  over  all  layers  (g  S  m"^). 

continued . . . 
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If  NSs  <  35 

then  y  =  MIN  (1 ,  3.5  -  0.1  *  NSs) 
otherwise  y  =  0 

where  y  is  the  sulphur  limiting  effect  on  shoot  photosynthesis  and  growth. 

For  sulphur-conserving  shoot  death, 

CSs  =  (CNs  +  5)  •  NSs 

fll  =  (Fn/CSs)  -  Fii   •  cj/CSsTS 

CSsM  =  (d-cj)  •     Fll)  /  U^ 

where  CSs  ratio   of   senescing   shoots   after  conservation,  CNs 

ratio  of  these  shoots,   6  is  the  increase  in  the  C/N  ratio  in  tissues  where  protein 

hydrolysis    has   occurred,   a  is   the   fraction   of   structural    shoot  carbon,  NSs 

N/S  ratio  of  the  shoots,  Fn  is  the  total  rate  of  flow  of  carbon  due  to  shoot  death  (g 

C    m~2    day~^),    fn    is    total    metabolic    sulphur    flow    in    dying    shoots    (g    S  m"^ 

day    "1),    CSsTS    is    the    shoot    structural    C/S    ratio    (fixed),    and    CS^s    is    the  C/S 

ratio  of  the  dying  metabolic  components  of  dying  shoots.     Thus,  structural  sulphur  in 

dying  shoots  is  given  by 

fl2  =  (Fn/CSs)  -  fll 

and  total  remobi lization  of  sulphur  by 

fl3  =  Fll   •  •  ((1/CSmos)  -  (1/CSms)) 

which  is  distributed  among  roots  in  all  layers;  CSf/iQs  is  the  C/S  ratio  of  the 
metabolic  component  before  protein  hydrolysis. 

Root  sulphate  uptake  is  calculated 

U  =  Vmax  •  [S04]/(Km  +  [SOa]) 

fg  =  Cr  •  Er(T)   •  Er(W)   •  Eu(CSp)   •  U 

where  U  is  the  basic  Michael  is-Menten  uptake  rate  (g  S  m"^  day"^  g~^  root 
carbon),  V^ax  is  the  maximum  uptake  rate,  Kfj,  the  half -saturation  constant,  [SO4] 
is  the  concentration  of  sulphate  in  solution  (ppm  S),  fg  is  the  rate  of  uptake  (g  S 
m"2  day~^),  Cr  is  root  carbon  content  (g  C  m"^),  Er(W)  and  Er(T)  are 
the  effects  of  water  and  temperature  on  root  metabolism,  and  Eu(CSp)  is  the  effect 
of  the  whole-plant  S/C  ratio  on  uptake. 

Sulphate  adsorption-desorption  by  soil  is  formulated  by  the  Freundlich  isotherm 

AS  =  b  •  [SOnll/n 

where  [SO4]  is  the  concentration  of  sulphate  solution  (ppm  S),  AS  is  the  quantity 
adsorbed  (g  S  g"^  soil),  and  b  and  n  are  empirical  constants. 

The    solubility   product    for   the    ions    Ca^^   and   SO42-   is   used   to  determine 
gypsum  precipitation 

Ksp  =  [Ca2+]  •  [S04^-]  •  y2 

where  Ksp  is  the  solubility  product  (2.4  x  10"^  at  25°C;  and  Y  is  the  mean  activity  co- 
efficient. When  Ksp  is  reached,  gypsum  may  precipitate.  Although  Ksp  is  actually 
temperature-dependent,  this  feature  is  not  included  in  the  model.  The  concentration  of 
gypsum  can  be  estimated 

[CaS04]  =  [Ca2+]  •  [504^-]  •  y/Kp 

continued . . . 
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where  Kp  is  the  dissociation  constant  (4.9  x  10  ^j.  jhe  mean  activity  coefficient 
for  the  soil  solution  is  estimated  from 

In  Y  =  0.509  z2  v~  /  (1  +  v^) 

where  z  is  the  valence  of  the  ion  and  u  is  its  ionic  strength,  given  by 

n 

u  =  1/2    I  C-jZi 
i=l 

where  C-j  is  the  concentration  of  one  of  the  n  ions  in  solution.  The  ions  considered 
here  include  Na+,  K+,  Ca2+,  Mg2+,  and  SO^^-. 

Diffusion  of  sulphate  is  described  by  Pick's  law,  where  the  diffusion  rate  is 
proportional  to  the  concentration  gradient.  Thus, 

Df  =  D  •  (A  [SO4]  /AX)  •  e 

where  Df  is  the  diffusion  rate  (g  S  cm"^  day"^),  A  [SO4]  is  the  concentration  difference 
between  layers,  A  X  is  the  distance  (cm)  between  the  midpoints  of  adjacent  layers,  D 
is  diffusivity  (cm^  day~^),  and  e  is  the  volumetric  soil  water  content,  neces- 
sary here  to  express  diffusion  rate  per  cross-sectional  area  of  soil. 

The  rate  of  sulphate  uptake  is  then  modified  by  the  effects  of  soil  water 
content  and  temperature  on  general  microbial  metabolism,  and  by  the  demand  for  sulphur 
as  expressed  in  the  microbial  C/S  ratio 

33  =  Em(T)  •  Em(W)   '  £^(05^,,) 

and 

^3     ^3  *  ^m  ■  ^max 

[SO4]  /  (Km  +  [SO4]) 

where  a3  is  the  combined  effect  of  temperature  Em(T),  water  Efn(W),  and  C/S  ratio 
(Ey(CSrp)).  Vfpax  is  the  maximum  uptake  rate,  [SO4]  is  the  concentration  of  sulphate  in 
soil  solution  (ppm  S),  0^  is  the  microbial  biomass  (g  C  m"^)^  ^^id  K^,  is  the 
half-saturation  constant. 

In  the  model,  mineralization  and  sulphate  uptake  are  affected  by  the  same 
functions  of  temperature  and  moisture.  The  model  also  includes  a  microbial  C/S  effect, 
that  mirrors  the  effect  of  the  C/S  effect  on  sulphate  uptake,  given  by 

^4  -  Em(T)  •  Em(W)   •  Emn(CSm) 

f4      a4  •  Sff,  •  V^ax 

where  Em(T)  and  Em(W)  are  as  in  Eq.  [4]  (Coughenour  1981),  E^nCCSm)  is  the  effect  of  the 
C/S  ratio,  a4  is  the  combination  of  these  effects,  f4  is  the  mineralization  rate  (g  S  m"^ 
day~^),  Vmax  'i^  '^^^  maximum  uptake  rate  (day~^),  and  S^  is  microbial  sulphur  content 
(g  S  m~2).  Maximum  rates  were  chosen  as  0.8  day~^  for  bacteria  and  0.165  day"^  for  fungi. 

The  oxidation  of  sulphide  to  sulphate  is  controlled  by  soil  aeration  and 
bacterial  metabolism.  Bacterial  metabolism  is  expressed  as  a  function  of  temperature 
and  soil  water  content.    The  oxidation  relationship  is 

35  -  Em(T)   •  Em(W)   •  Eox(W) 

^5  "  35  •  Vox  •  Sp 

where  Em(T)  and  Em(W)  are  the  effects  of  temperature  and  water  on  the  bacteria, 
Eox(W)  is  the  effect  of  water  on  aeration,  35  is  the  combined  effect,  fs  the 
oxidation  rate  (g  S  m^  day"^),  Sp  is  the  quantity  of    reduced  forms  of  sulphur  (g  S  m'"^), 
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and    Vox    ""^  maximum   oxidation    rate    (day"^)    based   on    results    from   a   silt  loam 

soil  in  which  55%  of  the  reduced  sulphur  was  oxidized  in  three  weeks.  Thus,  this 
expression  assumes  that  the  rate  of  oxidation  is  proportional  to  S^. 

The  rate  of  volatilization  of  sulphur  during  decomposition  is  represented  as  a 
proportion  of  the  rate  of  decomposition  of  metabolic  litter,  increasing  slightly  as 
aeration  declines 

H  =  H  ■  Pv(W)  •  SCme 

where  ff^  is  the  volatilization  rate  (g  S  m"^  day~i),  F5  is  the  rate  of  metabolic  litter 
decomposition  (g  C  m"^  day~^),  P\/(W)  is  the  variable  proportion  of  F5  as  a  function  of 
soil  water  based  on  data  of  Banwart  and  Bremner  (1976)  as  cited  by  Coughenour  (1981),  and 
SCfne  is  the  S/C  ratio  of  metabolic  litter.  We  assume  that  this  equation  applies  to 
all  soil  layers. 

Sulphate  esters  are  not  explicitly  treated,  because  of  the  uncertainty  of  their 
origin,  but  are  assumed  to  be  associated  with  humads  and  resistant  soil  organic  matter. 
A  regression  equation,  based  on  the  data  of  Cooper  (1972)  and  Bettany  et  al.  (1973, 
1974)  as  cited  by  Coughenour  (1981),  predicts  that  Hl-sulphur  (S^,  g  S  m"2)  is  a 
function    of    total     sulphur    (Stotal'    9    ^  total    carbon    (Ctotal»    9  ^ 

m-2) 

Shi  =  0-71  •  Stotal  -  0-0029  •  Ctotal 

This  equation  can  be  used  to  calculate  HI  sulphur  at  any  point  in  the  simulation. 

When  microbial  death  is  calculated  in  the  decomposer  submodel,  sulphur  flow  is 
controlled  by  carbon  flow.  Thus 

f"!  =  Fi  /  CSstm 

where  f]  and    F-|  are  the  flows    (g  day"^-),    respectively,    of  structural    sulphur  and 

structural  carbon  due  to  death,  and  CSs-j-m  ''s  the  fixed  C/S  ratio  of  the  microbial 
structural  component. 

The  flow  from  metabolic  components  varies  with  the  C/S  ratios  of  the  microbes, 

as 

f2  -  (F2  /  CSm)  -  fi 

where  f2  and  F2  are  the  flows  of  metabolic  sulphur  and  metabolic  carbon  due  to 
death,  and  CSf^  is  the  C/S  ratio  of  the  microbes. 

f7  =  RS  •  p  •  a 

fe  -  Rs  •  P  •  (1  -  a) 

where  f-j  is  the  rate  of  sulphate  release  from  humads  (g  S  m"^  day"^),  and  fQ 
is  the  rate  of  release  from  resistant  organic  material  (g  S  m"^  day"^-). 

a  =  \  •  Sh  /  (X  •  Sh  +  Sre) 

a  =  SCf^/SCpg 

where  X  is  the  lability  of  humad-associated  esters  compared  to  that  of  resistant 
organic-associated  esters,  SCh  and  SC^e  ai^e,  respectively,  the  S/C  ratios  of  humads  and 
resistant  organic  material,  a  is  the  proportion  of  sulphate  originating  from  humads,  S^ 
is  humad  sulphur  content  (g  S  m"2),  Spg  is  resistant  organic  sulphur  (g  S  m~^) . 

If  Bq  is  the  sum  of  root  and  microbial  biomasses  per  cm^  of  soil  in  the 
surface  layer,  and  B]  is  the  equivalent  measure  in  a  deeper  layer,  then 

13  =  Bi  /  Bo 

is  the  proportional  decrease  in  activity  for  layer  1. 
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Flow  of  S  uptake  by  nematodes: 
fni  =  Ic/46 

Flow  of  S  excreted  by  nematodes  and  to  metabolic  litter 

fno  ml  =  0.25  E/46 
Flow  of  S  excreted  by  nematodes  and  to  structural  litter 

fno  si  =  0.75  E/46 
The  first  approximation  to  intake  is 

Ic  =  Mr(T)  •  0.02925 

where        is    the  intake  rate    (g  C  nematode    C  day~^),    Mr(T)    is  the  temperature- 

dependent  metabolic    rate  (ml  O2        fwt  h""i),  and  feces    production  £  =  (1-AE).  AE 

is  assimilation  efficiency  (0.5). 
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Symbols 

Description 

Vd  1  Uc   U  jcU 

In  Examples 

^L 

LUiivcibiuii  TdLuur  TUT  (. rid iiy  1  iiij  lhx   lu   icdT  ury  wciynu  \,^'>»  ) 

Jo 

31 

LUciiiLiciiL  ucuciiii  III  Illy    uraiioiULaLiuii  ui  piiuuuoyiiLiiaLc 

1 

Coefficient  determining  amount  of  storage  to  be  held 
1 n  rese  rve  ^  y  j 

1 

Ao 

Total  concentration  of  PGA  and  RuDP  (mol/L) 

see  k3 

A 

Ratio  of  total  root-surface  area  to  area  of  ground 

bl 

Coefficient  determining  translocation  of  photosynthate 

1 

b4 

Parameter  determining  slope  of  function  that  cuts  off  flow 
of  photosynthate  out  of  main  storage  as    main  storage  level 
approaches  a  critical  value 

0.1 

Be 

Converts  g  CO2  cm"^  s"^  to  g  dry  wt  m"^  da~^ 

4.69  X  108 

ci 

Reduction  in  photosynthesis  for  unit  increase  in  pollutant 
with  one  unit  of  diffusion  resistance  (s  cm~^  ppm~^) 

0.7125 

Inverse  slope  of  function  cutting  off  flow  of  photosynthate 
from  leaf  storage  to  structural  leaf  as  x2/xl  -»0 

0.05 

c 

Concentration  of  CO2  in  air  (g/cm^) 

310  ppm 

Cc 

Concentration  of  CO2  at  chloroplast  (g/cm^) 

Cia 

Concentration  of  CO2  in  pore  (g/cm^) 

Cp 

Concentration  of  CO2  in  protoplasm  (g/cm^) 

Csoil 

Concentration  of  salts  in  soil  (mol/cm^) 

0 

d 

Inverse  of  square  of  standard  deviation  of  foliage 
density  (m"^) 

d7 
d8 

Thickness  of  rooting  zone  layer  in  soil  (m~^) 
Thickness  of  surface  zone  layer  of  soil  (m"^) 

D 

Diffusivity  in  the  canopy  (cm^/s) 

Dd 

Time  of  dawn  (da) 

Dl 

Leaf  width  (cm) 

7 

Do 

Turbulent  diffusivity  in  air  (cm^/s) 

300 

Ds 

Time  of  sunset  (da) 
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Symbols 

Description 

Value  Used 
In  Examples 

e 

Vapour  pressure  (mbar) 

es 

Saturation-vapour  pressure  (mbar) 

ei 

Leaf-growth  constant  (da~^) 

0.137 

64 

Growth  constant  for  stem  and  fruit  (da~i) 

2.44 

ee 

Growth  constant  for  root  (da""^) 

0.15 

ET 

Evapotranspi ration  (m/da) 

f(Z) 

Density  of  leaves  in  the  canopy  as  a  function  of  Z 

Angular  distribution  of  leaves 

F 

Leaf  area  index  LAI  (area  of  leaves/area  of  ground) 

Maximum  allowed  LAI  in  simulation 

5.6 

Ftot 

Total  LAI 

Net  flux  of  radiant  energy  out  of  layer  i  in  the  downward 
direction  (cal  cm"^  s"^) 

h' 

Ratio  of  Fjoj  to  total  canopy  depth  Zjgj  (cm~^) 

0.025 

H 

Unit  step  function  (Heaviside  function);  rate  of  sensible 
heat  transfer  from  the  leaf 

I 

Flux  of  light  or  light  intensity  normal  to  the  leaf  at  canopy 
depth  F  (nanoeinstens  cm"^  s"^  or  cal  cm"^  s'^) 

II 

Light  flux  reaching  F  from  the  direction  of  the  sun 

Im 

Light  intensity  at  photosynthesis  half  maximum 

In 

Net  radiation  into  the  canopy  (cal  cm"^  s"^) 

lo 

Light  flux  at  top  of  canopy  in  the  direction  of  the  sun 
(cal  cm~2  s"^) 

.02 

kl 

Equilibrium  constant  for  conversion  of  CO2  and  RuDP 
to  PGA  (cm3/g) 

2.28  X  108 

k2 

Equilibrium  constant  for  light  conversion  of  NADP  to  NADPH 
(cm2  s/cal  or  cm^  s/yeinsteins) 

21 .9  cm2 
s  cal-i 

k3 

Coefficient  for  maximum  photosynthesis  for  high  CO2 
and  high  light  (g  cm"^  s"^  mole'^Lz) 

k3AoNo  =  500 
mg  dm~2 

kd 

Relative  rate  of  growth  during  day 

1 

kn 

Relative  rate  of  growth  during  night 

1 

k 

Decay  constant  for  reduction  of  windspeed  in  canopy 

0.44 

•^soi  1 

Maximum  value  of  soil  hydraulic  conductivity  (weight  basis, 
m/da) 
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Symbols 

Description 

Vu  lUc:  UoCU 

In  Examples 

"^(Xl) 

nyuraUMC  LuriuucLiviLy  ui    trie  iuii   yiw/  ucx ) 

In 

CALinCtiUll    CUCll    ILICllL    UI       II^MU     III     tllC  UQllUpy 

L 

Niimhpr  nf  l;i\/pr<;  in  rr^nonv 

10 

"d 

Niimhpr  nf  nl;in1"^   npr  iini1"  Pirp;i  ^fn~2^ 

iiuiiiuci     v^i     piaiiL-j               uiiiu    uicu    V'li  / 

5 

n 

iiuriiiui    LU    Liic   pioiic  ui     Liic  leal 

Total  concentration  of  NADP  and  NADPH  (mol/L) 

see  k3 

p 

Phni"  n<;  vyn1"  hp   i  <;   (n  rm~2  c~iN 

rilWL-l./Oj^llUIICjl^                              Lrlll                       O  J 

UI           |J  1  lu  uu  o  y  1 1  L  M  c  ^  1  J           L.in        j  ) 

f^m 

1  lO Aiiiiuiii   vaiuco    UI    piiuuujyiiuiicoio   au   iii^ii    ii^iiu  icvcrio 

Pn 

Net  photosynthesis 

Pr 

Precipitation  (m/da) 

Pi 

Precipitation  intercepted  by  leaves  (m/da) 

rrcLi|jiLauiuii   oUi  ii^iii^   ^ruuiiu  wiuiiuuu    iiiuciLcpuiuii  \iii/ua^ 

PL 

Avpt^r^np  rn<;inp  nf   nrnipr+inn  nf   <;iin1inh1"  nntn   Ippif   fnr  ;i  1 1 
FAv^iuvj^    vuoiiic    \j  I     piujwL-uiwii    \j  I     ouiiii^iiv    \jiiwv.f     icai     I  \j  i  ail 

leaves 

PLa 

Average  cosine  of  projection  of  sunlight  onto  leaves  at  angle  a 

^10 

rvatiu   UI    rcb(JiraLiuii   at    \  '      ^    iu^    (\    LU   icbpiraLlUil   au  Ij^ 

°K,  where       is  a  reference  temperature 

2 

•"a 

Leaf-boundary-layer  resistance  (s/cm) 

Carboxylation  resistance 

Turbulent  diffusion  resistance  from  canopy  layer  i  to  top 
of  canopy  (s/cm) 

Cell-wall-diffusion  resistance  (s'/cm) 

0.126 

Specific  growth-respiration  rate 

0.35  at  250c 

Stomatal  resistance  due  to  all  factors  except  water 
stress  (s/cm) 

Mesophyl 1 -eel  1  resistance  (s/cm) 

Specific  maintenance-respiration  rate 

0.006  at  250c 

•"P 

Protoplasm  resistance  to  CO2  diffusion  (s/cm) 

1  .0 

i^s 

Stomatal-dif fusion  resistance  (s/cm) 

^soil 

Resistance  to  hydraulic  flow  of  water  in  the  soil  (da,  weight 
basis) 
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Desc  ri  pt i  on 

Value  Used 
In  Examples 

Resistance  to  flow  of  water  between  soil  and  the  plant 
(da,  weight  basis) 

Parameter  proportional  to  slope  of  curve  of  stomatal  resistance 
as  a  function  of  water  stress  (soil-water  potential) 
(bar  s  cm"i) 

5.004 

Ti  lO 

Value  of  »|»5  at  which  stomatal  resistance  is  singular  (bar) 

-15 

Soil-water  potential  at  a  soil-water  content  of  0  (bar) 

-267 

Logarithm  of  the  ratio  of  the  soil-water  potentials  with  one 
unit  of  soil-water-content  difference  between  them  (m~^) 

85.67 

Critical  or  maximum  stomatal  resistance  attained  when 
is  maximum 

•^1 

Leaf-boundary-layer  resistance  when  Dl  =  1  cm  and 
U  =  1  cm/s  (s/cm) 

1  .3 

Respiration  (usually  leaf)  (g  cm"*^  s~^) 

Rr 

Gas  constant  (bars  ^K"-*-  mol~^  cm^,  or  cal 
mol-i  oK-i) 

8.314  J 
mol-i  °K-i 

Rx 

Dark  leaf  respiration  at  25°C  (g  cm~2  s"^) 

2.5  mg 
(jni-2  h"^- 

s 

Unit  vector  in  direction  of  the  sun 

s 

Absorption  coefficient  of  CO2  in  water 

0.756  at  2501 

s 

Air-pollutant  concentration  (ppm) 

Sn 

Threshold  for  depressant  effect  of  pollutant  (ppm) 

20  pphm 

L 

Ratio  of  leaf  storage  to  leaf  structural  tissue  at  which 
photosynthesis  is  at  half  maximum 

Sp 

Reserve  concentration  of  leaf  storage  photosynthate  that  the 
leaf  retains 

0.1 

t 

Time  (da) 

T 

Temperature,  usually  of  air  (oK) 

Ti 

Transfer  flux  to  compartment  i  from  storage  (g  m"^  d"^) 

'  r 

Root  or  soil  temperature  (°K) 

Transfer  coefficient  for  movement  between  leaf  storage  and 
main  storage  (da~^-) 

3.0 

Tx 

Reference  temperature  for  calculating  respiration  (°K) 

250c 

u 

Inverse  of  exponent  that  transforms  soil  content  to  soil 
potential 

Inverse  of  max[Pn( I) ]/Ca  at  high  light  levels  (ppm  g"^  cm^  s) 
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Symbo 1 s 

uesc  r 1  PL  1  on 

Value  Used 
in  LXamp 1 es 

U 

Windspeed  (cm/s) 

V 

txponeriL  or  soii  pouenxiai  useu  xo  caicuiaLe  soil  conuucuivity 
from  soil  potential 

VI 

Inverse  of  slope  of  PpCI)  at  low  light  levels  (nanoeinsteins/g) 

w 

Value  of  soil  potential  at  which  soil  conductivity  has  fallen  to 

llu  1  1     lllaA  MMUIIl 

W-| 

^1 

Rnnm^cc    I          w+ ^    nf        riir  +  11  r"3  1    IpAf   m;?  +  p  r  i  ;i  1    I  n\ 

DIUIMu^^     ^Uiy    WU^     Ul      oLlULuUlul      ICal     lliuLCI    lul  Vy^ 

y  r> 

Diuiiiu^o    ^uiy  WLy    ui     icui    iiiuuci  lul     III    oLUiayc  Vy^ 

^3 

D  1  uiiid ib   ^ury  wu^  ui   iiia  in  buurayc  ^y^ 

X4 

OLclll   UlUMIabo    ^Uiy   w  I,  J  \\iJ 

y  r 

r  1  u  1  u   u  1  uiiiu jj    ^uiy  wuy  vy/ 

Dnn+    Knnmacc     1  i\         1.1+  \     /  n  ^ 

Koot  Diomass  v^ry  wt,^  vy/ 

X7 

Soil  water  in  the  root  zone  (m) 

Soil  water  in  the  soil-surface  zone  (m) 

xg 

Water  in  the  leaf-interception  storage  (m) 

Maximum  water  content  of  the  root  zone  (m) 

0.12 

Z 

Vertical  distance  usually  measured  from  the  top  of  the 

Canopy  ^ITi; 

Thickness  of  the  canopy  (m) 

a 

p 

Angle  of  the  leaf  to  the  horizon 
Elevation  of  the  sun 

T 

Psychometric  constant  (mbar/oK) 

0.66 

0 

Art  nip    Hpi'ufPPn  +Kpciinan/4+Kplpaf 
Hliy  1  c    Uc  LWccd    trie    bull    dllU    Lllc  Ical 

A 

ucrivaLive  or  saxurauion  vapour  pressure  wiun  respecx.  xo 
temperature  (mbar/°K) 

n 
e 

Angle  between  the  sun  and  the  normal  to  the  leaf 

Volumetric  water  content  at  which  matrix  potential  is 
equivalent  to  a  pressure  of  1  m  of  water 

\ 

Latent  heat  of  vaporization  (cal/g) 

558 

y 

Power  of  X4  used  to  calculate  biomass  of  X4  capable  of  growth 

1  .0 

Coefficient  proportional  to  slope  of  rs(I)  (s^  cm  cal"^) 

1.3  X  10-3 
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Symbol s 

Description 

Value  Used 
In  Examples 

no 

Coefficient  inversely  proportional  to  square  root  of  slope 
of  rs(I)  at  I  =  0  (nanoeinsteins  cm"^        or  cal 
cm~2  s~i) 

3.33  X  15-^ 
cal  cm"2 

Slope  of  rs(T)  (s  cm~i 

0.5 

at  0.0°C,  infinite  light  levels,  and  no  water  stress  (s/cm) 

-7.0 

V 

Water-stress  component  to  stomatal  resistance 

El 

Relative  growth  constant  for  stem 

2.0 

lo 

Relative  growth  constant  for  fruit 

0.0346 

PI  <: 

Saturation  concentration  of  water  vapour  at  the  leaf 

PI 

Concentration  of  water  vapour  at  reference  height  Z 

pc 

Volume  heat  capacity  of  air  (cal  cm"^  °K~i) 

290  X  10& 

a 

Function  of  a  and  |3  used  in  calculating  PL^^ 

TO 

Physiological  time  for  initiation  of  fruit  growth  (degree-da) 

800 

Physiological  time  for  termination  of  fruit  growth  (degree-da) 

2500 

Critical  water  potential  at  which       attains  its  maximum 
value  (m,  weight  basis) 

I's 

Soil-water  potential  (bars) 

'J'Li 

Leaf  potential  for  layer  i  (m,  weight  basis) 

^1 

Soil-water  potential  for  root-zone  layer  (m,  weight  basis) 

'I'B 

Soil-water  potential  for  soil-surface  layer  (m,  weight  basis) 

Azimuthal  angle  between  sun  and  leaf 

continued . 
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We  define  growth  as  the  transfer  of  photosynthate  from  a  storage  compartment 
(X2  leaf  storage  or  X3  main  storage)  to  a  structural  compartment  (X]  leaf,  X4 
stem,  X5  fruit,  or  X5  root). 

The  equations  describing  the  rate  of  change  of  the  state  variables  x-|  through 
X5  with  time  are  described  mathematically  as: 


and 


where 


dx-| 
dt 

dX2  SlPtH(Pt) 


PtH(Pt) 

+  H(-Pt)P7  -  [1  +  rg(T)]Ti  -  T3  -  rn(T)X2, 


dt  Sl  +  X2/X1 
dx3 


dt 


T3  -  (l+rg)(T4+T5+T5)  -  rn(Tp)X6  -  rn(T) ( X3+X4+X5) 


dx4 

  =  T4  -  rn(T)X4  H(a4X4iJ-X3) , 

dt 

dx5 

  =  T5  -  rn(T)X5  H(a4X4W-x3), 

dt 


dX5 
dt 


T6  -  <"n(Tr)X5  H(a4X4V-x3) 


F 

Pt  =  Be  J  Pp(F)  dF/Md. 
0 

rg(T)      =  growth-respiration  rate  with  a  dependence  on  T  as  R  above, 
rn(T)      =  maintenance-respiration  rate  with  a  dependence  on  T  as  R  above, 
T  =  temperature  of  above-ground  parts  (°K), 

Tp  =  temperature  of  roots  (°K), 

y  =  exponent  of  X4  for  calculating  the  fraction  of  X4  used  in  growth, 

a4  =  empirical  parameter, 

$1  =  ratio  of  leaf  storage  to  leaf  structural  tissue  at  which  photosynthe- 

sis is  at  half  its  maximum  value, 
Pj  =  added  biomass  (g  dry  wt  da~^  plant"^). 

Be  =  conversion    factor  from  g  CO2    cm"^  s"^  to  g  dry    wt  m'^  da"^    (4.69  x 

108),  and 

=  transfer  flux  to  compartment  i,  i=l,3,4,5,6. 
Tl  =  eiXi{kn[H(Dd-t)  +  H(T-Ds)] 

X2/X1 

+  kdH(t-Dd)  H(Ds-t)}    (1  -  MdXi/aLFo), 


=  Trns  X2 


aiX2 

X]  +•  X2      X3  +  X4  +  X5  +  X5, 
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X3  -  a,x^' 


b4  ^  X3  -  a^x^' 


and 


T5  =  H(t  -  tq)  e^^^x^Hix^  -  a^x^^)  H(t^  -  t) 


b4  ^  X3  -  a^x^ 


-  H(X3  -  a4X4  ) 


Pp 


S  >  Sq, 


,  0  <  s  <  s. 


where 


C]      =  reduction  in  photosynthesis  for  a  unit  increase  in  pollutant  at  one  unit 

of  diffusion  resistance,  (s  cm~^  ppm~^), 
rs     =  stomatal  resistance  (s/cm), 
S       =  concentration  of  pollutant  (ppm),  and 
rb     =  stomatal  resistance  at  which  Sq  was  determined. 


-b  -Vb2  -  4ac 


-  R 


2a 


where 

a  =  (k2l  +  1)  k-i[rp  +  r^^  ^  s(rs+ra)]. 

b  =  -(k2l  +  1)  {sCaki  ^-  Rki[rcw  +  s(rs+ra)]  +  1} 

-  k-ik2k3AoNoI[rp  +  r^^  +  s(rs+ra)]  ,  and 

c  =  kik2k3AoNoI{R[rcw  +  s(rs+ra)]  +  sCa)  . 
Following  Sinclair,  we  use  for       and  r^: 
ra  =  r-i  (Dl/U)^/2^ 

where  Dl  =  leaf  width,  U  =  windspeed,  and  r-\  is  a  constant,  and 
VI 
I  +  P2 


^s 


y3T  +•  v4  + 


where  y] ,  y2«    v3»         ^4  constants    determined  by    fitting  rs(I,T)  to  data  obtained 

by  varying  light  intensity  and  leaf  temperature  while  measuring  rs  •  yi  is  in  units  of 
either  nanoei nstei ns  cm"^  or  cal  cm"^,  ^2  in  either  nanoei nstei ns  cm"^  5"^  or  cal  cm~2 
s"^,  P3  in  s  cm"i  C~^,  and  y4  in  s  cm~^.  v  is  a  water  stress  term  to  be  discussed  below. 
The  biochemical  constants,  k-j ,  k2,  k3,  Aq,  and  Nq  can  be  easily  determined  by  examining 
Eq.  [21]  in  Kercher  1977  at  various  limits  and  using  experimental  data  to  fit  the  equa- 
tion at  these  limits  (see  Sinclair  in  Kercher  1977  for  a  full  discussion).  Our  expression 
for  respiration  is  that  of  Waggoner: 
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R  =  Rxexp[Tx(Tx  +  10)  lnQio(1/Tx  -  1/T)/10] 

where 

Rx        =  base  respiration  rate, 

Tx        =  base  temperature  (°K),  and 

OlO       =  ratio  of  R(Tx  +  10)  to  R(Tx)  . 

dX7      i\iQ  -  i(»7  -  iifLi 

  =    H(X7max  -  X7)  -  E    , 

dt        rsoil  i  rsp 

dX8 


dt 


[P2  +  PiH(X9  -  xgmax)]  H(X8max  "  ^s) 
^8  ~  '♦'7 

-    -  ET(G,O.T,e.O,rgnd)  . 

•'soil 

dxg  L 

  =  PiH(X9niax  -  xg)  -  E  ET(  Ini  ,Ini+i  .T.e.O.ra^ '  +  r^')  . 

dt  canopy 

layers 

i  =  1 

«(f8       =  soil-water  potential  of  surface-soil  water  (m), 

v|;7       =  soil-water  potential  of  rooting-zone  water  (m), 

ij/L-j      =  leaf-water  potential  of  leaves  in  layer  i  (m), 

•"soil  =  resistance  of  soil  to  water  flow  (da), 

r^p     =  resistance  to  flow  between  soil  and  plant  leaves  (da), 

P2       =  precipitation  striking  ground  without  interception  (m/da),  and 

P-|       =  precipitation  intercepted  by  leaves  (m/da). 


't'7  -  fLi 

  =  ET  (Ini,  Ini+i,  T,  E.  r^^' ,  rai '  +  Tci')  . 

•^sp 

Stomatal  resistance  (of  H2O)  of  the  i'^h  foliage  layer  is  given  by 

('"c  ~  '•'Li)'" 

where 

ro     =  critical  or  maximum  resistance  attained  when  vi»|_-j  =  tc(s/cm) 

=  critical  potential  (mbar),'and 
m       =  slope  of  resistance  with  respect  to  change  in  ^\_^  at  (fnbar"^) 


ci  DoK^h 

da  d7 


and 


•"soil 


sp 


2K(X8)  2K(X7) 
d7         d7  6 


2Kr  fK(x7)A 
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where 


d-j  =  thickness  of  soil  layer  i  (m), 

6  =  ratio  of  cross-sectional  area  of  roots  in  X7  to  area  of  X7, 

f  =  fraction  of  root  surface  area  through  which  water  is  absorbed, 

A  =  ratio  of  total  root  surface  area  to  area  of  ground, 

Kp  =  conductivity  of  plant  tissue  (m/da); 

and  we  calculate  soil  potentials  using 


and 


Soil  conductivities  are  calculated,  using 


K(xi)  = 


where  v  and  u  are  empirical  constants  easily  determined  by  log-log  fits  of  soil  water 
potential  versus  soil  water  content  and 

l^soil      =  soil  conductivity  at  soil  potential  equal  to  0, 

e  =  volumetric  water  content  at  which  matric  potential  is  equivalent  to  a 

pressure  of  1  m  of  water,  and 
w  =  value  of  soil  potential  at  which  soil  conductivity  has  fallen  to  half 

maximum. 


dX7  L 

  =  Pr  H(x7max  "  X7)  -  E  ETdpi  ,Gi  ,Ti  ,e,rsi  •  ,rai '-t-rci ' ) , 

dt  i=l 


where 


Ip-j       =  total  incident  radiant  energy  into  layer  i  (cal  cm'^  s~^), 

G-j         =  total  incident  radiant  energy  emerging  from  layer  i  (cal  cm"^  s~^), 

T-j         =  temperature  of  layer  i, 

e  =■  vapour  pressure  of  atmosphere  above  the  canopy  (mbar), 

rs-j '      =  stomatal  resistance  of  layer  i, 

rg-j '      =  boundary  layer  resistance  of  leaves  in  layer  i, 

r^i  '      =  canopy-diffusion  resistance  from  layer  i  to  a  reference  level  where  e 

is  measured, 
Pp        =  precipitation  rate  (m/da),  and 
H(x)      =  1 ,    X  >  0 
=0,    X  <  0. 


rs'  =  n'  +  v'. 
Curry  gives  v'  as  a  function  of  soil-water  potential 


i-wl  i^wl 

v'    =    -   

'^w2  +  1's(x7.Cs)  +  ts(x7max.O) 

where  the  soil-water  potential  is  a  function  of  soil-water  content  and  soil  salt 
concentration. 
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where 


where 


Rq       =  gas  constant  (83.14  x  10^  mbar  cm^  mol-^-  , 

Tp       =  temperature  of  soil  (°K),  and 

^soil  =  concentration  of  salt  in  soil  (mol  cm"^  solute). 

Adn  -  G)  +  pc[es(T)  -  e]  /  (rg-  +  rc ' ) 

ET  =    864  , 

\{A  +  Y  [1  +  rs'/(ra'  +  rc')]} 


ET  =  transpiration  (m/da), 

es(T)  =  saturation  vapour  pressure  at  temperature  T  (mbar), 

pc  =  volume  heat  capacity  of  dry  air  (290  x  10"^  cal  cm~3 

Y  =  psychometric  constant  (0.66  mbar  °K~^), 

\  =  latent  heat  of  evaporation  (585  cal/g), 

A  =  des/dT; 


and  es(T)  and  A  are  calculated  by  the  computer  code  using  the  Clausius-Clapeyron  equation. 
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Kercher,  J.R.  and  M.C.  Axelrod.  1981.  SILVA:  a  model  for  forecasting  the  effects  of 
SO2  pollution  on  growth  and  succession  in  a  western  coniferous  forest. 
UCRL-53109,  Lawrence  Livermore  National  Laboratory,  Livermore,  CA.  72  pp. 


H  =  137  +  b2D  -  530^, 

where  b2  and  b3  are  constants. 

G  D  [1  -  (DH)/(DmHfn)] 

5D  =    r(A)  Q(DEGD)  S(BAR)  W(^)  GR 

274  +  3b2D  -  4b3D2 

where: 


a 

growth  constant  for  each  species  at  the  optimum 

L 

leaf  area  of  tree  (cm^) 

D 

diameter  of  tree  (dbh  in  cm) 

H 

height  of  tree  (cm) 

maximum  height  species  attains 

Dm 

maximum  dbh  species  attains 

r 

response  to  light  (dimensionless) 

A 

available  light  at  tree 

Q 

response  to  climate  (dimensionless) 

DEGD 

growing  degree  days  (F-days) 

S 

response  to  soil  quality  (dimensionless) 

BAR 

total  basal  area  of  stand  (cm^) 

response  to  water  stress  (dimensionless) 

T 

ratio  of  actual  to  potential  evapotranspi ration 

GR 

response  to  pollution 

A 

=  Aq  e(-Kl  LAI) 

i 


where   k   is   an   extinction   coefficient  for  coniferous   canopies,   LAI   is   the   leaf  area 
index,  and  the  sum  is  over  all  trees  taller  than  the  tree  under  consideration,    Aq  is 
the  light  at  full  sunlight.    The  constant  k  is  taken  to 
be  0.47. 

ZLAI  =  (ZL)/PLTSIZ 
where  PLTSIZ  is  the  ground  area  of  the  plot. 

r(A)  =  1  -  e[-4-64  (A-0.05)]^    shade  tolerant 

=  2.24{1  -  e[-1  •''3&(A-0.0B)]},    shade  intolerant. 
The  leaf  area  of  each  tree  as  a  function  of  tree  diameter  is  given  by: 

L  =  c.D°'-' 

where  c-j  and  ai  are  constants  determined  for  each  species  for  logarithmic 
regressions.    Usually  a-j  lies  between  1.5  and  2.5. 

Q  =  4  (DEGD  -  DMIN)  (DMAX  -  DEGD)/(DMAX  -  DMIN)2 

DEGD  =  J  (T  -  40)  B(T  -  40)  dt 
growing 
season 

where  T  is  the  temperature  in  Fahrenheit,  and  B(T  -  40)  is  1  for  T  greater  than  40  and  0 
for  T  less  than  40. 

S(BAR)  =  1  -  BAR/(S0ILQ  •  10.000  •  PLTSIZ) 


continued . 


103 


where  BAR  =  total  basal  area  of  all  species  on  stand.  The  parameter  SOILQ  is  the 
maximum  basal  area  which  the  plot  can  carry  and  is  dimensionless  (m^/m^).  This 
parameter  is  determined  by  soil  nutrients.  The  factor  of  10,000  converts  cm^  of  BAR 
to  m2.  The  function  S  reduces  growth  for  all  trees  as  the  total  basal  area 
approaches  the  maximum  allowed. 


WSM  -  T 


NWS 


W(t)  =  1  - 


^WSM  -  WSO  / 


to  quantify  W,  where 


T  =  annual  actual  evapotranspi ration  /  annual  potential  evapotranspi ration 

WSM  =  value  of  t  at  which  W  is  maximum 

WSO  =  lower  limit  of  range  of  tolerance  in  t 

NWS  =  exponent  relating  linear  function  of  t  to  tree  response. 

Average  overall  M  age  classes  to  get  the  total  response  of  the  tree: 


GR 


M 

^  Z  w^GRk 
k=l 


where  w^  =  weight  of  each  needle  age  class  to  total  productivity. 

Let   GRPk   (SAj)    be   the   response  of  the  productive  capacity  of  the  k'^'^  year 

needles  to  fumigations  during  the  j'^h  year.     The  current  year  is  1;  last  year  is  2; 

etc.  Then  the  accumulated  effect  of  all  fumigations  on  the  k'^'^  year  needles  is  given 
by  either: 


GRk  =  GRPkl  Z  SAj 


or 


GRk 


GRPk(SAj) 


SA  =  seasonal  average  concentration  (ppm). 

GRPk(SAj)  =  1  -  PL  (SAj  -  CR) 

CR  =  threshold  concentration  for  growth  reduction  (ppm). 
PL  =  fraction  growth  potential  lost  per  ppm  of  SO2  (ppm~^). 

PL  and  CR  vary  from  species  to  species. 

For  a  series  of  successive  episodes  during  a  growing  season,  it  is  necessary  to 
assume  that  the  concentration  y-j  and  duration  t-j  for  each  episode  are  known.  The 
accumulated  effect  of  all  fumigations  on  the  k'^h-year  needles  is  given  by 


GRk  =  GRPk^Z  J  yj  dt 

or 


GRk 


GRPklZ_  yij  Atij 


Equation  31a  or  31b  is  used 
response  by  averaging  over  all 


in  Equation  19  (Kercher 
M  needle-year  classes. 


1977)   to  get  the  total  growth 
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The    response   of    the    productive   capacity   of    the    k'^^^-year  needles   to   the  N 
episodes  during  the  j'^h  year  is  given  by 

/n    q     \  1       N  q 

GRPk(  I  Jyjdt)  =1-<C      -  In  zjy  (t;^! j      j )dt  -  In  mg  /  In  Sg 

\i      /         q  i 

or  ^  ^ 

N    q  1      N  q 

GRPk  z  y-jjAt-jj  =^-<l>    -In  I  y-jjAt-jj  -In  rrig  /In  Sg 

T  q  i 

where  in  Equation  30b  (Kercher  and  Axelrod  1981), 

At-jj  =  the  duration  of  the  i'^'^  episode  in  the  j'^h  year,  and 

y^j    =  the  concentration  of  SO2  in  the  i'^^  episode  in  the  j'^h  year,  depending  on 
the  pollutant  episode  data  set. 

The  <^  refers  to  the  fractional    loss  of  the  productivity  of  the  plant  in  the 
presence  of  pollution  compared  to  productivity  in  the  absence  of  pollution. 

mg,  Sg,  and  q  are  parameters  determined  by  fitting  data. 

The  parameter  v  is   a  measure  of  the  skewness   of  the  pulse.     For  v  near  1, 
the  pulse  is  maximally  skewed;  as  v  increases,  the  pulse  becomes  more  symmetric. 

The  parameter  oj  is  the  width  of  the  pulse. 

Only  establish  seedlings  for  species  i  if  the  condition 

DMIN(i)  <  DEGD  <  DMAX(i) 

is  satisfied.    Total  number  of  trees  of  each  species  J  added  to  the  plot  is 

N(J)  =  SDLINGS(J)   •  Ps- 

J 

P  =   


J(l  -  h)  +  N 

where 

P  =  probability  of  a  successful  seed  production  year, 

J  =  number  of  successful  seed  production  years  in  a  sequence  of  N  years, 
h  =  1   +  number  of  blocked  years  following  a  successful  seed  crop  year  for  the 
species . 

10  seedlings  established  per  (10m  x  10m)  plot  (SDLINGS)  for  each  species  for  its  good 
seed  year. 

As  a  first  order  approximation  for  the  effects  of  shade  on  seedling  survival  of 
species  intolerant  to  heat  and  moisture  stress, 

P3(I)  =  probability  of  survival  (moisture-stress  intolerant) 
=  1  _  e(-Yl  LAI). 

On  the  other  hand,  some  species  are  adapted  to  exposed  sites,  are  shade  intolerant,  and 
exhibit  increased  mortality  under  shade  as  seedlings.    For  these  species 

Ps(T)  =  probability  of  survival  (moisture-stress  tolerant) 
.  e(-Y2  LAI). 

A  quantitative  fit  for  y-]  and  y2  was  not  possible  from  the 
available  data  in  the  above  references.  The  authors  have  hypothesized  the  value  1.0  for 
both  Y")  and  y2  because  of  the  qualitative  descriptions  in  the  above  references. 
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PDj  =  PDe  +  PDq  -  PDeRDq 

where  the  T,  E,  and  G  subscripts  refer  to  total,  ecological,  and  growth-suppressed, 
respectively,  and 

PDj  =  probability  of  single  tree  death  during  year 
PDq  =  0.369 
PDe  =  4/AGEMX. 


King,  D.A.,  J.R.  Kercher,  and  G.E.  Bingham.  1983.  Modeling  the  effects  of  air  pollutants 
on  soybean  yield.  In:  Analysis  of  Ecological  Systems:  State-of-the-art  in 
Ecological  Modelling,  eds.  W.K.  Lauenroth,  G.V.  Skogerboe,  and  M.  Plug.  New 
York:    Elsevier  Scientific  Publishing,  pp.  545-555. 


The  change  in  photosynthetic  capacity  F,  caused  by  a  pollutant  acting  on  that 
component  with  turnover  time  d"^-. 

df 

_  =  B(I,E)  -  (d  +  I)F 
dt 

where 

F  =  factor  multiplying  leaf  photosynthesis 

I  =  internal  concentration  of  phytoactive  agent 

8  =  regeneration  rate  of  F 

E  =  available  energy  for  regeneration. 

The  function  B  and  rate  d  are  hard  to  measure.     For  small  turnover  and  repair 
rates  and  chronic  fumigations,  a  simple  approximate  solution  to  Eq.  1  is  given  by 

t 

F  =  1  -  a  ;  Idt 
0 

where 

a  =  constant  determined  from  photosynthesis  measurements, 
I  =  internal  concentration  of  phytoactive  agent. 

Assuming  that  internal  loss  of  toxic  agent  is  proportional  to  its  concentration 
and  the  agent  is  replenished  by  diffusion  into  the  leaf,  then  Eq.  2  may  be  rewritten 

t  C03(t)] 

F  =  1  -  a'  J    dt 

0  r(t)  +  b 

where 

r        =  sum  of  resistance  to  O3  diffusion 

[O3]    =  external  ozone  concentration  (time  dependent) 

a';  b  =  constants  determined  from  photosynthesis  measurements. 
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where 


Mortensen,  P.    1984.    Modelling    ion    uptake  in    agricultural    crops.    RISO-M-2446.  Riso 
National  Laboratory,  Roskilde,  Denmark.    33  pp. 

Calculation  of  the  total  concentration  of  an  ion  is: 
^in  ^out 

CC  =  CCo  +    (mol  •  kg(DW)-M 

BM 

where 

CCq  is  the  concentration  of  an  ion  at  the  beginning  of  calculation 
(mol -kgCDW)-!) ,  and 

BM  is  the  biomass  of  the  crop  (kg(DW)-m^2) 
Content  of  a  certain  ion  in  the  crop  is: 
Cin  =  Cino  +  lupt  '  (mol-m-^) 

lupt  =  Wupt  •  Csoil  •  Kup       (mol-m-2  •  time'^, 

lupt  'i^  "the  ion  uptake  in  the  crop  from  the  interstitial  water, 

C-jno  is  the  internal  content  of  an  ion  at  beginning  of  calculation 
(mol -m'^) , 

Wupt  "^^  "the  water  uptake  (mm  or  l(w)-m"2), 

Csoil  is  the  calculated  soil  water  concentration  of  a  certain  ion  as 
described  in  Broderson  1984  (mol • 1 (w) ) ,  and 

Kyp  is  the  uptake  coefficient  dependent  on  the  crop  type  and  ion  (no  unit). 

The   introduction   of   the   uptake   coefficient   is   for  calibration   purposes  and 
shall  be  looked  upon  as  a  selectivity  coefficient  with  the  following  description: 

Kyp  >  1      active  ion  uptake 

Kyp  =  1      passive  ion  uptake 

Kyp  <  1      selective  ion  uptake 

If  the  total  water  uptake  per  month  from  a  soil  is  W-^up,  then  the  uptake  of  water  from 
compartment  (i)  is: 

^upi  "  W-tup  •  Fpoi        (mm  •  month~i) 
Fraction  of  roots  in  soil  layer: 


FpQ-j  =  I  RDC       (decimal  fraction) 
li 

RDC  is  the  root  distribution  coefficient. 

Calculation    of    the    content    of    a    certain    ion   on   the   surface   of   a   crop  is 
calculated  in  the  following  way: 

Cout  =  Couto  +  (UPair  -  ^wof)  "  DT  (mol-m^^) 
DT  is  the  timestep  (e.g.,  month). 

continued . . . 
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This  wash  off  rate  (Rwof)  described  in  a  simple  way  as  a  constant  wash  off 
coefficient  (C^of)  multiplied  with  the  surface  intercepted  content  (Cgut)  "the 
vegetation: 


where 


where 


Rwof  =  Cwof  •  Cout       (mol-m-2.time  M 

C^of       wash  of  constant  (0.037  -  0.054  day-^)  (Chamberlain  1970)  (time"!),  and 

Cout  ^5        content  on  the  surface  of  intercepted  ions,  (mol-m"^). 

This  description  shall  be  used  in  the  present  model  as  the  direct  interception 
of  air  pollutants  as  dry  and  wet  deposition  on  the  vegetation.    The  formulation  shall  be: 

UPair  =  F  •  DEP       (mol -m-^-time-^) 
F  =  1  -  e-MyW       (no  unit) 

F  is  the  fraction  intercepted  on  the  vegetation, 

DEP  is  the  total  deposition  of  the  pollutant,  calculated  from  the  air  pollution 
model  (Hojerup  1984)  (mol -m'^ -time"^) , 

My  is  the  filtration  parameter  with  values  around  2.3  -  3.3  (Chamberlain  1970) 
(m2  •  kg~^),  and 

W  is  the  herbage  areal  density  (dry  weight  •  m"^) 

Calculation  of  the  root  and  above  ground  content  of  ions,  e.g.,  calcium  shall 
be  calculated  in  the  following  way: 

CAtot 

CApoot  =    (mol (Ca) -m'^) 

(BM/(DClca  •  ORB  •  ILAG))  +  1 

where 

CApoot       calcium  in  the  root  biomass  (mol (Ca) •m"^) , 

CAtot  ^5  "the  total  calcium  in  root  and  overground  biomass  (mol  (Ca) -m'^) , 
BM  is  the  overground  biomass  of  the  crop  (kg(DW) -m'^) , 

DClca  is  the  distribution  coefficient  between  root/top  for  calcium  in 
(mol(Ca) •kg(DW)~i  root  biomass/(mol (Ca) • kg( DW) above  ground  biomass). 

The  organic   root  biomass   (ORB)    in  dry  weight  for  a  crop   in  a   certain  soil 
compartment  shall  be  calculated  in  the  following  way: 

ORB  =  0.07  •  Froi   '  BM/ILAG  ( kg( DW) • 1 ( s ) ) 

BM  is  the  above  ground  biomass.  The  use  of  0.07  as  exchange  coefficient  for  all  crops 
is  a  rough  estimate.  This  shall  be  changed  when  better  data  is  available.  ILAG  is  the 
depth  of  the  soil  layer  (mm). 

Calculation  of  the  above  ground  calcium  content  CA^eg       calculated  as: 

CAveg  =  CAtot  -  CApoot       (mol (Ca) -m-^) 
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Distribution  coefficient  for  calcium  DClca  is  0.5  if  the  soil  water  concentration  is 
less  than  10"^  mol-L(w)~^  and  0.1  if  the  concentration  is  greater  than  or 
equal  to  10"=  mol-L(w)"^.  The  same  distribution  coefficients  are  used  for  all 
ions  except  potassium  (K)  where  a  value  of  0.5  is  used  for  all  concentration  levels 
(Baker  1983). 

For  cadmium  a  value  of  0.1   is  used  for  potatoes  and  0.2  for  beets  (Coughtrey 
and  Thorne  1983).    For  other  crops  the  values  for  calcium  are  used. 

CAveg 

CAstraw  =   (mol ( Ca) -m-^) 

(BMg/(DC2ca  '  BMst))  +  1 

where 

CAyeg  is  the  total  Ca  content  of  the  overground  vegetative  biomass 
(mol -m'^) , 

BMg  is  the  biomass  of  grain  (l<g( DW) -m'^) , 
BM3-t  is  the  biomass  of  straw  (kg(DW)-m~2)  and 

DC2ca  is  the  distribution  coefficient  between  straw  and  grain  for  calcium 
(mol(Ca-kg(DW)-i  straw/(mol (Ca) -kgCDW)"!  grain). 

Calcium  content  in  grain  is  then: 

CAgrain  =  CAyeg  -  CAstraw       (mol ( Ca) -m'^) 


Concentrations  are  found   by  division  with  the  grain  biomass.     It  is  assumed 
that  50%  of  the  above  ground  biomass  of  cereals  is  grain. 
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A  listing  of  the  program  of  the  model  SUCROS. 


TITLE  SUCROS  -  A  SIMPLE  AND  UNIVERSAL  CROP  GROWTH  SIMULATOR 

***  DRY  WEIGHT  OF  PLANT  ORGANS,  GROWTH  RATES  AND  PARTITIONING 


WLV  =  INTGRL  (WLVI,  GLV  -  DLV)  101 

WST  =  INTGRL  (0.,  GST)  102 

WSO  =  INTGRL  (0.,  GSO)  103 

WRT  =  INTGRL  (WRTI,  GRT)  104 

*  WEIGHTS  OF  LEAF  BLADES,  STEMS  (TRUE  STEMS  AND  LEAF  SHEATHS),  STORAGE 

*  ORGANS  AND  ROOTS  RESP.,  IN  KG/HA 

INCON  WLVI  =  25.,  WRTI  =  25.  105 

*  WEIGHT  OF  LEAVES  AND  ROOTS  AT  EMERGENCE 

GTW  =  (GPHOT  -  MAINT)  *  CVF  106 

*  GROWTH  RATE  OF  ALL  ORGANS  COMBINED,  IN  KG/HA/DAY 

GRT  =  GTW  *  (1 .  -  FSH)  107 

GSH  =  GTW  *  FSH  108 

GLV  =  GSH  *  FLV  109 

GST  =  GSH  *  FST  110 

GSO  =  GSH  *  FSO  111 

*  GROWTH  RATES  OF  ROOTS  AND  SHOOTS  (LEAVES,  STEMS,  STORAGE  ORGANS) 

*  IN  KG/HA/DAY 

DLV  =  WLV  *  RDR  112 

*  DEATH  RATE  OF  LEAVES,   IN  KG/HA/DAY 

RDR  =  AFGEN  (RDRTB,  DVS)  113 

FUNCTION  RDRTB  =  0.,  0.,  1.,  0.,  1.01,  0.03,  2.,  0.03  114 

WLVD  =  INTGRL  (0. ,  DLV)  -  .  115 

*  DEAD  MATERIAL  (LEAVES)  AT  THE  FIELD  IN  KG/HA 

FSH  =  AFGEN  (FSHTB,  DVS)  116 

FUNCTION  FSHTB  =0.,  0.5,  0.3,  0.5,  0.45,  0.775,  0.7,  0.825,  1.,  1.,  2.,  1.  117 

*  FRACTION  OF  GROWTH  OCCURRING  IN  SHOOTS  AS  FUNCTION  OF  DEVELOPMENT  STAGE 

FLV  =  AFGEN  (FLVTB,  DVS)  118 

FST  =  AFGEN  (FSTTB,  DVS)  119 

FSO  =  1  .  -  FLV  -  FST  120 

FUNCTION  FLVTB  =  0.,  1.,  0.45,  1.,  0.85,  0.,  2.,  0.  121 

FUNCTION  FSTTB  =  0.,  0.,  0.45,  0..  0.85,  1.,  1.,  1.,  1.01,  0.,  2.,  0.  122 

***  CARBON  BALANCE  PROCESSES 

LAI  =  WLV  *  SLFA  201 

PARAM  SLFA  =  0.0020  202 

*  LEAF  AREA  INDEX  IN  HA/HA  AND  SPECIFIC  LEAF  AREA 

*  IN  HA  LEAF/  KG  LEAF  WEIGHT 

GPHOT  =  DTGA  *  30.  /  44 .  203 

DTGA  =  FOV  *  DGAO  +  (1.  -  FOV)  *  DGAO  204 

*  GROSS  PHOTOSYNTHESIS  IN  KG  (CH20  AND  C02  RESP.)  PER  HA  PER  DAY, 

*  CALCULATED  FROM  LEAF  CHARACTERISTICS  (HMAX,  EFF),  LAI  AND  ACTUAL 

*  DAILY  RADIATION  (AVRAD),  AND  CORRECTED  FOR  DAYLENGTH  (DL  AND  DLE): 

DGAC  =  INSW  (LAI  -  5.,  PHCL,  PHCH)  205 

DGAO  =  INSW  (LAI  -  5.,  PHOL,  PHOH)  206 

PHCH  =  0.95  *  (PHCHl  +  PHCH2)  +  20.5  207 

PHCHl  =  SSLAE  *  AMAX  *  DLE  *  X  /  ( 1 .  +  X)  208 

X  =  ALOG  (1.  +  0.45  *  DRC  /  (DLE  *  3600.)  *  EFFE  /  (SSLAE  *  AMAX))  209 

PHCH2  -  (5.  -  SSLAE)  *  AMAX  *  DLE  *  Y  /  (1.  +  Y)  210 
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Y  =  ALOG  (1.  +  0.55  *  DRC  /  (OLE  *  3600.)  *  EFFE  /  ...  211 

((5.  -  SSLAE)  *  AMAX))  211 

SSLAE  =  SIN  ((90.  +  DEC  -  LAT)  *  PI  /  180.)  212 

PHCL  =  AMINl  (PHC3.  PHC4)  *  (1 .  -  EXP  ((AMAXl  (PHC3,  PHC4)  /  ...  213 

AMINl  (PHC3,  PHC4))))  213 

PHC3  =  PHCH  *  (1.  -  EXP  (-0.8  *  LAI))  214 

PHC4  =  DL  *  LAI  *  AMAX  215 

PHOH  =  0.9935  *  PHDHl  +  1 .1  216 

PHOHl  =  5.  *  AMAX  *  OLE  *  Z  /  (1.  +  Z)  217 

Z  =  DRO  /  (OLE  *  3600.)  *  EFFE  /  (5.  *  AMAX)  218 

PHOL  =  AMINl  (PH03,  PHC4)  *  (1.  -  EXP  (-  (AMAXl  (PH03,  PHC4)  /  ...  219 

AMINl  (PH03,  PHC4))))  219 

PH03  =  PHOH  *  (1.  -  EXP  (-0.8  *  LAI))  220 

EFFE  =  (1  .  -  REFLC)  *  EFF  221 

PARAM  EFF  =  0.5,  AMAX  =  3.0,  REFLC  =  .08  221 

*  INITIAL  LIGHT  USE  EFFICIENCY  AND  LIGHT  SATURATED  C02  ASSIMILATION 

*  RATE  OF  INDIVIDUAL  LEAVES.     UNITS:  KG  C02/  HA/  HR  /  (J/  M2/  S)  AND 

*  KG  C02/HA  LEAF/HR 

FOV  =  (DRC  -  AVRAD)  /  (0.8  *  DRC)  222 

*  AVERAGE  FRACTION  OF  PERIOD  OVERCASl  DURING  A  DAY. 

*  CALCULATION  OF  DAILY  RADIATION  OF  A  CLEAR  AND  AN  OVERCAST  SKY  (DRC 

*  AND  DRO,  P.A.R.,  IN  J/M2)  AND  OF  DAYLENGTH  (IN  HR)  AS  A  FUNCTION  OF 

*  LATITUDE  (LAT,  IN  DEGREE),  DECLINATION  (DEC,  IN  DEGREE)  AND  DATE: 

DRC  =  0.5  *  1300.  *  RDN  *  EXP  (-0.1  /  (RDN  /  (DL  *  3600.)))  223 

DRO  =  0.2  *  DRC  224 

RDN  =  3600.  *  (SINLD  *  DL  +  24.  /PI  *  COSLD  *  SQRT  (1.  -  ...  225 

(SINLO  /  COSLD)  **  2))  225 

SINLD  =  SIN  (DEC  *  PI  /  180.)  *  SIN  (LAT  *  PI  /  180.)  226 

COSLD  =  COS  (DEC  *  PI  /  180.)  *  COS  (LAT  *  PI  /  180.)  227 

DEC  =  -23.4  *  COS  (2.  *  PI  *  (DAY  +  10.)  /  365.)  228 

DL  =  12.  *  (PI  +  2.  *  ASIN  (SINLD/COSLD))  /  PI  229 

DLE  =  12.  *  (PI  +  2.  *  ASIN  ((-SIN  (8.  *  PI  /  180)  +  SINLD)  /  ...  230 

COSLD))  /  PI  230 

DLP  =  12.  *  (PI  +  2.  *  ASIN  ((-SIN  (-4.  *  PI  /  180.)  +  SINLD)  /  ...  231 

COSLD))  /  PI  231 

CONSTANT  PI  =  3.1416  232 

PARAM  LAT  =  -15.  233 

*  MAINTENANCE  RESPIRATION 

MAINT  =  AMINl  (GPHOT,  MAINTS  *  TEFF)  234 

MAINTS  =  WLV  *  0.03  +  WST  *  0.015  +  WSO  *  0.01  +  WRT  *  0.01  235 

TEFF  =  010  **  (0.1  *  TMPA  -  2.5)  236 

PARAM  QIO  =  2.  237 

*  GROWTH  EFFICIENCY 

CVF  =  (FLV  *  0.72  +  FST  *  0.69  +  FSO  *  CVFSO)  *  FSH  +  . . .  238 

(1.  -  FSH)  *  0.72  238 

PARAM  CVFSO  =0.73  239 

***  DEVELOPMENT  OF  THE  VEGETATION 

DVS  =  INTGRL  (0.,  INSW  (DVS  -  1.,  DVRV,  DVRR))  301 

FINISH  DVS  =  2.  302 

DVRV  =  0.0252  *  AFGEN  (DVRTTB,  TMPA)  *  AFGEN  (DVRDTB,  DLP)  303 

DVRR  =  0.0477  *  AFGEN  (DVRRTB,  TMPA)  304 

FUNCTION  DVRTTB  =  10.,   .63,  15.,   .83,  20.,   .92,  25.,   .96,  30.,   .98,   ...  305 

35.,   .99  305 

FUNCTION  DVRRTB  =  10.,   .08,  15.,   .38,  20.,   .575,  25.,   .71,  30.,   .80,   ...  306 

35.,   .865  306 

FUNCTION  DVRDTB  =  10.,  0.223,  11.,  0.425,  12.0,  0.575,  13.,  0.685,   ...  307 

14.,  0.767,  15.,  0.828,  16.,  0.872,  17.,  0.906  307 
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***  WEATHER  DATA 

DAY  =  AMOD  (TIME,  365.)  401 

AVRAD  =  0.5  *  41820.  *  AFGEN  (AVRADT,  DAY)  402 

FUNCTION  AVRADT  =  1.,  523.,  15.,  526.,  46.,  532.,  74.,  575.,  105.,  ...  403 

557.,  135.,  509.,  166.,  466.,  196.,  482.,  227.,  545.,  258   403 

611.,  288.,  644.,  319.,  556.,  349.,  521.,  365.,  523.  403 

TMPA  =  0.5  *  (AFGEN  (TMAXT,  DAY)  +  AFGEN  (TMINT,  DAY))  404 

FUNCTION  TMAXT  =  1.,  28.6,  15.,  28.3,  46.,  28.1.  74.,  28.7,  105.,   ...  405 

28.9,  135.,  27.2,  166.,  24.9,  196.,  24.7,  227.,  27.4,  258.,   ...  405 

29.9,  288.,  33.6,  319.,  31.3,  349.,  28.9,  365.,  28.6  405 

FUNCTION  TMINT  =  1.,  18.2,  15.,  18.2,  46.,  18.2,  74.,  16.3,  105.,   ...  406 

13.9,  135.,  10.1,  166.,  8.5,  196.,  7.4,  227.,  9.7,  258.,   ...  406 

13.4,  288.,  16.6,  319.,  18.2,  349.,  18.2.  365.,  18.2  406 

***  SIMULATION  RUN  SPECIFICATIONS 

TIMER  FINTIM  =  1000..  DELT  =  2.,  PRDEL  =  2.,  OUTDEL  =  2.,  TIME  =  300.  501 
*              INITIAL  VALUE  OF  TIME  INDICATES  STARTING  DAY  OF  SIMULATION 

METHOD  RKSFX  502 

PRINT  WLV.  WLVD,  WST,  WSO,  WRT,  LAI,  DVS,  MAINT,  DT6A,  CVF  503 

NWRT  =  -WRT  504 

WLVT  =  WLV  +  WLVD  505 

WVEG  =  WLVT  +  WST  506 

TADRW  =  WVEG  +  WSO  507 

PRTPLOT  NWRT,  WLV,  WLVT,  WVEG,  TADRW  508 

PAGE  GROUP,  NPLOT  =5  509 

END 

STOP 

ENDJOB 
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An  explanation  of  the  abbreviations  used  in  the  mpdel  SUCROS  as  listed. 


NAME  DESCRIPTION  UNIT 


AMAX         C02  ASSIMILATION  RATE  OF  A  LEAF  AT  LIGHT 

SATURATION  KG  (C02)/  HA  (LEAF)/  H 

AVRAD        ACTUAL  DAILY  RADIATION  (400-700  NM)  J/M2/D 
AVRADT      TABLE  MEASURED  GLOBAL  RADIATION  (CAL/  CM2/  D) 

VS.  DAYNUMBER 

CVF  CONVERSION  EFFICIENCY  FOR  GROWTH  OF  PLANT  DRY 

MATTER  KG  (DM)/  KG  (CH20) 

CVFSO        CONVERSION  EFFICIENCY  FOR  FORMATION  OF  STORAGE 

ORGANS  KG  (DM)/  KG  (CH20) 

DAY  NUMBER  OF  DAYS  IN  THE  YEAR  FROM  1ST  OF  JANUARY  DAY 

DEC  DECLINATION  OF  SUN  WITH  RESPECT  TO  THE  EQUATOR  DEGREE 

DGAC         DAILY  GROSS  C02  ASSIMILATION  -CLEAR  SKY-  KG  (C02)/  HA/  D 

DGAO         DAILY  GROSS  C02  ASSIMILATION  -OVERCAST  SKY-  KG  (C02)/  HA/  D 

DL  ASTRONOMICAL  DAYLENGTH  H 

DLE  EFFECTIVE  DAYLENGTH  H 

DLP  PHOTOPERIODIC  DAYLENGTH  H 

DLV  DEATH  RATE  OF  THE  LEAVES  KG/  HA/  D 

DRC  PHOTOSYNTHETICALLY  ACTIVE  RADIATION  -STANDARD 

CLEAR  SKY-  J/  M2/  D 

DRO  PHOTOSYNTHETICALLY  ACTIVE  RADIATION  -STANDARD 

OVERCAST  SKY-  J/  M2/  D 

DTGA         ACTUAL  DAILY  GROSS  C02  ASSIMILATION  KG  (C02)/  HA/  D 

DVRDTB      RELATION  BETWEEN  RATE  OF  DEVELOPMENT  AND 

DAYLENGTH 

DVRR         RATE  OF  DEVELOPMENT  IN  REPRODUCTIVE  PHASE  IN 

RELATION  TO  TEMPERATURE  1/D 

DVRRTB      TABLE  OF  DVRR  AS  FUNCTION  OF  TEMPERATURE 

DVRTTB      RELATION  BETWEEN  RATE  OF  DEVELOPMENT  AND 
TEMPERATURE 

DVRV         RATE  OF  DEVELOPMENT  IN  VEGETATIVE  PHASE  IN 

RELATION  TO  TEMPERATURE  AND  DAYLENGTH  1/D 

DVS  DEVELOPMENT  STAGE  OF  THE  CROP  FRACTION 

EFF  EFFICIENCY  OF  USE  OF  ABSORBED  VISIBLE 

RADIATION  FOR  C02  ASSIMILATION  AT  LOW  LIGHT 

LEVELS  KG  (C02)/  3/  HA/  H  M2S 

EFFE         EFF  BASED  ON  INCIDENT  RADIATION  KG  (C02)/  J/  HA/  H  M2S 

FLV  FRACTION  OF  LEAVES  IN  SHOOT  BIOMASS 

FLVTB        TABLE  FLV  VS.  DEVELOPMENT  STAGE 
FOV  FRACTION  OF  TIME  THAT  SKY  IS  OVERCAST 

FSH  FRACTION  OF  SHOOT  IN  TOTAL  PLANT  BIOMASS 

FSHTB        TABLE  FSH  VS.  DEVELOPMENT  STAGE 
FSO  FRACTION  OF  STORAGE  ORGANS  IN  SHOOT  BIOMASS 

FSOTB        TABLE  FSO  VS.  DEVELOPMENT  STAGE 
FST  FRACTION  OF  STEMS  IN  SHOOT  BIOMASS 

FSTTB        TABLE  FST  VS.  DEVELOPMENT  STAGE 

GLV  GROWTH  RATE  OF  THE  LEAVES  KG  (DM)/  HA/  D 

GPHOT        DAILY  GROSS  C02  ASSIMILATION  KG  (CH20)/  HA/  D 

GRT  GROWTH  RATE  OF  THE  ROOTS  KG  (DM)/  HA/  D 

GSH  GROWTH  RATE  OF  THE  SHOOT  KG  (DM)/  HA/  D 

GSO  GROWTH  RATE  OF  THE  STORAGE  ORGANS  KG  (DM)/  HA/  D 

GST  GROWTH  RATE  OF  THE  STEMS  KG  (DM)/  HA/  D 

GTW  GROWTH  RATE  OF  TOTAL  PLANT  BIOMASS  KG  (DM)/  HA/  D 

LAI  LEAF  AREA  INDEX  M2/M2 

LAT  LATITUDE  DEGREE 

MAINT        MAINTENANCE  RESPIRATION  OF  THE  VEGETATION  KG  (CH20)/  HA/  D 

MAINTS      MAINTENANCE  RESPIRATION  AT  STANDARD 

TEMPERATURE  (25oC)  KG  (CH20)/  HA/  D 
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NAME         DESCRIPTION  UNIT 


NWRT         NEGATIVE  WEIGHT  OF  ROOTS  (OUTPUT  VARIABLE)  KG  (DM)/  HA 

PI  CIRCUMFERENCE  OF  A  CIRCLE.  DIVIDED  BY 

ITS  DIAMETER 

QIO  INCREASE  IN  RATE  OF  MAINTENANCE  PROCESSES  PER 

10  DEGREES  C 

RDN  AVERAGE  LEVEL  INCOMING  PHOTOSYNTHETIC  ACTIVE 

RADIATION  J/  M2/  S 

RDR  RELATIVE  DEATH  RATE  OF  THE  LEAVES  1/0 

RDRTB        TABLE  RDR  VS.  DEVELOPMENT  STAGE 

REFLC        REFLECTION  COEFFICIENT  OF  THE  CANOPY  (FRACTION) 
SLFA         SPECIFIC  LEAF  AREA  HA  (LEAF)/  KG  (LEAF) 

TADRW       TOTAL  ABOVE-GROUND  BIOMASS  KG/HA 
TEFF         EFFECT  OF  TEMPERATURE  ON  RATE  OF  MAINTENANCE 
RESPIRATION 

TIME         SIMULATED  TIME  DAY 
TMAXT        TABLE  MAXIMUM  TEMPERATURE  VS.  DAYNUMBER 
TMINT        TABLE  MINIMUM  TEMPERATURE  VS.  DAYNUMBER 

TMPA         AVERAGE  AIR  TEMPERATURE  DEGREE  C 

WLV  WEIGHT  OF  THE  GREEN  LEAVES  KG/HA 

WLVD  WEIGHT  OF  THE  DEAD  LEAVES  KG/HA 
WLVI  INITIAL  WEIGHT  OF  THE  LEAVES  KG/HA 

WLVT         WEIGHT  OF  THE  GREEN  PLUS  DEAD  LEAVES  (OUTPUT 

VARIABLE)  KG/HA 
WRT  WEIGHT  OF  THE  ROOTS  KG/HA 

WRTI  INITIAL  WEIGHT  OF  THE  ROOTS  KG/HA 

WSO  WEIGHT  OF  THE  STORAGE  ORGANS  KG/HA 

WST  WEIGHT  OF  THE  STEMS  KG/HA 

WVEG         WEIGHT  OF  THE  VEGETATIVE  PARTS  (OUTPUT 

VARIABLE)  KG/HA 
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